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1 Collision Dynamics of Stretched Atoms 

The abstract of the original proposal awarded 6/99 is: “The theoretical study of collisions 
with stretched (excited and Rydberg) atoms is proposed. Classical and semiclassical expres¬ 
sions for the inelastic form factor will be developed in such a way that they will form the 
basis of new theories for the cross sections of (n£ —► n't') collision processes involving excited 
atoms. Stark £-mixing [nt -► nt') -will also be investigated by classical techniques. Collisions 
involving planetary atoms will be studied by adiabatic invajiioiice methods. 

This report will ftirnish the huge successes and acheivements accomplished under AFOSR 
funding. 

2 Importance 

Collision Dynamics of Stretched atoms is of key importance in several classes of problems, 
such as the following 

(A) Three-Body Recombination of positron with antiprotons p, 

e++p + e++ (1) 

to form antihydrogen H at positron (cryogenic) temperatures 4K. 

(B) Electron - ion recombination process 

e + Xe'^ + e-* Xe + e (2) 

at electron temperatures of Te ~ 5 mK. The theory to be developed for the above 
processes is important in that it iTiUSt be derived from fc st principles and will therefore 
provide a fundamental understanding of two quite different examples of recombination at 
ultracold energies. 

The theory and computational techniques under development for (A) and (B) are also 
important for understanding current experiments at NIST on recombination and for trapping 
experiments to be performed at CERN in the near future. 

2.1 Antihydrogen Formation 

At the Low Energy Antiproton Ring (LEAR) at CERN, antiprotons have been captured in 
specially designed Penning traps, cooled to meV energies and stored for hours in a small 
volume of space. Positrons have also been similarly accumulated. 

The current problem is how to arrange them to recombine efficiently into H and then 
controlling H in such a way that permits spectroscopic mesurements. Trapped atom tech¬ 
nology can then be applied to H. Since the sharp (2s —► Is) double photon transition in H 
can be measured to an accuracy of 1 in 10‘®, then, by comparing the measured values of the 
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Rydberg and anti-Rydberg wavenumbers, H becomes a fundamental physical system ideal 
for new tests of fundamental symmetries under combined operations of C (charge conjuga¬ 
tion) P (parity) and T (time reversal). CPT invariemce is minimed condition for the existence 
of antiparticles within quantum field theorj^ Moreover, by comparing the gravitational red 
shift of the (2s Is) transition for H and H, the first test of the gravitational Weak Equiv¬ 
alence Principle (WEP) for antimatter would then be possible. New basic theory on (1) at 
ultracold temperatures is required. 

2.2 Electron-Ion Recombination 

The theory under development for anti-hydrogen recombination (1) will also be applicable 
to the electron-^ ion recombination process 

e + Xe-^ + e^ Xe + e (3) 

and 

e -f Rb'^ + e Rb + e (4) 

at electron temperatures of Tg ~ 5 mK. Accurate experiments bn these systems are currently 
under investigation by S. Rolston at NIST and by P. Gould at University of Connecticut. The 
eventual comparison with our theory under development is of key significance and importance 
since it will be the first accurate comparison of ultracold recombination. As has been shown 
in our current investigation [1], the mechanisms of recombination at ultracold energies are 
quite different from those for recombination at room and higher temperatures. 

2.3 Basic Collision Problems 

Within the general theory of the three-body recombination processes above, there axe sev¬ 
eral key collision mechanisnos, operating at ultralow energies, which are -addressed in this 
proposal. On changing antimatter signs to matter signs, these sire: 

1. Theory of Stark Mixing collisions 

A+ -K H{n£) ^ A+ -t- Hint) (5) 

2. Theory of nC — n'i' ultracold collisional transitions 

P + H{nl) H{n'£') (6) 

where projectile P can be electrons or ions. These are important to (1), (3) and (4). 
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3 Present Results from AFOSR support 

3.1 XJltralow energy Three-body electron-ion recombination 

Three-body electron-ion recombination is described [1] at ultralow electron temperatures T*. 
At the initial stage involves extremely rapid collisional capture into high Rydberg states 
n > 100 - 400 with high angular momentum / w n - 1 at a rate ~ This is followed 

by extremely slow collisional-radiative decay. The key coilisional mechanism appears to be 
collisional /-mixing of the Rydberg atoms A(n) by ions and electrons until sufficiently low Vs 
are attained so as to permit relatively rapid radiative decay to the lowest electronic levels. 
This sequence is in direct contrast to the sequence of much slower collisional capture at 
higher T* followed by the much faster decay oi A(n) by electron collisions to lower levels 
where radiative decay completes the recombination. At ultra-low temperatures, the rate 
limiting sequence is therefore collisionad /-mixing followed by radiative decay in contrast to 
recombination at much higher energies and electron densities Ne ~ 10® cm"®, where the rate 
limiting step is the initial collisional or radiative capture at intermediate re(~ 1 eV) and 
higher T* (~ 10 eV), respectively. 

3.2 Sub-Projects 

We therefore had to investigate (a) dynamical properties of Rydberg atoms, (b) Stark Mixing 

jil _» jii' collisional transitions and (c) nL —*• n'V collisional transitions in Rydberg atoms 

in levels n « 100 - 400, for which fully quantal calculations were prohibitively difficult. 
Highlights of these investigations are: 

1. We defined, formulated and developed the concept of a classical inelastic form factor. 

2. By appeal to the dynamictd SO(4) symmetry of we were able to provide the 

first exact classical and quantal solutions for the Stark Mixing problem for the whole 
array of nl —► nfl' transitions in Rydberg atoms. This has been a long out.®i anding 
problem in Atomic Physics for the past 40 years. 

3. By exploiting the classical form factor, we were able to formulate the first classical 
impulse theory for nl -* n'V collisional transitions in Rydberg atoms. 

All of this work resulted in the following seven publications in Peer Reviewed Scientific 
Journails 


4 




i 


4 List of Publications 

1. Classical Atomic Form Factor, by D. Vrinceanu and M. R. Flannery, Physical Review 
Letters, 82, (1999), pp. 3412-3415. 

2. Classical and Quantal Atomic Form Factors for arbitrary transitions, by D. Winceanu 
and M. R. Flannery, Phys. Rev. A 60, (1999), pp. 1053-1069. 

3. Quantal-Classical Correspondence Impulse Theory, by M. R. Flannery and D. Vrinceanu 
Physical Review Letters, 85, (2000), pp. 1-5. 

4. Quantal Stark Mixing at ultralow energies, by D. Vrinceanu and M. R. Flannery, J. 
Phys. B^: At. Mol. Opt. Phys. 33 (2000), pp. L721-8. 

5. Classical Stark Mixing at ultralow energies, by D. Vrinceanu and M. R. Flannery, 
Physical Review Letters, 85, (2000), pp. 4880-3. 

6. Analytical quantal collisional Stark Mixing Probabilities, by D. Vrinceanu and M. R. 
Flannery, J. Phys. B: At. Mol. Opt. Phys. 34 (2001), pp. L1-L8. 

7. Classical and Quantal Collisional Stark Mixing at ultralow collision energies, by D. 
Vrinceanu and M. R. Flannery, Phy. Rev. A 63 (2001) 032701-(l-22). 

4.1 Referee Reports 

1. The referee’s report on paper ^^1, Classical Atomic Form Factor, Phys. Rev. Letts, 
82, (1999), 3412 reads; 

' “.... This is a very nice example of how classical mechanics can be used to understand 

the behavior of quantum mechanical observables, even for comparatively small quan¬ 
tum numbers. The present work is new and simple. It is important because it makes 
a significant contribution to advancing our understanding of the relation between clas¬ 
sical and quantum mechanics. This is a subject of substantial and growing interest in 
almost all fields of modem physics, so there is no doubt that the criterion of ’’broad 
interest” is fufiUed. The manuscript is concisely and well written; it should be easily 
understandable by a broad audience.I do not hesitate to recommend publication in 
Physical Review Letters.” 

2. The referee’s report on Paper #2, Classical and Quantal Atomic Form Factors for 
arbitrary transitions, Phys. Rev. A 60, (1999), 1053 reads: 

“This is a remarkable paper with a great deal of physical insight combined with sys¬ 
tematic analyses of the numerical results. The unified classical/quantal treatment is 
both interesting and useful. I reconunend publication in the Physical Review A.” 
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3. In the referee report on paper #3 Quantal-Classical Correspondence Impulse Theory 
Phys. Rev. Letts, 85, (2000).!. the referee wrote: 


“ In this paper a new way of formulating the quantal impulse approximation is pre¬ 
sented which allows one to obtain its classical limit in a natural and rather simple way. 
The presented approach which involves the VVigner functions of initial and final states 
exhibits in a lucid way the quantitative relationship between this quantum impulse 
approximation and the binary encounter approach. So far this connecti<^n has been 
rather obscure. As far as inelastic collisions are concerned, the presented results have 
a great potential in inspiring new semiclassical approaches. 

In my t5pinion,the main result of this paper is of topical significance in two respects. 
Firstly, the quantal impulse approximation is of central importance in general scatter¬ 
ing theory. Secondly, recently there has been a general strong interest in the power of 
classical dynamics and in its connection with quantum dynamics (Compare with Refs. 
10-15). As a minor remark, in this context it might be worth not only to refer to time- 
independent problems but also to some semiclassical work dealing with explicitly time 
dependent problems (e.g. G. Alber and 0. Zobay, PRA 59, R3174 (1999)). Puther- 
more, the presentation of the material is clear and scientifically sound. I recommend 
publication.” 

4. In a referee report on paper #4, “Quantal Stark Mixing at ultralow energies”, J. Phys. 
B: Atom. Mol. Opt. Phys. 33 (2000) L721, the referee wrote: 

“The article presents an elegant derivation, based on group theoretical arguments, of 
angular mixing transition probabilities between £ and t states in hydrogen induced by 
low energy collisions with charged particles...” 

5. In a referee report on paper ^5 “Classical Stark mixing at ultralow collision energies , 

Phys. Rev. Letts, 85, (2000), 4880, the referee wrote: _ 

“The paper presents an elegant solution to a long-standing problem in atomic physics, 
but the method of approach is of interest to a much wider community because the 
authors show how to exploit the classical/quantum correspondence in a system with a 
rich dynamical symmetry. It is rare to find so comprehensive a solution to a difficult 
problem....” 

6. In a referee report on paper ?^6, Analytical quantal collisional Stark Mixing ProbabUi- 
ties, J. Phys. B: At. Mol. Opt. Phys. 34 (2001), LI, the referee writes: 

“The authors present interesting analytical and numerical results on collision-induced 
Stark mbdng probabilities for hydrogen and, more generally, for Rydberg atoms. These 
results are more compact and easier to use than the authors corresponding findings in 
previous publications..” 
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7. In two referee reports on our paper 4^7, “Classical and Quantal Collisional Stark Mixing 
at low energies”, to be published in Phys. Rev. A 63 (2001) 032701-(l-22), the referees 
wrote: 

(A) ’’The paper under consideration makes an important advancement in the theory of 
collisional 1-mixing. These new results are given by Eqs. (28) and further by Eq. (34). 
The authors succeeded in presenting the formulae in rather compact form reflecting 
mathematical beauty of the problem, in addition to its pragmatic value. This material 
and analysis of the results certainly deserve publication.” 

(B) It is quite clear that this manuscript meets the requirements for publication in 
Physical Review A. The authors have presented a uniform approach to the problem of 
collisional Stark mixing in slow collisions between Rydberg-atoms and charged parti¬ 
cles, and have obtained analytic solutions, both classical and quantal, for the full array 
of intrashell nl-nl’ transitions. Previous work on this problem was restricted to initial 
s-states; the present work achieves a significant generalization for arbitrary initial val¬ 
ues of the angular momentum. This manuscript could well have formed the basis of 
three separate papers, one on the classical solution, one on the quantum mechamcal 
treatment and a third on the Monte Carlo simulations; However, by using the SO(4) 
isomorphism to exploit the rich symmetry of the problem, the authors have achieved 
a uniform approach that points out important insights and connections between the 
classical and quantum results. 

The present uniform approach evidently allows an elegant and complete solution to 
the problem that had not been previously achieved. 

5 Classical and Quantal Form Factors 

THj first breakthrough came from our definition of a “classical” form factor, in classical 
cf'despondence to the quantal form factor 

The physical significance of the form factor is that Pi/(q) = is the probability [2] of an 
internal transition arising from any external impulsive perturbation (whether due to collision 
with particles or exposure to electromagnetic radiation) which induces a sudden change q 
to the internal momentum of the target system. This work resulted in two published papers 
(1) and [2] above. 

5.1 Publications; Classical and Quantal Form Factors 

1. Classical Atomic Form Factor, by 

D. Vrinceanu and M. R. Flannery, Phys. Rev. Letts, 82 (1999) pp. 3412-3415. 

Here, the general trends exhibited in the variation of the inelastic form factor in collisional 
transitions nl n'l', when I' is changed and n, I and n' are kept fixed were explained solely in 
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terms of classical mechanics. Previous quantal results were reproduced from purely classical 
mechanics principles. Our conclusions are valid not only for large quantum numbers (which 
provide the usual classical correspondence) but also for other cases, which, up to now have 
only been described by quantal or semiclassicai methods. The interesting trends exhibited in 
the form factor are directly reflected in experimental and theoretical treatments of collisions 
involving excited atoms. 

2. Classical and Quantal Atomic Form Factors for arbitrary transitions, by 
D. Vrinceaiiu and M. R. Flannery, Phys. Rev. A, 60 (1999) pp. 1053-1069. 

Here, the classical form factor was deduced from exact correspondence with a phase 
space representation of the quantal form factor. Analytical expressions were provided for 
nl -* n'l', nl^ n' and n n' transitions in hydrogenic systems and for n n'm the 
one dimensional harmonic osciUator. An eflicient procedure for calculation of quantal form 
factors as analytical functions of momentum transfer, for arbitrary quantum numbers, was 
presented. The classical approach has the ability to explain quite succinctly interesting 
trends and various important aspects which remain hidden within the quantal treatment 
of form factors. The classical/quantal comparison ranges from being qualitatively good for 
nl —► n'l' transitions to close agreement for rd -* n' and n -* n' transitions. Excellent 
agreement is obtained for the integrated form factor for all transitions. 

6 Collisional Stark Mixing in Rydberg Atoms 

In the original proposal awarded 6/99, we mentioned in the section on Theoretical Method 
that “ By appeal to scattering theory, theories for Stark mixing of Rydberg atoms initially 
in high Estates must first be developed. The mixing wiU occur by coUisions mainly with the 
heavy particles as in 

p-^H{n,l)-*p + H{n,l'). 

and by micro fields in plasma. The theory must account for large angular momenta changes 

A/ = 1 — i' and occurs from large £ « n — 1.” 

During the past 24 months (12/15/98-12/14/2000) of the grant, research was conducted 
by the PI (M. R. Flannery) and Daniel Vrinceanu, a Ph. D graduate student in the School 
of Physics, working under the Pi’s supervision as Thesis Advisor. It was discovered [l] 
that the first stage in recombination at ultralow temperatures Te is a very rapid collisional 
capture into high (circular) Rydberg states n£, with high angular momentum £ ~ n - 1, 
ft nH large radiative lifetimes, at a rate proportional to T7^ ®. Stark Mixing occurs when the 
electron of a Rydberg atom (in a state with principal quantum number n) changes its angular 
momentum, without changing its energy, as a result of a collision, at large impact parameter 
6 with a slow massive psirticle of charge Z\e moving with velocity v. The £-mixing collisions 
al-e essential in producing the low angular momentum states required to radiatively decay 
at relatively high rate to low n-levels, thereby stabilizing the recombination. This sequence 
is in direct contrast to the sequence of much slower collisional capture at higher T* foUowed 
by the much faster decay of A(n) by electron collisions to lower levels where radiative decay 
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completes the recombination. The £-mLxing was caused mainly by collisions with ions rather 
than by light electrons. A full theoretical treatment of nl —> nl' transitions by collisions with 
ions at ultralow energies was therefore initiated. 

6.1 Collisional Stark Mixing: Research Performed 

On considering the Ryaoerg atom in a frame rotating with the internuclear a.xis. the Stark 
Mixing problem can be reduced'to the problem of the Rydberg atom in mixed static fields: 
electric, provided by the projectile ion and magnetic, produced by the non-inertial (Coriolis) 
forces. In this way, the well known equations, in both classical [4] and quantum [5, 6] 
mechanics, fonthe problem of interaction between weak fields and an atom can be adopted 
to provide, in principle, a solution for the Stark Mbdng problem. Both quantal and classical 
[7] versions of this approach have succeeded only for ^ = 0 to higher angular momentum i' 
transitions, appropriate to the experiments[8]. 

Since three-body collisional capture produces the Rydbergs in high i states, a new theory 
for the full array ni — nt of collisional transitions in H{nJ) was required. The present 
new treatment [9] developed during the grant period is not an extension of any previous 
theory and is capable of providing the first comprehensive classical and quantal solutions 
for Stark Mbdng. Our research has resulted in formulating and developing a unified theory 
for the general time dependent solution of Collisional Stark Mbdng. The exceptional rich 
dynamic symmetry of the hydrogen atom provided the key foundation which enabled both 
exact classical [9, 10] and exact quantal [9, 11] solutions to be constructed in a unified way 
by using group representation theory. This classical-quantal correspondence transcends the 
well known Ehrenfest’s theorem just because of the SO(4) dynamical group symmetry of 
the energy shell of the hydrogen atom. New classical and quantal solutions applicable to 
transitions between arbitrary angular momentum states have been derived. A new expres¬ 
sion [10] for the classical transition probability Pvt is defined in a language which exploits 
the dynamical symmetry _The derived classical probability for the general array £ 
of transitions has a v^iy simple fimctional form, can be easily calculated for any principal 
quantum number and provides physical insight and simple geometrical explanations for the 
behavior of the transition probabilities. Monte Carlo simulations [9] were also performed to 
yield results in agreement with both the exact quanta! and classical probabilities. Represen¬ 
tative quantal, classical and Monte Carlo probabilities Pvt are shown in Figs. [1-5]. Many 
similar results are contained in Re£s.[9, 10, 11]. Here, the Stark parameter a is 3Za„Un/2u6 
and the classically allowed and forbidden regions correspond to positive and negative values, 
respectively, of the parameter defined in Ref. [9]. A step discontinuity in the clas¬ 

sical Pt't and an exponential behavior in the quantal Pfi occurs as i' increases through this 
region. Within the classical region, the quamtal Pvt oscillate about the classical background. 
Also, cusp-type singularities occur in the classical Pvt when the parameter A{i, t') of Ref. 
[9] passes through zero. By revealing essential characteristics which remain obscured within 
the quantal treatment, the classical results complement the quantal results. 

This work on Collisional Stark Mixing has resulted in four papers [9, 10, 11, 12]. 

Summaries are as follows: 
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Figure 1: The Monte-Carlo simulation (step-like lines), the classical (solid line) and quantal 
(dots) transition probabiUties Pet{oc) for a given Stark parameter a = 0.2 and initial £ = 12 
within the n — 20 shell. ^ < 0 for all £^ 



Figure 2: The Monte-Carlo simulation (step-Uke lines), the classical (solid line) and quantal 
(dots) transition probabilities Pee{a) for a given Stark parameter q = 0.6 and initial £ = 14 
within the n = 20 shell. B > 0 for all £'. 


10 










Figure 3: The Monte-Carlo simulation (step-like lines), the classical (solid line) and quantal 
(dots) transition probabiUties Pfi{a) for a given Stark parameter a = 0.6 and initial ^ = 18 
within the n = 20 shell. Across the first dotted line B ch^ges sign and A is negative in 
both sides of this line. A changes sign across the second line while B remain positive. 



Figure 4: Stark Mbdng probability for transitions within the energy shell n = 100 from the 
initial state i = 80, with angle tj given by COST 7/2 = 0.1, as a function of the final angular 
momentum Exact quantal result is represented by the dots, Classical limit is represented 
by the solid line and the initial state is indicated by the dotted line. 
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Figure 5: Stark Mixing probability for transitions within the energy shell n = 100 from the 
initial state C = 60, with angle v given by cos7?/2 = 0.9, as a function of the final angular 
momentum t. Exact quanta! result is represented by the dots, classical limit is represented 
by the solid line and the initial state is indicated by the dotted line. 

6.2 Publications: Collisional Stark Mbcing 

3. Quantal Stark mixing at ultralow collision energies, by 

D. Vrinceanu and M. R. Flannery, J. Phys. B 33, L721 (2000). 

Here, a new exact solution of the time*dependent quantal equation was obtained for the 
full array of angular momentum mixing transitions n£ —* n£ in atomic hydrogen induced by 
collisions with charged particles at ultralow energies. Based on this new solution, efficient 
numerical procedures were devised. It was proven that the present (fixed frame) solution is 
equivalent to the rotating frame approach described by Kr..zansky and Ostrovsky [6, 7) and 
that it overcomes the difficulties therein. Analytic expressions for low quantum numbers n 
were presented. Numerical results for the transition array with n — 28 were reported. 

4. Classical Stark Mixing at ultralow collision energies, by 

D. Vrinceanu and M. R. Flannery, Phys. Rev. Letts (85, 4880 (2000)) 

Here, exact solutions of the time-dependent classical equations were obtained for the full 
array of an gular momentum mixing transitions n£ —* ni* in atomic hydrogen induced by 
collisions with charged particles at ultralow energies. A novel classical expression for the 
transition probability Pet was presented. The exact classical results for Pee{ot) as a function 
of i, a and the Stark parameter a agreed exceptionally well with (exact) quantal results. 
They complement the quantal results by revealing essential characteristics which remain 
obscured in the quantal treatment. 

5. Analytical Quantal Collisional Stark Mixing probabilities by 

by D. Vrinceanu and M. R. Flannery, J. Phys. B: At. Mol. Opt. Phys. 34 (2001), pp. 
L1-L8. 
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Exact expression for the probability P/,"’ of Stark Mixing transitions between arbitrary 
angular momentum states i and i' within the same energy shell n, as a result of a collision 
with a slow charged projectile at large impact parameters is presented. The formula obtained 
is compact and easy to use for numerical evaluations even for very large quantum numbers 
(n ~ 100). A classical approximation is directly obtained and compared with the exact 
quantal result in the limit of large n. Two distinct sets of quantal oscillations are predicted. 

6. Classical and Quantal Collisional Stark Mixing at low energies, by 

D. Vrinceanu and M. R. Flannery, Phys. Rev. A 63, (2001, March issue) 

Here, exact classical and quantal solutions were presented for the full array of intrashell 
transitions n£ — nf, between any angular momentum states, induced by slow distant colli¬ 
sions with a clfarged particles. The collisions considered are adiabatic with respect to orbital 
frequency and sudden with respect to Stark precessional frequency. The rich symmetry of the 
problem allows a unified approach and is the source of the excellent agreement, beyond the 
usual Ehrenfest’s correspondence principle, between the classical and quantal treatments. 
A classical transition probability is defined. Probabilities for transition between any an¬ 
gular momentum states within a high Rydberg energy level are derived in exact analytm 
forms and are analyzed for a large number of numerical examples. The transition probabili¬ 
ties obtained from the three methods - quantal and classical formulations and Monte-Carlo 
classical simulations - are directly compared to provide excellent agreement. The common 
50(4) symmetry provides this classical-quantal correspondence at a level, more fundamental 
than Ehrenfest’s theorem and the Heisenberg correspondence. The classical method is abo 
complementary in that it reveals very succintly essential and valuable characteristics which 
remain obscured within the quantal treatment. This reflects the essential power of classical 
dynamics when applied to collision problems. 

6.3 Publication: nl n'l' Collisional Transitions 

A successful formulation [13] of a classical version of the quantal irup-*lse treatment was 
made and applied to nl — n'l' Collisional Transitions in Rydberg acorns. The summary is 
as follows: 

7. Quantal-Classical Correspondence Impulse Theory, by 

M. R. Flannery and D. Vrinceanu, Phys. Rev. Letts, 85 (2000) pp. 1-5. 

Here, the quantal impulse cross section was derived in a novel form appropriate for 
direct classical correspondence. The classical impulse cross section is then uniquely defined 
and yields the first classical expression for ni - n'i' collisional transitions. The derived cross 
sections satisfy the optical theorem and detailed balance. Direct connection with the classical 
binary encounter approximation is also firmly established. The unified method introduced is 
general in its application to various collision and recombination processes and enables new 
directions of enquiry to be pursued quite succinctly. 
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7 Invited Papers Presented at Scientific Conferences 

The following four papers were invited papers presented at the following conferences. 

1. The Recombination Era from the Big Bang to Modem Astrophysics, Aeronomy and 
the New State ( Bose-Einstein Condensate) of Matter, M. R. Flannery, Invited Dis¬ 
course,presented at The Royal Irish Academy, Dublin, Ireland, April 26. 1999. 

2. Dissociative Recombination, M. R. Flannery, Invited Paper, presented at the Fourth 
International Conference on Dissociative Recombination, Nasslingen, Stockholm, Swe¬ 
den, June 16-20, 1999. 

3. Three-^dy Collisional Ion-Ion Recombination, M. R. Flannery, Invited Paper, pre¬ 
sented at the Gordon Conference on Simple Systems in Chemistry and Physics, New- 
port, RI, July 11-16, 1999. 

In addition, various research papers were presented at: 

1. DAMOP in 1999 

2. 17th International Conference on Atomic Physics, Florence, June 4-9, 2000. 

3. Atoms, Molecules and Quantum Dots in Laser Fields: Fundamental Processes, Pisa, 
June 12-16, 2000. 

Also an Invited Colloquium: ’’Recombination”, was presented to the Department of 
Physics and Astronomy, University of Nebraska, Lincoln, September 30, 1999. 

8 Conclusion 

During the period (6/12/99-present) of our current AFOSR grant, we have been §i-ia:essful, 
as described herein, in providing (a) an exact classical and quanta! theory of collisional 
Stark Mbdng and (b) the definition and development of a classical form factor which we 
have shown is basic to our formulation of (c) a classical impulse theory oird — n't' collisional 
transitions. Much more work needs to be done along the lines proposed in my Renewal 
Proposal to AFOSR and which, is felt will lead to many more fundamental breakthroughs 
as those described in the present report. 
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The general trends exhibited in the variation of the inelastic form factor in collisionaJ transitions 
nl n'V, when T is changed and n, L and n' arc kept fixed, are explained solely in terms 
of classical mechanics. Previous quantal results are reproduced from purely classical mechanics 
principles. Our conclusions are valid not only for large quantum numbers (which provide the 
usual classical correspondence) but also for other cases, which, up to now have been described 
only by quantal or semiclassical methods. The interesting trends exhibited in the form factor are 
directly reflected in experimental and theoretical treatments of collisions involving excited atoms. 
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With the advent of new technology which facilitates 
the accurate measurement [1] of electron-excited atom 
collision cross section there has also been renewed interest 
in the theory [2] of collisions involving Rydberg atoms. 
Recent experiment [1], in particular, has confirmed that 
the cross section for the quadrupole 2^5 —>3^0 transition 
in e - He(2^5) collisions is much higher than that for the 
pure dipole 2^S 3^ P transition at low and intermediate 

energies, in accord with the theoretical predictions of 
Ref. [3] (Bom and multichannel eikonal approximations). 
Flannery and McCann [4] have noted that this unexpected 
behavior is only part of a more general systematic 
trend in that (a) the 2^5 — n^D collisional transitions 
are predominant over all other transitions to the same n 
value, even for transitions to the electronic continuum, 
and (b) there is a unique value of the final angular 
momentum P that is preferentially populated in nl — n'/' 
transitions (n' » n) in collisions between Rydberg atoms 
and e>..ctrons or atoms. 

i ri'^ origin of this general behavior was traced [4] to the 

-lation with /' of the quantum mechanical inelastic form 
factor 

jF/<(q) “ - <^/(p + q)l<^j(p)> 

( 1 ) 

for i(n,l) —» f(n',l') transitions between atomic states; 

are the wave functions in position space and 
d>,./(p) - J ilfij{r)exp{-ipr/n)dT, the 

wave functions in momentum space. 

When an instantaneous impulse applied at f * /o trans¬ 
fers momenmm q to an atomic electron, the exact solution 
of SchrOdinger’s equation under Hamiltonian 

//(p.r.f) — p^/2m - e^jr - r ' qfi(/ - fo) (2) 

is 

'l'(r,r) - [I + - l)»(t - /o)]d',/«(r). 

where 9 is the Heaviside step function. The probability 
for I “ In/) —• / ■ In'/') transitions from the (2/ + 1) 
initial sublevels is then 


PniMq) - - X Kn'/'m'k‘’'^/*|n/m)P 

m,m* 

(3) 

also deduced in [5]. The probability of any impulsive 
i -»/ transition, whether due to particle collisions or 
electromagnetic field, is therefore 

PiM) - IJ/,(q)l'. (4) 

which provides physical significance to the inelastic form 
factor, a fundamental property of the atom. For impulsive 
collisions between a particle 1 and a Rydberg electron 2 
bound to a core 3, the overall transition matrix element T 
decomposes as [5] 

Ifi(q)Tn{q). ( 5 ) 

where the matrix element for (1-2) free-free elastic 
scattering in the (1-2) center of mass, is a function only of 
q, as for Coulomb scattering In’" or for 

Bom’s approximation, Tn * / V'(^l 2 )exp(/q^/^)d^l 2 • 
The probability of transition in the target atom per 
each 0-2) impulsive encounter is P,/ — ITi/l^/irizP. 
agreement with (4), 

The cross section is obtained by the following integra¬ 
tion of the form factor (5) over momentum change, 

“ il^] ( 6 ) 

where k-,,/ are the initial and final wave numbers of rela¬ 
tive motion of the projectile-target system of reduced 
mass W and q — - k/| is the momentum change. 

The scattering amplitude for (1-2) collisions of reduced 
mass Mi 2 is f\i - {2Mn/^7rR})Tn. For (1-2) slow col¬ 
lisions with scattering length a, the Fermi interaction 
V(ri 2 ) - [47ra(«*/A/i2)]5{ri - r 2 ) also yields decom¬ 
position (5) with /i 2 “ o- 

The inelastic quantal form factor therefore not only 
exerts primary importance in collision studies, but also 
has a deep physical reality. In recent experimental studies 
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of excitation of Rydberg atoms by short unipolar half¬ 
cycle electromagnetic pulses the transition amplitude is 
determined directly by the inelastic form factor [6]. 

Analytical quantal [7,8] and semiclassical [9] form 
factors are available, although general systematic trends 
cannot be easily extracted from them. A classical form 
factor for n n' has been deduced [10] from binary 
encounter impulse theory and from a microcanonical 
distribution in en<*ri>y space. A key point of this paper 
is that a complementary classical approach for nl — n'V 
transitions can also be developed in a way which reveals, 
quite succinctly, important aspects which remain hidden 
within the quantal treatment. 

Consider a Rydberg atom in a stationary (n,/) state 
with energy E and Aguiar momentum L. If the atom 
is perturbed by any general impulsive field [as in Eq. (2) 
or the Fermi interaction], then the transition probability to 
the final state (n', /') (of energy £' and angular momentum 
U) is the inelastic form factor. 

The quantal probability density for finding the electron 
in the radial interval (r, r + dr) is 

pUr) - (7) 

where Rni is the hydrogenic radial wave function ex¬ 
pressed in terms of the generalized Laguerre polynomial. 

The phase space of a classical atom, with Hamil¬ 
tonian W(r,p) - p^/2m + V(r), angular momentum 
L(r.p) ” r X p, and period t„i * v„i in stationary 
state (n, /) is populated according to the microcanonical 
distribution [8,10] 

p‘,drdp = {hvnihS{H - £)5(|L| - £)} 


normalized to (2/ l).sutes in all of phase space. On 
integrating over the momentum space p and angular p^ 
f of the configuration space r, the classical distribution is 


p‘niir)dr 


where the radial speed is given by mr^l2 
£ _ v(r) - (/ + For the Kepler 

atom (t„/ “ 2irn^ a.u.) and p,/ (in a.u.) is 
1 r-) I + 

Pniir) 


TTfl dliU « 

.lFI - -L - ILLL^T 

L r »■* J 


I 1 


TT/I- 


Hr)- 


(9) 


The quantal (7) and classical (9) radial probability den¬ 
sities are illustrated in Fig. I. As in the textbook example 
of the harmonic oscillator, the classical distribution has 
singularities at the corresponding turning points given by 
the radii (in a.u.) 

/j^{l ± e} - n^{l ± [1 - (/ + l/2)Vn^]''^}- 

( 10 ) 


R- 



FIG. 1. Classical and quantal radial densities of probability 
of localization for the stationary state of the hydrogen atom 
[£ - -1/(2 X 20^) and! - 8]. 


The classical distribution is zero outside the accessible 
region, bounded by /?*. 

By using definition (1) the transition probability (3) can 
be converted to the new form 

/’i/(q) “ O-ttR? j Pn/(r.p)Pw-(r-P + q)drdp. 

( 11 ) 

where the quantal distributions in phase space are given 
by p’(r,p) - ( 2 irfi)~^/^tl'(r)exp(-ip • r//i)<^*(p). 
This form is now suitable for classical correspondence 
obtained by replacing densities p^ by the phase space 
distributions (8). The basic definition of the classical 
form factor is therefore given by (8) and (11). The 
physical significance is that the initial and final sutes 
correspond to definite regions in phase space, populated 
according to the microcanonical distribution (8), and 
that the transition probability is given, in a geometric 
sense, by the amount of overlap of these regions. In 
configuration space alone, the regions are spherical 
shells with inner and outer radii given by (10), 
the pericenter (£") and apocenter {R*) of the Kepler 
orbit 

Analytical expressions with explicit dependence on 
q for quantal and classical probabilities for nl —* n'l', 
nl n', n n transitions are developed in a separate 
paper [8]. Rather than examining the 1' variation of (11) 
for a given q, the key results are more readily deduced 
and are easily transparent by investigating the probability 
for all momentum transfers 

/ R,/(q)dq-'(2nH)^f^drpUr)Pn'r(r). 

( 12 ) , 

where !R is the overlapping region in configuration 
space defined by intersection of (R,~,R*) and (/?/,/?/) 
intervals. 
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Inserting p(r) = An(r)r‘ with (9) in (12) gives the 
classical form factor (CFF) 


ill' + 1) drlr 

•'Rm.n rM'rjir)' 


where = Tmx(R7,Rj) and R^^tx “ niinfR, ,Rj) 
define the bounds of the overlapping region 31. Different 
overlap situations are illustrated in Fig. 2 for a represen¬ 
tative case. The gray region is the accessible region for 
the initial state and the curves ’ possible final state tra¬ 
jectories. Transitions occur only when the final state tra¬ 
jectory penetrates the initia' :.ate accessible region. The 
longer time spent by the electron on the final state trajec¬ 
tory within the initial state accessible region, the bigger is 
the transition probability. 

As I increases from, zero to its maximum value for 
circular orbits. R~ increases from zero to n^, while R* 
decreases from Iri^ to the same value For final states 
n' > Jin, then “ R* for all values of Three 
regions of overlap are then apparent and are, respectively, 
accessed as /' is increased. 

Region /, RJ < Rf —Here the overlap region 31 ■ 
(Ri'.R*) is determined solely by the initial state and has 
spatial extent which remains constant as /' is varied from 
zero to some value l\ where Rf ” Ri ■ There is always 
an orientation of the final orbit which will then intersect 
the initial orbit, as exhibited in Fig. 2, for (n 3, / ” 2) 
and («' =• 8, /' = 0 - 2) orbits. The /' variation of (^3) 
is contained solely within the increasing integrand (r/) '. 

Region II, Rf < R/.—Here the overlap region 31 ■ 
(RJ, R*) includes the / pericenter and has spatial extent 
which decreases, as I' increases, eventually to zero when 
Rj = R *. In this region, the initial and final orbits can 
intersect each other, as for the (n' * 8, I' 4) orbit in 
Fig. 2. The /' variation of (13) results from variation 
of both the increasing lower iimit R/ and the increasing 
integrand (r/)"'. 


Region III RJ > R* ■ —Here the initial and final 
trajectories no longer intersect, since the pencenter of the 
final state is greater than the apocenter of the initial state. 
This region where (i — /) transitions do not occur, as 
illustrated by (/»' = 8, /' = 5,6.7) orbits in Fig. 2. is the 
classically inaccessible region. 

The boundaries between regions 1 and II and between 
regions II and III occur, respectively, at /' = l\ where 
Rj(n',l') = R7(n,l) and at /' = h where RJin'J') = 
R*{n, /). Thus /| and /j are giy;en by 

^/i.2 y) " ~ (14) 

where e is the eccentricity [!-(/ + l/2)*/n‘]' * of the 
initial orbit. 

Variation of the CFF (13). with final angular momen¬ 
tum /' is then determined both by the lower integration 
limit (which is a constant /?," in region I and m- 
creases as R/ in region II) and by the integrand (r/) 
Figure 3 illustrates the general pattern. As 1' is increased 
from 0 to /i (region I). the increase in CFF originates 
purely from the increasing integrand As /' is var¬ 

ied from /] to lit th* increasing integrand is offset by the 
decreasing range (R/.R,*) of integration (region II). For 
I 2 < I' < n - I, CFF is zero because transitions are not 
classically allowed in region IE. 

At /' — /i the trajectories touch only at their corre¬ 
sponding pericenters and CFF has a turning point sin¬ 
gularity characteristic of classical descriptions. The zero 
radial speed of the electron at the contact point of both 
initial and final orbits causes the infinite CFF (transition 
probability). 

As is evident from Figs. 3-5, the agreement between 
the classical and quantal results is excellent in region I. 
even for small quantum numbers. In region II, the quantal 
results oscillate about CFF.. Since classical motion is 
confined to a definite region, the dramatic fall for large 
/' is steeper than that for the quantal case where states 



no. 2. Various final state in' - 8. - 1-7) trajectories 

and the initial accessible region corresponding to (n - 3, 

/ - 2 ). 



no. 3 . Characteristic dependence of the inelastic form factor 
on the final angular momentum for fixed n (- 3). / (- 2), 
and n' (- 8). Classical calculations: solid line; quantal results: 

dots. 
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FIG. 4. Classical (solid line) and quanta! (dots) inelastic form 
factor for transitions from state (n * 35, / “ 30) to (n' “ 55, 
/' = 0 — 54) sutes. ^ 


have exponential tails within the classical inaccessible 
region III. As expected from correspondence principles, 
for the larger quantum numbers, the quantum form factor 
tends to CFF, even in the regions 11 and HI, as shown in 
Fig. 4. For quasielastic transitions nl —► nl' the classical 
and quanta! results are in excellent agreement for all 
angular momenta (Fig. 5). The quantal results exhibit 
maxima in the neighborhood of /' * l\,l 2 where CFF has 
the classical singularities. The position of l\ defined by 
(14) in the limit of large /, where the eccentricity e — 0, 
is /,(/ - n - 1) - /iy2[l - l/2(n/n')^]‘/* - 1/2. an 
exquisite result for initial circular orbits. For n' ^ n, 
1] tends from the bottom to /i(/ —* n — l,n' » 
rty/l - 1/2, a key result in detailed agreement with 
that previously derived from consideration of the quantal 
momentum-space overlap (4). 

For small initial angular momentum /, < ^ 1 and 
/i is then zero so that the maximum CFF is-given by 
l^{l — 0) - 2nn - (n//i')^]‘'^ - 1/2. appropriate 
to highly eccentric initial orbits. In the n' » « limit 
then /:(/ 0. n' » n) - 2n - 1/2. As the initi^ I 

increases, there is therefore a slow variation (2n —• v2n) 
in the position h of the maximum of CFF, which is 



FIG. 5. Classical (solid line) and quantal (dots) inelastic form 
factor for quasielastic transitions from n “ 20, / " 10 state. 


pushed slightly to lower values. This theoretical predic¬ 
tion is also confirmed by the quantal results [4]. 

When the energy £' of the final orbit is not sufficient to 
accommodate the value of I 2 deduced above {n' < -Jin), 
the peak in CFF (as in Fig. 5) is given by li, provided the 
initial / is large enough. When the final n' is sufficiently 
small so that the lower l\ cannot be accommodated, 
i.e., li > n' - 1, CFF and the quantal result exhibit a 
monotonic increase confined to region I, which is always 
characterized by excellent agreement between quantal and 
classical results. 

In summary, the pattern exhibited by the /' variations 
(Figs. 3-5) is essentially identical with the quantal pat¬ 
tern. The positions of maxima of the /' variation of CFF 
depend strongly on the initial n and only weakly on the 
initial /. in agreement with the quantal calculations [4], 
which were restricted to certain cases. Excellent quan¬ 
titative agreement between classical and quantal results 
makes the classical form factor a very useful tool particu¬ 
larly at large quantum numbers (Rydberg atoms) where 
exact quantal results are not easy to obtain (either analyti¬ 
cally or numerically) and to use. due to the highly oscilla¬ 
tory nature of the wave function. Although the emphasis 
here is on the electron form factors, the present analysis 
is applicable also to form factors for transitions between 
rovibrational states of molecules. 

This work was supported by AFOSR Grant 
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The classical form factor is deduced from exact correspondence with a phase-space representation of the 
quantal form factor. Analytical expressions are provided for nl—n'V. nl-*n' . and n—»n transitions in 
hydrogenic systems and for n—n' in the one-dimensimal harmonic oscillator. An efficient procedure for 
calculation of quantal form factors as analytical funcuons of momentum transfer, for arbitrary quantum num¬ 
bers. is presented. The classical approach has the ;<h;!i.y to expldn quite succinctly interesting trends and 
various important aspects which remain hidden within the quantal treatment of form factors. The classical- 
quanta! comparison ranges from being qualitatively good for nl—m'l' transitions to close agreement for nl 
—n' and n-*n' transitions. Excellent agreement is obtained for the integrated form factor for all transitions. 
[51050-2947(99)09107-6] 

PACS number(s): 32.80.Cy. 34.50.-s. 3l.l5.Gy 


I. INTRODUCTION 


example, the generalized oscillator strength / is written in 
terms of the inelastic form factor as 


The inelastic form factor 




is a very basic quantity. It can accurately describe the overall 
response and dynamics of an atom or molecule involved in 
various processes or external interactions. It is also common 
in various schemes of approximation. 

The study of the hydrogen atom form factor serves as a 
pivotal starting point for the general study of electronic tran¬ 
sitions between highly excited states of atoms and molecules. 
The additional effects of the nonhydrogenic core may be 
incorporated via use of quantum defect theory. Inelastic tran¬ 
sitions between the ro-vibrational states of molecules can be 
snidied by appeal to the inelastic form factor for the har¬ 
monic oscillator. 

Inelastic scattering of incident (neutral or charged) par¬ 
ticles or of photons (or short bursts of electromagnetic radia¬ 
tion [1]) by a structured target can be decomposed into an 
internal structure part, provided by the form factor of the 
target, and a dynamic part which depends on details of the 
external projectile-target interaction. Bound-bound, bound- 
continuum . continuum-continuum, and ionization transitions 
are treated on the same footing by using the form factor. 

The quantal impulse, scmiclassical impact parameter, first 
Bom approximation, and binary encounter methods of colli¬ 
sion theory [2] focus on the transition probability as a dy¬ 
namic response of the target in the field of the projectile. The 
physical significance of the form factor T is that 
= 1.^^ is the probability of an internal transition arising from 
any external impulsive perturbation (whether due to collision 
with particles or exposure to electromagnetic radiation) 
which induces a sudden change to the internal momentum 
of the target system [3]. Since any interaction can be decom¬ 
posed as a scries of sudden interactions, the scattering cross 
section or other observables arc determined by integrating 
the form factor, multiplied by some weighting factor charac¬ 
teristic of the interaction, over momentum transfer q. For 


/= 






where K=q(ao/^) dimensionless, and where is the 
change in energy between the initial and final states. The first 
Bom approximation for the inelastic scattering of structure¬ 
less ions of charge Ze, of speed v, by a hydrogenlike ion of 
charge Z'e, is 


(T( v ) = 


SttZZ Qq 
(v/Vo)^ 


JX™. 


-which can emphasize small momentum transfers K occurring 
at high energies. The inelastic scattering of an ultraslow neu¬ 
tral "particle by a Rydberg atom is ^ 

( vIVqYJK^ 


where a is the scattering length of the projectile-Rydbcrg- 
electron interaction and VQ = e^lh is the atomic unit for ve¬ 
locity. This expression emphasizes intermediate and larger K 
occurring within the interaction distance a. 

Classical mechanics provides a good quantitative descrip¬ 
tion of excited states via the correspondence principles [4]. 
Classical mechanics also promotes physical insight into the 
process by transparent causality, and provides scaling laws 
and elucidation of the dynamics. A phase-space description 
combined with statistical properties (microcanonical distri¬ 
bution in most cases) are the basis for an alternative or 
complementary view of quantal phenomena. 

Based on the recognition of the above fundamental as¬ 
pects of the form factor, this paper presents results for quan¬ 
tal form factors, and defines the classical form factor for the 
highly excited hydrogen atom and harmonic oscillator. Clas¬ 
sical mechanics is advantageous here (a) in revealing essen¬ 
tial details of the dynamics for inelastic transitions, (b) in 
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explaining the interesting trends in the behavior of the form 
factor (as a function of the various variables), and (c) in 
predicting quantitative results inaccessible to rigorous quan- 
tal calculations because of the formidable numerical restric¬ 
tions imposed by the highly oscillatory wave functions. 

Some analytical quanii [5,6], semiclassical [7.8] and 
classical [9.6] form factors are available, but general system¬ 
atic trends cannot be easily extracted from them. A key point 
of this paper is that a complementary classic'll approach for 
general nl^n'V transitions is developed in such a way that 
reveals quite succinctly important aspects svhich remain hid¬ 
den within the quantal treatment The consistency of this 
approach is verified by applying it to nl^n* and n^n' 
transitions. The known results [6] arc then rederived in a 
unified way. 

The simple example of the one-dimensional harmonic os¬ 
cillator is treated in Sec. II. The correspondence between the 
quantal and classical form factor for inelastic transitions in 
this system is apparent. On this ground, a generalization for 
three-dimensional systems (like the hydrogen atom) becomes 
feasible. A definition of the classical form factor interpreted 
as a transition probability between two states described by 
microcanonical distributions in the phase space, is proposed 
in Sec. III. The classical form factor for n/—►n'/' is intro¬ 
duced in Sec. IV. and compared with the quantal counterpart 
Various summations, over the momentum transfer and initial 
and final quantum numbers are obtained within the same 
theoretical framework. The classical calculation are based on 
the microcanonical distributions presented in Appendix A, In 
Appendix B the classical calculations are detailed. An effi¬ 
cient algorithm for calculation of the various quantal form 
factors is introduced in Appendix C. 


This function provides an easy way to compute the form 
factor by using 


1 

d"*”'Ht.z\Q) 

vn!/i'! 

rtr'z 


and reveals the following analytical structure of the quantal 
form factor for the harmonic oscillator: 

where is a polynomial of order n + n'. 

The square of the absolute value of the form factor can be 
interpreted as the transition probability of the harmonic os¬ 
cillator when an impulsive interaction imparts momentum q 
[3], Based on this observation, a classical analog of the form 
factor can be defined. Consider the phase space to be popu¬ 
lated according to the microcanonical distribution. The tran¬ 
sition probability is then given, in a geometric sense, by the 
volume of that region in phase space where both initial and 
final states can coexist 

The density of probability in phase space for a given state 
of the harmonic oscillator is 

Pii()“("“•“ i) ^ ^)» 

where the normalization factor A^=a>/2ir ensures one par¬ 
ticle in all of the phase space. The transition probability is 
then given by the conditional probability that the the system 
in the initial state n is in the phase-space volume element dr 
dp, and the ^splaced final state has quantum number n' 
within the same volume element Hence 


IL FORM FACTOR FOR THE HARMONIC OSCILLATOR 

The simple example of the one-dimensiopui harmonic os¬ 
cillator with the Hamiltonian 


hw f 

drdpS{H{r.p)-(.n + {)hu) 
X Sim r,p + q)-in' + j)hu) 


H—p^f2m+mm^r^f2 

is considered in this section. The quantized energy levels are 
£, = (/j+l/2)fta», and the corresponding wave functions 
M,(r) have the generating function 


Phase space integration yields 

-[-C* + 2(n + n' + l)Q*-(n-n')^]"''*. 

TT 

( 2 ) 


f(t.r) 


2 

.“O ijn\ 


_-l/4 -w 
IT Fq 


Xexp[-(r/ro)*/2+ ^rr/ro-f^/2]. 


where ro= y/Khrua is the natural length of the oscillator. 

The inelastic form factor for the transition from level n to 
level n' when the momentum q is transferred is defined as 


and has the generating function, in terms of the dimension¬ 
less variable Q-qr^lyjlh 


where the dimensionless parameter Q^-{q^l2m)lhu) is, 
again, q^r\l7ti^. Then 

P,AQ)^ 

which shows that q is restricted to those values for which the 
square root is a re^ number, otherwise the transition is clas¬ 
sically forbidden, and is zero. The limiting values of 
Q = Q. are given by 

el=(n-l-n' + l)±>/(2n-t-l)(2n' + l). 


l{t,Z'.Q)— 2 ^ni>' I—. -T 

/i.m»0 yjnlmi 


=exp(-e^/2)exp[rz+j(»+2)C]. 


where the probability P„„‘ exhibits characteristic classical 
singularities. The classical transition probability (2) satisfies 
detailed balance. The transition probability for elastic colli¬ 
sion in-n') is 
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no. 1. Classical and quantal transition probabilities for the (n 
= 10)—►(/:' = 20) transition in the harmonic oscillator vs momcn* 
turn transfer Q. 




(3) 


The advantage of using the classical transition probability 
is illustrated by the following example. Consider the 10 
—►20 transition. The quantal expression deduced from Eq. 
( 1 ). 

-^10^2o(2)” 


g-e'/iQio 


=(-670442572800 
5225472000^323323 

+ 609493248000 - 228559968000 Q* 

+ 46884096000 Q* - 5860512000 Q* 


1.5 


0.5 


ii 

■i 




10 


15 


20 


FIG. 2. Classical (solid line) and quantal (dots) transition prob¬ 
abilities for any momentum transfer q from the state n = 4 to states 
n'=0-20. 


which reduces to 
2 


vQ. 


I IT -s\ rn I [s 

^4 — 


where s-\-Q\fQ\ 


in terms of the incomplete elliptic function F. 

Figure 2 shows that there is an excellent agreement be¬ 
tween the quantal and classical transition probabilities for 
general transitions in the harmonic oscillator for any 

momentum transfer. The singularity at n' = n arises from Eq. 
(3). The characteristic agreement is displayed in Fig. 2. 

in. PHASE-SPACE EXPRESSION FOR THE FORM 
FACTOR 


+ 468840960 8800 Q 

+ 820800e‘^“ I7l00q^^+200C‘®-Q^°) 

is rather large, and is numerically inefficient due to the os¬ 
cillations in the wave functions. The classical transition 
probability (2), however, has the simpler form 

Pio.2oiQ)=^/ V-Q‘‘+62C?"-l(X) 

within the classically allowed range 31“ V861^C^<31 
+ v861. 

The classical and quantal are compared in 

Fig. 1. The comparison exhibits excellent “background** 
agreement within the classically allowed region of Q, the 
characteristic classical singularities at Q*, and the charac¬ 
teristic exponential quantal tails in the forbidden region. The 
number of quanta! oscillations in is given by 

min(n,n*)+l which occur within the extent [(2/i-F l)(2n' 
-H 1 centered on the median value {n + n* + 1) of the 
range. 

The classical transition probability for any momentum 
transfer (or the integrated form factor) is 

P,f= Q)^C 


The quantal amplitude 

is expressed as the abo^ vintegrations over either configura¬ 
tion space or over momentum space, where the momentum 
wave functions are defined as 

4)(p) = (2irt)"^| 'I'(r)exp(-iq'r/ft)<ir. 

The transition probability 

can therefore be expressed as 

(q'/r)|<r“'"*|'I',(r)),(<l>/p+q)|<I>/(p))p 

= J <frdp['I^;(r)tf'<''*’>''*<l»/p+q)] 
x[q'i(r)<-‘'-''*<I>r(p)]. 

The quantal phase-space distribution may be defined by 
p( r, p)=( 2 irft )■ “'I'(r)e ■'P''*<!)•( p) 


— )(Q^ probability densities p(r) in configuration space 

ir Jq. ^ ■ * ’ and p(p) in momentum space (obtained by integrating the 
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quantal phase-space distribution over momentum or configu¬ 
ration space)*yield p{r)='/p(r,p)dp=|<^(r)|^ and p(p) 

= |4>(p)|^, respectively. This distribution is the standard or¬ 
dered version [10] of the Wigner distribution [11]. 

The transition probability is therefore 

/^.i<q) = (27rfi)^ I <irifpp,(r,p)p/(r.p+q). (6) 

This expression for the transition probability is now in a 
form appropriate for classical correspondence, obtained by 
replacing the quantal densities Pij (even though pi f have no 
direct physical interpretation) by the classical phase-space 
distributions p‘^(r.p). Thus 

is the basic expression for the classical probability for itnpul- 
sive transitions. The number of initial states in the phase- 
space element dvdp is pfdtdp^ and {2Trfi)^pf is the prob¬ 
ability that the final state is in the same phase-space element 
Two fundamental properties, corresponding to similar 
properties of the quantal result (5), can be readily proven for 
the classical transition probability (7). The classical distribu¬ 
tions satisfy 

p(r.p) = 2 Pn(r.p)-(2irft)~^, 

n 

which means that the total number states in the phase volume 
clement is drdp/(2irfi)^, the number of elementary phase- 
space cells. The probability of transition from initial state i to 
all states / is then 

2 Pi(<q)= J Pi{r,p)dTdp=gi, 

the statistical weight g, of the initial state. The second prop¬ 
erty provides the transition probability for all momenmm 
transfers q. Integration of Eq. (7) over all possible values of 
the momennim transfer q involves 

J p(r,p-i-q)dq-|p(r.p')dp'J 5(p+q-p')dq 

= p(r). 

where p(r) is the classical distribution in configuration 
space. Then 

I P,Kq)dq=(2irfi)^| p,(r)p/(r)t/r. (8) 

This classical property is again in correspondence with the 
quantal result 

j l’I'i(r)|*|'I'/r)|*dr. (9) 

which can be derived from Eq. (4). 

In the action-angle representation for bound states, the 
classical distribution is 


p^{J,w) dJdw== S(Jfh-n) dJdw/(2Trhr 

for a general £)-dimensional system with a set of action- 
angle variables (y.w) in a state specified by the set of quan¬ 
tum numbers n. The classical probability for n —► n’ tran¬ 
sitions is therefore 

P,f(q) = (27r^)“^ J J S(J/h-n) S(r/h-nl 

xS(p’^q-p') S{r-T') d^dw dJ' dw\ (10) 

which provides a more general classical correspondence with 
the quantal expression (6), rewritten as 

/»iKq)=(2irfi)^ J J drdpdr'dp’p,{r.p)p/<r'.p') 
X5(r-r')5(p+q-p'). 

These expressions emphasize the impulsive namre of the 
momentum transferred. 

IV. FORM FACTOR FOR THE HYDROGENIC ATOM 
A. Form factor for nl-^n'V transitions 

The classical distribution for an atom with given energy £ 
and angular momentum L in (r.p) phase space is (Appendix 
A) 

p(E,L-,T,p)dEdLdTdp 

drdp , ^ 

=‘dES{H-E„,)dLS(\U-L)^^^^, ( 11 ) 

where both the Hamiltonian ff(r,p)=p^/2m +V(r) and an¬ 
gular momennim |L(r.p)| = rpsintf^ are constants of mo¬ 
tion. The classical transition probability for the case of a 
central field one-electron atom, between the states, or bands 
of states centered at (E,L) and (E',L')» momentum 

transfer q, is then 

P(E,L;E'.L';q)dEdLdE'dL' 

=dEdLdE'dL' f d^^S{H(t,p)-E) 

J {2irfiy 

X 5(1 L( r. p) I - L) 5(H( r, p +q) - £) 
x5(|L(r.p+q)|-Z.')- (12) 

The quantity /»(r;r':q) is the transition probability density 
(per unit intervals dV dt '). When E, E', L, ot L are quan¬ 
tized. the transition probability between corresponding states 
arc obtained by the formal replacements dE-*h v„, and dL 
-*k on the right-hand side of Eq. (12). The transition prob¬ 
ability between bound states with given quantum numbers 
(rt.O and(n',/') is then 
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xP{E„,Al+m)h;E„.r .W + iri)h). 

(13) 

where r"/ is the radial frequency of the classical orbit. 

Since the densities used in Eq. (13) are already normal¬ 
ized to (2/ + 1) particles in all of the phase space (see Ap¬ 
pendix A), Eq. (13) represents the basic definition of the 
classical form factor, in direct correspondence with the (sym¬ 
metrical) quanial form factor 

F„,„-,'(q) = S 2 \{nlm\e‘<^-^'*\n'rm')\\ (14) 



0.06 0.08 0.10 0.12 


The physical significance of the basic expression (7) is that 
the initial and final states correspond to definite regions in 
phase space populated according to the microcanonical dis¬ 
tributions (11). Transitions can only occur if these two re¬ 
gions overlap, and the amount of overlap is a measure of the 
transition probability. The classical form factor (13) which 
has been developed in detail in Appendix B, will be directly 
compared with the quantal result of Eq. (14), developed in 
Appendix C as a function of q for arbitrary nl^n*I' transi¬ 
tions. 

The result of the classical calculation (13) (sec Appendix 
B for details) is 


’2(2 /+l){2/' + l)ft^ 


X 


f drlr^ 

G*,{r,q) + C~f{r,q) 

Jk rr' 

R 


Qir.q). 


(15) 


where ©(r,(y) is the step function having value unity within 
r<r* (given in Appendix B). and zero otherwise, and the 
function 





(16) 


must be real, so that q must be within the classically acces¬ 
sible range 11 given by 

= m^(r-r')^ + (Z. + Z.')^/r^ 

The radial velocity r (r') is a fiinctioa only of r for a given 
nl (or n'/') state. For a given momentum transfer q, the r 
integration proceeds over the radial ranges within which the 
square roots in Eq. (16) are real. This siniation is illustrated 
in Fig. 3 for the (4,3)-»(8,2) transition in the hydrogen 


FIG. 3. Integration region for a typical (4,3)—►(8.2) transition. 

atom. The dotted curve is the boundary (r^-r-O) of the 
region within which the function © is unity, and which en¬ 
compasses the physical accessible region 
When q is smair(below the A- curve) or large (above the 
B 4 , curve), Gjf has complex values and the transition prob¬ 
ability is necessarily zero. In the shaded regions only G J (a) 
or G J (b) or both G^f (c) can contribute to the integral for a 
given q. The range 71 of radial integration always lies within 
the region specified by A + =A-= real and = 

= real. The boundaries of this region are then given by 
/?^„=max(/fr Jf/) and R^=min(/?* Jl}), where and 
are the pcricenter and the apocenter of the Kepler orbit. 
The three situations possible (details in Ref. [3]) for the over¬ 
lap of the initial and final orbits are illustrated in Fig, 4 as L' 
of the final orbit is increased. Region I gives the maximum 
overlap when region 11-(R^ ,R^) is specified only by the 
initial state. In Region II the overlap is partial, because Jie 
lower limit of U is given by the pcricenter of the final .t 
RJ . In Region III the pcricenter RJ has moved outside R^ , 
so that there is no overlap, the transition is classically for¬ 
bidden. 

The quantal transition probability for bound-bound transi¬ 
tions, for hydrogcnic atoms, is a rational function in the mo¬ 
mentum transfer q. The proof of this statement and the algo¬ 
rithm for quantal calculations are presented in Appendix C. 
The form factor for the (4,3)-+(8,2) transition, as a function 
of q, is 

(I) totMl overlap (17) partial overlap (Ul) no overlap 



I' 


no. 4. Overlap situations for fixed initial n, /, and n' and 
varying / 
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34359738368<7 j977QQg755367_g4352288228992^*- 52980973020131328<7'‘ + 914957411105636352^?* 

{9 + (Aq-)^* 

-(-557017745580707807232<?* + 491446003263768782438409'°+ 1304678925402985186983936^'^ 

- 12446991865892540818391040<7 '* + 839683843672479677596827648^ 

+ 5269768456130999660417384448^'*+226149481324737121139390676992(7“ 

-3490013454414646748315148877824(7“+2'2442701774022980630594896003072^^* 

-634317556573974483528854336962569“+66458636923717615551358326276096^“) 


The results of calculations for the quantal and classical prob¬ 
abilities (form factors) are compared in Fig. 5. Tne four sin¬ 
gularities in the classical transition probability, wnich indi¬ 
cates maxima in quantal results, correspond with those 
values of the momentum transfer for which the (q 
= const) line is tangent to one of the curves A±=q or B 
i = q in Fig. 3. One of these equations has then a double 
root in r, which eventually yields after integration a logarith¬ 
mic singularity. 

In Fig. 6 the quantum numbers n,l, and n' are fixed and 
the transition probability versus the momentum transfer q is 
plotted for various final angular momenta I '. As I' increases, 
the quantal and classical transition probabilities increase, and 
attain a maximum for some value of V . This value is roughly 
given by in agreement with the results derived in 

Ref. [3]. Further increasing V produces a sha^ decreaw in 
the quantal transition probability. The classical transition 
probability is forbidden for I' —6 and 7, since there is no 
overlap at all between the phase-space regions occupied by 
the initial and the final states, for any momentum transfer q. 
This situation corresponds to region ED of F\g. 4. The quantal 
results for / ' = 6 and 7 are therefore classically inaccessible. 

When the final principal quantum number n' is varied, 
keeping n, I, and I' fixed, the shape of the transition prob¬ 
ability versus momentum transfer is preserved and the inag- 
nitude rapidly decreases as n' increases. This observation, 
valid for both quantal and classical cases, is demonstrated in 
Fig. 7 for a specific case. Because the transition probability 
(14) contains the factor I/t'- l/n'^ the classical form factor 
provides an explanation for this behavior. 



FIG. 5. Quantal and classical transition probabilities for the 
(4.3)—*(8.2) transition. 


When the final state is in the continuum with energy be¬ 
tween E' and E' + dE', result (15) is still valid provided 
hv„.,. is replaced by dE' . The probability for a bound-free 
transition is then 

f(2/+l)(2/' + l)fi^d£' 

P,,,,(E';q)t/£'=[--1 


f drtr^ 

C*f(r,q) + G~f(.r,q) 

JH r|r2 

9 


B. Integrated nl^n*V form factor 

The integrated form factor or the transition probability for 
any momentum transfer is the integral of the ^-dependent 
transition probability over the q space. The quantal calcula¬ 
tion gives 



no. 6. Quantal (a) and classical (b) transition probabilities for 
the (4.3) (84*) transitions, with / ' = 0.7. 
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= (2nh)^ f dq Z l'{',/«(r)|*2 

J m 

which, with 'i'(r) = {R^ilr)Y,„(,r). reduces to 

= 2ff-{2/+l){2/' + l)ft^| (17) 

where p'^{r)dr = Rli(r)dr is the radial probability. 

Integration of the corresponding classical transition prob¬ 
ability (8) gives 


/ 


PnU'l'(<l) <^q = - 


8-n-^{2/+I)(2/'4-l)ft^ r (dr/r^) 


nri 


/. 


n riTj 


(18) 



Upon integration, Eq. (18) yields 


no. 7. Quantal (a) and classical (b) transition probabilities for 
the (4,3) (n’,5) transitions, with n' = 6.12. 


g 

Pnin'i' = ~, —:F(arcsinVjll/i)[(-*3“-*2)(-*4~-*i)]''^ 

(X 3 -X 2 )(X 4 -X,) 

with 5 = 7 - -r, -r, 

(X3-X,)(X4-X2) 

v Vcfc F is the incomplete elliptic function and x, (i 
= 1,2,3, and 4) is the sorted set (/?," ,RJ ,Rp. When 
there is no overlap between the initial and final states (/?/* 
<RJ) the transition is of course classically forbidden (situ¬ 
ation III in Fig. 4). Comparison between the quantal and 
classical expressions reveals the definition of the classical 
radial probability: p^(r)dr^ldtlr, in agreement with the 
customary correspondence (deduced in Appendix A). The q 
integrated transition probabilities for fix^ initial quantum 
numbers n and / and final n* as function of the final angular 
momentum V are shown in Fig. 8. There is excellent agree¬ 
ment between the quanta! and classical calculations before 
the first singularity in which marks the transition from 
region I to region 11 in Fig. 4. For larger the quantal 
transition probabilities oscillate about the classical transition 
probabilities. As proven in Ref. [3], there is a limiting value 
of after which the quantal transition probability txpch 
ncniially decays while the classical form factor is zero. This 
situation corresponds with region in in Fig. 4, where the 
transition is classically forbidden. If this special value of 
cannot be accommodated, because — I, the transition 

is classically always permitted and the quantal transition 
probability has no exponential tail. This is the case of quasi 
elastic transitions (between the same principal quantum num¬ 
bers), as presented in Fig. 9. For this case, the agreement 


between quantal and classical calculations is excellent for 
any The q-integrated form factors were discussed in Ref. 

[3]. 

C. Form factor for ' transitions 

Summation over the final angular momentum number I* 
provides the form factor 

P.,.,>(q)= 2 2 (19) 

1*^0 m,m* 

The basic definition (7) gives the classical analog for this 
form factor 

/*,/.«' (q) = (2 IT* )* J dr dp p„( r. p)p,. (r, p+ q). 

where the densities p,, and /.averaged p,. are described in 
Appendix A. Using V{r)=‘E-p^l2m, the final distribution 
is rewritten in the r independent form as 

p,.(^,p+q) = (/j»',.)^(p+q)^/2m-p^/2m 

+ E-E'V{2irh)\ 

where t means “not r," The classical transition probability 
for impulsive nl—*n' transitions is then 
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RG. 8. Quantal (dots>^ classical (solid line) transition prob¬ 
abilities for the (35,30) (55,r) transitions, with /' = 0.54. 



no. 9. Quanta! (dots) and classical (solid line) transition prob¬ 
abilities for the (20,10) — (20,r) transitions, with r = 0.19. 


(summing) integrating the original probability Pn^r over 
all possible final angular momenta. The probability for nj 
—transitions in hydrogen is 


which is the classical overlap only of the momentum space 
distributions, rather than the full phase-space distributions, as 
in Eq. (12) for transitions. Since 

j S{pq cos tfp ,Im + q^l2m + E-E') dp-2irm/pq 

for p>po=\2m{E'-E)-q^\/2q, and zero otherwise, this 
transition probability reduces to 

2iTm f" / V * 

(*«'»') PMpdP' 
q jpo 

which involves or'y the momentum distribution of the initial 
sutc. The sa.he result is also obtained in Appendix B by 


2(21+1) , ._j 

f, /L(l+nV)\n"‘^ 


where all quantities are in atomic units, and L — 1/2. The 

limiB of integration are given by 

-\ln'^\l2q,p-]. where p. = (l *(1are the 
extreme values (at pericenter and apocenter) of the momen¬ 
tum of the electron on a given orbit. As a specific example, 
consider the 4,3—>8 transition. The quantal form factor is (cf. 
Appendix C) 


^^^^^^^^5^^^(9721215 + 4498478208 q^ + 471603326976 q* + 6554684489728 q® - 451062079160320 q* 
5(9 + 64q*)'’ 

+ 6344423684177920q ‘®- 12676750592966656 q 12899470417068032 q '^+ 535928355657089024 q “ 


+450359962737049600 q'*). 



The classical and quantal form factors, are compared in Fig. 
10. The insets show the p integration range as a function of 
the momentum transfer q. Three special cases arc presented: 
(a) excitation, (b) quasielastic transition, and (c) de- 
exciution. For q sufficiently small that Po>P "*■ • transi¬ 
tion is classically forbidden [cases (a) and (c)]. As q in¬ 
creases, the integration limits of Eq. (20) are (po.p*)- The 
increase of P in cases (a) and (c) is due to the effect of 
increasing range of integration overwhelming the back¬ 
ground q " * decrease. With a further increase in q, the inte¬ 


gration limits change to (p _ .p .) which are independent of 
q. Here P decreases purely as q"'. as in (a) and (b). For 
larger q, the more rapid decrease in P results firom a q 
variation combined with the effect of a decreasing range of 
integration, as exhibited for all cases. The transition is agan 
classically forbidden, for all cases, in the limit of Iwge q 
when po>P+- These three overlap situations deKribed 
above are illustrated in Fig. 11. For deexcitation (£ 

Po has always one minimum value pQ=\2m(E~E )] 
When pj>p_. then the pattern of case (c). with momentum 
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limits (pQ,p.), is established. This occurs for de-excitation 
to levels = De- 

excitation to levels n'>n^ is characterized by the panem of 
case (c). When transitions are classically forbid¬ 

den. This occurs^ for de-excitation to final states, n' 
<n.(n.L) = (L/v 2 )[H“(l-L-// 2 ^)^^]"‘^, whose orbits 
are fully within the orbit of the initial (n,L) state. The /t_ 
limit therefore delineates the classically allowed from classi¬ 
cally forbidden de-excitation tr? isitions. The n* demarca¬ 
tions are illustrated in Fig. 11. For excitation, po can be zero 
at q*-[2m{E'-E)y^. so ^hat there is always a range of 
transition momenta q for which Po<p^ i.e. n'>n^(n,L). 
Excitation is therefore always characterized by the pattern of 
case (a). The quantal-classical agreement for transi¬ 

tions is overall very good. 

D. Form factor for n ' transitioas 

The classical probability of transition between states 
specified only by their principal quantum numbers as func¬ 
tion of the dimensionless parameter Q — qa^IZh (as derived 
in Appendix B) is 



which is the classical correspondence of 

® /' *0 m.m' 

( 21 ) 

This quantal form factor (as derived in Appendix C) is again 
a rational function in the momentum transfer q since is a 
summation of Pnin*i* form factors. The classical result is in 
agreement with the result deduced by Vriens [9] from com¬ 
parison of binary-encounter and Bethe treatments of 
electron-atom collisions, and by Borodin [12] from the mi- 
crocanonical phase space distribution 

Exact quantal and classical form factors for 6—► 40 tran¬ 
sitions are compared in Fig. 12. The expression for the quan¬ 
tal form factor for this specific transition 


J 


P (o) = 382205952 10^ (289+14400 177345298871 8867471792636791 

(529+14400<7^)*‘ 

-P 1411176845994965835001817764524075792X 184318171941496097624317093441846656 

X 10* 9^ -I- 93670997716818370857325330958800896X lO’ 9 *-2400432080403014981637748114489344 
X 10‘°<?*-i-34865881226093843259112916256817152X10‘®9‘°-2927694744198493018588901567102976 
X 10'2<7'^+ 14387567274612996874680196399104X 10‘* 9 '^-395454399288571654223098281984 
X 10'’9'®+5' J0746343048968862171136X 10“9‘*-396176923859529631650545664X 1()^°9“ 

16561843 i6737309252780032x 10“ q^) 


is an application of the general expression for transitions 6 
presented in Table L These results are derived in Ap¬ 
pendix B (classical form factor) and Appendix C (quantal 
form factor). Due to the correspondence principle, there is 
excellent agreement between the quantal and more compact 
classical expressions The agreement is also ex¬ 

pected because the characteristic classical singularities in the 
form factor are “smoothed” after the /, /', m, and m' sum¬ 
mations. 


V. SUMMARY AND CONCLUSIONS 

Based on the phase-space description of an atomic sys¬ 
tem, classical expressions for the inelastic form factor have 
been derived. The formulas obtained are the exact classical 
correspondences of the quantal form factors. The classical 
methods quite succinctly reveals important aspects which re¬ 
main hidden in the quantum treatment. An efficient algo¬ 
rithm for calculation of quantal form factors as analytical 


I 

functions of momentum transfer q, for arbitrary initial and 
final states, has also been developed. 

For transitions, the classical method provides 

both the qualitative behavior of the quantum results and its 
physical interpretation. The classical-quantal agreement is 
particularly noteworthy for the integrated form factors (cf. 
Figs. 8 and 9) for inelastic and quasielastic transitions. This 
is because both classical and quantal form factors depend 
only on the overlap of the initial and final distributions in 
configuration space [cf. Eqs. (17) and (18)], so that the clas¬ 
sical singularities apparent in Fig. 5 are averaged to produce 
the smooth results in Figs. 8 and 9. 

The increasing accuracy obtained upon I* integration is 
due to the absence of the multiple delineation of the phase 
space associated with transitions (sec Fig. 3). 

Again the classical picture not only provides the physical 
explanation for the quantal behavior when the momentum 
transfer q and quantum numbers are varied, but also identi¬ 
fies the patterns associated with each type (excitation, quasi- 
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elastic, and de-excitation) of transition (cf. Fig. 10). The 
agreement between classical and quantal integrated form fac¬ 
tors is again excellent In the limit of summing over all final 
states, the total transition probability is 
same result for both quantal and classical cases, which en¬ 
sures full agreement in this limit 

On integrating over angular momentum quantum number 
I for transitions, the agreement is excellent for all q 

even for small qm mtum numbers. This is due to die fact that 
the phase-space region common to the both initial and final 
states (a sphe*^ in configuration space with of radius r*) is 
densely and continuously popular^. 

The classical form factors represent an attractive approach 
for classical collision theory. The form factor is a collision 
kernel to be convoluted according to the dynamics of the 
external interaction causing the transition. Due to the oscil¬ 
latory nature of the wave functions, quantal calculations for 
processes involving highly excited states arc still compuU- 
tionally expensive (in terms of precision, memory, and/or 
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RG. 10. Quantal and classical transition probabilities for (a) 
(4,3)-^8, (b) (4,3)-^4. and (c) (8,3)-»4 transitions as a function 
of the momentum transfer. Insets: the gray area is the integrauon 
range. 






RG. ll. The three overlap situations in momentum space: in 
region (a) there is a value of q for which • in region (c) 

Po>P- for any <1* evon though Po<P^ for some q; in shaded 
region transitions are classically forbidden, Pq>p^ 

time), while classical models are capable of exact results, 
according to the correspondence principles [4], Although 
classical-quantal" comparisons have been made Xo the one¬ 
dimensional hannonic oscillator and to hydrogcnic systems, 
classical form factors can be useful for other atomic and 
molecular systems. The present method would also be valu¬ 
able in determining the response of the three-dimensional 
Rydberg atom to a train of unidirectional short pulses of 
electromagnetic radiation [1]. The classical form 'factor 
methods would be also useful for excited-atom collisions 
[13]. 
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APPENDIX A: MICROCANONICAL DISTRIBUTIONS 

The basic classical probability density for a particle mov¬ 
ing in a symmetrical potential V(r) is given by the microca- 
nonical distribution 



RG. 12. Quantal and classical transition probabilities for the 6 
^40 transition as a function of the momentum transfer. 
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TABLE 1. Quanta! form factor [Eq. (21) from text] for n'—n transitions. 


n ' 

2^ nrqr[{nJn 

1 - 1 n‘q‘ 

2 

3<7‘'n‘-5q^n"-rq=n’ + 9n=/4-23nVl6+15n*/64-9n*<j=/2+65n‘qVl6-llnV/'*- 1 

3 

3 11 14q‘n‘ + 6 < 7 "n*-( 7 -n= + 7l/i=/27-62». /27+I790n‘*/2187 - 243l n*/l9683 

+ 377 n'®/59049-4 n‘< 7^/9 + 290n‘<jV81- l460n*q^/729+ l''43n‘®(j^/6561- l22n®<7*/9+310n*q*/27 
- 442 n '®q*/243 - 412 n*<j‘/27+518 n '®<?‘/81 - 43 n '®q*/9- 1 

4 

3q'''n'‘+I7q‘^n'^+39q‘®n'®+45 9 */i‘ + 25<j‘n® + 3qV-3<7-n- + 45/i^/16-213U"/768+5093n®/4096 

- 167011 /i‘/589824+ 35231 n '®/1048576- 32899 n '2/16777216+ 106183 n “/2415919104+ 41 n*(i^/S -269 n®(j2/256 

- 3475 n''q2/3072+"107449 n‘®(j2/i96608- 128581 n‘2(72/1572864+ 201047 n'*(72/50331648 - 217 /i‘(7*/16 

+ 2691 nV/ 128 -59239 n V/6144+333157 n'2,7^196608- 290087/.'V/3145728-197 16273 n'®<7‘/384 

- 34543 n‘29®/3072+188555 fi'*(?®/196608- 917/I'®9*/16+8027/i'2q‘/256-16035 n'*9*/4096- 239 n‘2q'®/8 

+ 2119 n'^(7'°/256-95 n'"9'2/16- 1 

5 

3<7'*n'*+23,7'*n“+769'^n'*+140<7‘2n‘2+154 9'®/j'®+98q'n'+289®n®-4q‘n‘-5q2„J+73n2/25 

- 28828 n^/9375 + 24084n®/15625-1475182 ii*/3515625 + 29028838«'®/439453I25-4495228n‘2/732421875 

+ 90535316n‘*/274658203125-129O6833n“/1373291015625+1251587n'Vll444091796875+ 56/i^q2/5 

- 13876 n‘92/i875+75224 nV/46875+319042 n'®,7’/3515625 - 7425752 n'2,72/87890625 + 9110956 n'V/732421875 
-8239208/1'‘92/10986328125 + 4510567 n'*<72/274658203l25-4n‘qV5+6836nV^75 - 673004n'®9V46875 

+ 3108484/1 '29‘/703125- 18475444/»“9*/29296875 + 2024596n“9*/48828I25- 546l4636/i"<7*/5493164O625 

- l9l2nV/25+ 178396/.'®9‘/1875-672976/I'29‘/15625+ 10191692/i'V/n71875-22581416/i'V/29296875 
+ 18069764/i'‘9‘/73242l875-914/.'®9*/5+59212/.'29'/375-2347228/i“9'/46875+1541158/.'V/234375 
-2802002 /i'«9''/9765625-5272/i'29'«/25+ 1233332/i'*9‘®/9375- 1272296/i“9‘®/46875 + 2099902/.'*9'®/1171875 
-676/i'‘9'2/5 + 6996/.“9'2/125-17868/.'»9‘2/3125-232/i'‘9'‘/5+ 6076/."9'"/625-167/1'•9'‘/25-1 

6 

3 922„J2+29 920„20+l259'‘/i“+3l59“/i“+5109'^n'‘+5469'2/i'2+3789‘®/»'®+1509‘/i‘-709‘/i*+l59‘n‘ 

- 15 9*/,‘-7 92„2 + 323/i 2/108- 1417/.V432+1230149/i‘/699840- 6667507/i*/l2597120+ 43800179/.‘®/453496320 

- 33636707 n '2/3023308800+ 805306763 n “/979552051200- 8183128999 n '‘/211583243059200 

+ 8488861159 n '•/7616996750131200-4894595653 /|2®/274211883004723200+ 398746151 /j22/3290542596056678400 

+ 35/i*92/2- 18907/i‘ 92/1296+157637/iV/29160-781087/!292/g39go8+1290871 /j'292/25194240 

- 167865689 n '‘ 92 / 27209779200 - 57381211 n '‘ 92 / 48977604560 + 5399280829 /i 'yn052n4m6400 

- 1 17424517/i 2 ® 92/50779978334208 + 2465858993/i2292r 1-^)3964001574400+ 905n‘9*/36+7 n*q*n0i 
-43949/.'®9‘/3888+4792471 n'yii39m- 188652173/i“9Vl51165440 + 479247l /i'29*/839808 

- 188652173 /i “ 9 */151165440+ 77931593 n “9*/544195584 - 32278661 n '‘9‘/3627970560 

+ 264468259 /.2®9‘/940369969I52- 594305017 /|22,‘/169266594447360+ 43385 n '®9‘/324- 1458407 n 'Y/m96 

- 60952999/I'‘ 9 ‘/2519424 - 82299863/I “9‘/22674816+ 470998955/I'’9‘/1632586752 

- 111918155 /i2®9‘/9795520512+ 166167173 n^q*/940269969l52-6225 n'°q*ni+25il55 n'^q*/64S 
-4222535 /:'*9'/23328+ 11088691 /i“9*/279936- 43657751 /.'‘9‘/10077696+ 1017419 n^q*mn916 

-58621799/i 229»/13060694016-6167/.'29'»/9+128209/i'‘9'®/216-2015063n'‘9'°/9720+141006431 n"'9'®/4199040 

- 61592369 n^q'^nSmiAO* 13034627 /i22,‘0/201553920- 14035 n '‘ 9 * 2 / 18 + 170017 n '‘9*2/324 
-23723233/:"9*2/174960+ 12337645/i2®9*2/839808-27182041 n229*2/50388480-5030n*‘9‘‘/9+ 9070l /i"'9'‘/324 
-93505/t2»9'‘/1944+ 737941 /|229*‘/279936- 8995/i"9'‘/36+36175/t2®9'‘/432 

- 37405 /i 229 *‘/ 5 184 - 3455 n^q'*/54+ 14045 n^q'*n296- 2S9n^q^n6- 1 


p{E.L.L:;T,p)dEdLdL^drdp 

. drdp 

= {S(H-E)dE S{\LhL)dL S{L^Z-L,)dQ 

(2irhy 

(Al) 

where the Hamiltonian H, angular momentum |L|, and the 
projection of the angular momentum on z axis L-z are con¬ 
served quantities, and specify the state of the system. Vari¬ 
ous other less restrictive [14] distributions arc directly de¬ 


duced from Eq. (Al) by dropping from Eq. (Al) those S 
functions which correspond with the restrictions on the state 
of the system to be relaxed. For example, when the projec¬ 
tion Lj of the angular momentum is arbitrary, the distribution 
is 

p(E,L;r,p) dE dLdrdp 

drdp 

= {<5(//-£)JEc5(|L|-L)<fL}—(A2) 

(2 7rn) 
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and describes a population of 2L states in all of phase space. 
The physical interpretation is that Eq. (A2) is the number of 
states (phase-space cells) compatible with energy and angu¬ 
lar momentum conservation. The { } factor is a fractional 
number of states in the interval dEdLdL. of about 
(E.L.L-). 

If the system is described in terms of discrete quantum 
numbers, e.g.. the motion is bounded within a finite spatial 
region, the classical distribution is defined in the action - 
angle representation by 

PnimdJdw^S(J^/h-n) SU 2 /.^-(/ + 1 / 2 )) 



^/2m+V(r)-E„) 


pdp 


/II t 


m 




= l/(rr) 


for pi = 2m(E„,-V(r))>pl^„. which ensures real radial 
speeds r, th^n 

j S(h‘-E„,)S(\L\-L)drdp^SiT^L <j> j = 


xSUilh-m) 


dJdw 


Upon action-angle variables (J.w) integration, this distribu¬ 
tion corresponds to a single particle in all of phase space. 
Also l„i„p„i„ = (2vA)~^ is the number of particles in all 
states occupying the unit phase-space element. The phase- 
space distribution for state nl is 


p„,dJdw=S(Ji/h-n)S{J2/h-(l+l/2)] 


X 


dJ\ dJj dJ-j dw^ dw2 dwj 
(2rrft)^ 


Since Jj - y 2 Cos £ z, the integration gives 2 J 2 , so that the 
above p„i distribution corresponds to a population of (2/ 
+ 1) states in all of phase space. The corresponding distri¬ 
bution in (r.p) phase space is then, by changing variables 


p„,drdp=h 


'mJi,J2) 

'0HJ2) 

1 

1 




X 


<5(|L(r,p)|-L> 


drdp 

(2rrft)'‘ 


Since dH/dJi = v„,= r ~,', the frequency (or inverse period) 
for radial bounded motion, then 


p,,(r,p) drdp={hv^,S{H-E^,)hS(\U-L)}j:^^^^- 

(A3) 

This result can be obtained, formally, from Eq. (A2), by 
replacing dE and dL by hv^, and ft, respectively. The sepa¬ 
ration between highly excited neighboring energy levels n 
and rt -1 is hp„i, the Bohr correspondence, and ft is the 
separation between neighboring angular momentum levels 
n,/ and n,/£ 1. 

By noting that 

f S{rp sia 0—L) d(cos $) =— {r^p^ — L^) 

J-i '‘P 


pfUr^po 


Distribution (A3) is thus confirmed as being normalized to 
(21 + 1) states. 

Since dTdp/{2irh)^ ij the total number of bound and 
continuum states, with all quantum numbers, in the phase- 
space volume element drdp^dJdw, then the { } factor in 
(A3) represents the fnu^ional number of states with specific 
quantum numbers. 

For the particular case of Coulomb attraction the energy 
levels are degenerate. The phase-space distribution for an 
hydrogenic atom in the energy level E,, corresponding with 
the principal quantum number n, is 


p,(r,p) dTdp=hv^S{H-EJ^^^^ (A4) 

for bound states of degeneracy n^. The same expression 
holds for states in the continuum if hv^ is replaced by dE. 

The classical distribution p«/(r)=/p^i(r,p)^fp in con¬ 
figuration space is 


p„,(r) r^drdr= 


gi 2 dr dr 
4ir 


rni 


where g/=2/+1 is the statistical weight of the ni level. For 
Coulombic attraction, V(r)* -Ze^lr, r„ = 27r/i^ a.u., and 


4irp«/(r) r^dr-Rniir) dr 


1/22 (/+l/ 2 )^\"''^ 


nit' 


^ \ r. 


dr 


for one nl state (g/* 1). 

Integration of Eq. (A2) over the configuration space 
yields the momentum space distribution p(E,L;p) 
= /p(£,L;r,p) dr. Then 


p(E,L-.p) p^dp dp 


gi 

2irh'‘ 


2 pdp 


r (p^-LVr?)‘'*|V"(r,)| 


4ir’ 


(A5) 


vhere r,- are the roots of p^=2m(E-V(r)] for a given p 
md V'=dV/dr. The radial momentum distribution, p^i 
=p(E,L)hv„,h, for bound hydrogenic states, reduces with 
j,-lto 


and 
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Pnlip)^T^P"^P^ 


2Ze"dp\ 

I 





X 


2miE!L“ I)“ \ 
(2mZe‘)" jc- I 


where x = p/(2mlEl)^^, Since £= -Z‘(e^/ao)/2/t^, then x 
= plp^ where p^^Zp^fn is the characteristic momentum in 
the Rydberg orbit n, po-hlao is the atomic unit for linear 
momentum, and gq the Bohr radius is the atomic unit for 
distance. With L = (/+l/2)A, the hydrogenic momentum 
distribution is 



FIG. 13. The basic geometry for calculation of . 


4 \ 

p„,ip)Anpdp = -n^^^n2p2)2y- 


/+I/2 

n 


X 


2nP / 


-in 


where P = plZpo. This compares formally with the quanta! 
results and is useful for the nl-^n' classical form factor. 
Another formulation of classical momentum distributions 
was recently presented in Ref. [15], 


APPENDIX B: CALCULATIONS OF CLASSICAL 
FORM FACTORS 

The classical form factor for transitions between energy 
and angular momentum bands + is 

P^^^)^P{EM\E\U\^)dEdE* dLdV 

--(2TThy^^M)dEdE^ dldU 

in terms of the phase-space integral 


When £jf is summed (integrated) over all L states, then 
SR,^dLdU^\ and 

J R,fdV = )^^ ^<5[rp|Sin^-L]d(cos^) = L/(mrVi'^). 

(B3) 


1. Form factor for EL^E*V transitions 

Evaluation of k\f for fixed L and U is facilitated by not¬ 
ing, from Fig. 13, that 

cos 5= cos $12 cos $+ sin $12 sin $ cos (f> 

in terms of ?($,</>) and the fixed angle $12 between p \ and 
Pa- 

On changing the <f> variable to 5 in Eq. (Bl), the 4^ inte¬ 
gral, for $ fixed, is 

fair , , 

S{rp2sind-U)d4>=^2 


AiKq) = J dTdpSmr,p)-E)S(\Ur,p)\-L) 

X ( 5 (//(r,p+q)-£')< 5 (|L(r,p+q)|-L') 

of S functions involving states i-(E,D and/=(£',£'). 
For transitions between bound states the transition probabil¬ 
ity is obtained simply by the replacements dE-*hv and dL 
The integral can be recast in terms of the radial integral 

1 fi 

4:;;J^^<5(rp,sin e-L) 

tilt 

xd(cos $) j S{rp 2 sin d-L')d<f>, 

(Bl) 


X <5(rp2 sin 5-L')d(cos 3). 

The factor 2 originates from the fact that, as 
-►27r, the range (1^-^ 12 !.^+^ 12 ) in 3 is covered twice. 
The function 

5 = sin ^12 sin ^sin ^ 

is expressed as a function of 3 and 6 variables by 
5^( 3, e\ $ 12 ) - sin^ 6 sin^3- (cos ^, 2 “cos 6 cos 3)^ 
Subsequent integrations arc facilitated by noting that 

S(Fix))=^——, f(x,) = 0. 

- |F'(x,)| 


where d and 5 are the angles between r and ^ and P 2 - 
respectively, as 

AjKq)= J <^Pi <^P2 ^(q+P|-P2)<5(W(r.p,)-£) 

(B2) 


Hence 

(sin d)S(rp 2 sm 3-L') = 




mp2r^r* 


2 <5(3-3,). 


i- I 


where the roots. 3 ,<:t/ 2 and ^ 2 = w- 5,. are given by 
sin 5, = L7r/72. cos 5i=-cos 52- 
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r2n 

S[rp2sin6^ L']d<t>=^ 
Jo 


mp2r"r' 


x[s:'(di,e)+SZ'{0\.9)l 

where 

5i(5,, ) = sin- d sin- 5i - (cos 5,2=cos 5 cos 5, 

Upon 6 integration. 

2L' 


The integral ^,KPi .r) is then 
ZLL' 


4iri?i,(li,;r) = 


m^r*rr' 


[C:f(r.q)^G-f(r.q)]. (B4) 


where 


4ir/?iKPi.P2:'’) = 


2L' n 

mpjr'r'J J-i 


Gif(r,q) = 


[s:'(5,.5) 


y/(q^-AiHBi-q‘)' 


Since 


+ 5Z'(52,5)](5(rp| sin d-L) d(cos 6) 
2LL’ 


mr^l 2 =E- V(r)-L~l 2 mr‘ 


m‘PxP2r*rr' 


[S:‘(5i.0i) 


-t-S:‘{5,.52) + ^*‘(5,.5,) 

+ 5:'(3,.52)]. 


holds for the initial and final states, the 5, A, and B functions, 
and hence integral (Bl), are ail independent of pi. The tran¬ 
sition integral (B2) is then 

AiKq)* ^TrR J,r,q)Il{r,q) r^dr, (85) 


where sin =sin /rpi.and cos $ 2 --cos Bi. From where the only p integral is 
these relations and from the above definitions of the func- 

n(r.?)= I S{HiT,p)-E)SiH(r.p+q)-E')dp. 

S-(5,.5,) = S*(5,.52). 

S-(5,.52)=S*(5,.5,). 
with the result that the radial integral is 
4LL' 


Hence 


4irRiKp,.P2:r) = 


m^P\Pir*rr' 


[s;'(3,.5,;5,2) 


+S_'(9i .5, ;5,2)]. 

Upon P 2 integration in Eq. (B2) 

^if<Pi:'•)=! J^u<Pi.P2:'’)<5(P2-'(Pi+qi)MPj* 
5. are evaluated using the substitutions 

sin5=Z./rp,, cos5=mf/p,, 
sind=L'/rp2, cosd*mr'/p2, 
and 

cos5i2=(pf-t-P2-<7^)/2piP2 
to give, simply, 

2p,P2S>(dx.ex)’‘^(q^-Ai)(Bi-qb. 
expressed in terms of the momentum-change limits 


n(f,^) = 2'7r[ p^dp S{p^/2m + V(r)-E] 

Jo 

f* - ‘ 

'< j S{pq cos 9/m -E'-E ' 7 -q^l 2 m] a{cos 9), 

where 9 is the angle between p and q. There is only one 9 
root provided that 

p»Po=‘m/q\E—E' +q^/2m\. 

The p integrations therefore yield 

2irm^ 

n(r,q)« —^— B(r,q), 

where 6 is the step function having a value I if V(r)^E 
- pjpm is satisfied and zero otherwise. The final expres¬ 
sion for the classical transition probability density (13) for 
E,L-*E',L' impulsive transitions is therefore 


P(£.L;£'.L':q) = 


1 2irm^ 


(2irt)^ R 

X J 4 wkj,r,q)B(r,q) r^dr. 


(B6) 


Ai = mHri:r')^ + {L-L')^/r^ 


where is given by Eq. (B4). Hence 
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167rLL' 

(27Th)^q 



p{r) = 2m[£- V{r)]^pQ--^\E* -E-q"^l2m\. 


^G-f{r.q)]@{r.q). 

For hydrogenic systems, the condition V(r)^E-pl/2m is 
satisfied for any q and all r within the radial region H (cf. 
Fig. 3). Thus the step function 0 is always unity. 

The probability oi n A—^nW transitions due to an impul¬ 
sive momentum change is then 

X P(£„.{/+ l/2)ft:£,..(/' + l/2)ft;q). 

2. Form factor for EL-*E* transitions 

On using Eq. (B3) in Eq. (B6) the probability density for 
E,L~^E' transitions is 


This probability density is exact for all V(r) and has a 
simple physical interpretation. The transition probability is 
given by VAp/(2ir^)^, which is the number of states in the 
“reaction** volume V=4/3'nrr^ multiplied by the volume of 
momentum space Ap consistent with energy conservation. 
The initial- final-state energy conservation equations a” E 
= p^/2m + V(r) and £'=£+q-p/m + <?V2m, respectively. 
Then 

f2v 2nm" 

Ap= dp—2vp^dpd{cosS)— - dEdE\ 

where the z axis is along q. 

The classical form factor for transitions is 


P(£,L;£';q) = 


— (,vr)~'e(r.q)dr. 
irh^qJo 




For Coulomb attraction V(r)= -ZcVr, then 


where v(r) is the speed along the initial trajectory. Since 
p^=2m[E- VCr)]. integration over r may be replaced by p 
integration and drlv=dplV’(r) so that 

/>(£.£;£';<?)=-^2 f \rjV'{rj)\-'dp, 

Trn^q j Jpo 

where rj is the root of p^ = 2m[£“ V(r)]. In terms of the 
momentum distribution (A5), then 

mL r* 

£(£,L;£^q)==->- p(£.L;p) 


r^{q)=i{Ze^m)q\q*-^m(E+E')q^ + dmHE-E')^]-' 

so that the transition probability, with the substitution q 
= Q(Zhlao) now becomes 


Pn.n’iQ) 


3ir{nn') 


’\y 


Q* + 2 






which for bound nl states is in agreement with previous re¬ 
sults [8.16,6]. 

The classical form factor for nl—*n' transitions in hydro¬ 
genic systems is, in atomic units. 


Pm.A<1) 


2 ( 21 + 1 ) 
irqn'^ Jpi^ 





£(l+nV) \^ 

2n^p I 


- 1/2 


dp. 


where the integration limits are mven by two conditions: the 
integrand is real and 1/"'^“ l/n^i/2q. 


3. Form factor for E~*E' tmaMoat 

Using Eq. (B6) with Air= 1. the probability density for 
£-*£' transitions is 


P(E,E'-.q) = 


2irm^ 




where is the largest r which satisfy the condition 


APPENDIX C: CALCULATIONS OF QUANTAL 
FORM FACTORS 

The quantal transition probability (5) for the hydrogen 
atom is in general a rational function of the momentum 
transfer q because Eq. (5) with ^it) = Rni(r)Yi„{r) can be 
decomposed as 


/+/' 

/^,/«-('(q)“(2/+l)(21' + l) 2 Ai(q)H2L+l) 

L 1 /']* 

^[0 0 0 

L 

where {• • •} is the Wigner 2j symbol, and fnifq) is the 
radial integral: 

Ai(q)=f Rn,{r)R„,,'{r)jd(ir)r^dr, (Cl) 

n't’ JO 

where is the modified Bessel function. Because the radial 
wave function has the simple structure e~''''dQ(r), 
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where 2 is a polynomial of order n~l - I, integral (Cl) where the polynomials 
contains only terms with the form 


[ e'"' djd<!r)dr, 

JO 

where a is !/«+ l/n' and /t is an integer number (k-l^V 

+ 2.n + rt') greater than L This integral is the following 

rational function in q: 

/t.t,(a-<?)-(2L+ IK2L- 1)!! 


f„{x)= g,(x)= ,f,(2-/t.3jc/n) 


are given by the degenerate hypcrgeometric function 
where the first argument is a negative integer. Finally, 
is the integral in jc, y variables: 


ri^f [ dxdysm(q(x-2)/2) 

24nn Jo Jo 

X(x^+4xy+y*)/*,(x,y)/*,.(x.y)/(x-y). 




k-L 

~T~ 



k-L 

~2~ 


+ 2-^+2’" 



since the hypergeometric function i is a polynomial when 
cither the first or second argument is an integer. This proves 
that integral (Cl), and hence the quanta! form faaor (14), are 
rational function in q. It also provides the practical procedure 
to calculate the quantal probability (14) in an analytical 
form. 

The quantal transition probability (21) can be written in 
terms of the density matrix element, p,(r,r') 


The observation that i>,(x.x)=0. means that x-y is a divi¬ 
sor for PJx.y). On writing sin in exponential form, the 
integral contains only primitive terms of the form 
x^i/ or g)(x)e"“=', with various positive integer powers 
m and a. = ( l/n+ Un'±iq)/2. The elementary integrals 



X 2 F'i(n+ l,l-m,2,l/an) 


PnAg)- j j t/rdr'eW^ ''p,(r.r')p:,(r.r'). 

The density p„ is the residue of the Coulomb Green’s func¬ 
tion [17] and can be calculated from 


and 


/: 


x"g,(x)e"“dx= 


/ an — 1 ^ 

"m! 

" 1 

1 an j 

' a ' 

[an-lj 


X ,(nT.1.2-m,3,l/an) 


I Pn 

p„(r.r')= lim (£-£,)2£(r.r')=-JTTT 

£-£, ‘•mao •* y 

in the spatial variables (x,y) given by 

x = Z/ao{r+r'+p), y=Z/ao(r+r'-p) 

where p = |r-r'|. The function P, has a simple structure. 
Because 

where M is Whittaker’s function, then P„ is simply 
exp[-(x+y)/2n3x polynomial in x and y. Thus 


Pn(x.y) 




2n^ 


-[2n(x-y)/;,(x)g«(x) 


'^{n-l)xy{fn{x)gn{y)-gnix)fniy))l 


arc rational functions in a which is linear in q. The form 
factor Pnn» is therefore a rational function of q. The proce¬ 
dure described is remarkably efficient since it reduces the 
multiple integrations to a finite number of symbolic opera¬ 
tions by (a) recognizing the primitive terms, and (b) replac¬ 
ing them with the appropriate elementary integrals. Based on 
this procedure, the results of calculations for the transition 
probabilities for /i' = 1.2,3,4,5, and 6 and arbitrary n are pre¬ 
sented in Table I. This illustrates the power of the method. 
The form factors for transitions from L, and Af shells were 
obtained by Bethe and Walske [18]. The form factors for 
transitions to continuum states with wave number k arc ob¬ 
tained by analytical continuation replacing n with i/k. The 
dipole oscillator strengths for n-n' transitions, 


/«„, = 2A£..«.lim 




q^o q^ 


can be readily deduced from the results presented in Table I. 
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The quantal impulse cross section is derived in a novel form appropriate for direct classical corre¬ 
spondence. The classical impulse cross section is then uniquely defined and yields the first general 
classical expression for nf - coilisional transitions. The derived cross sections satisfy the optical 
theorem and detailed balance. Direct connection with the classical binary encounter approximation is 
also firmly established. The unified method introduced is general in its application to various collision 
and recombination processes and enables new directions of enquiry to be pursued quite succinctly. 


PACS numbers: 03.65.Sq, 3l.l5.Gy. 32.80.Cy. 34.60.+ z 

The quantal impulse approximation (QIA) is an estab¬ 
lished and productive cornerstone in general scattering 
theory [1]. It is based on the assumption that an inelastic 
collision of incident particle i with a two-particle sub¬ 
system (2,3) results from each binary (i - j) scattering 
under the two-body interaction Vij alone, leaving the 
“spectator” particle k unaffected. The internal interaction 
Vjk between the “active” priticlc j and the spectator 
particle k is ignored duri the {i — j) collision except 
insofar as it generates a d* zibution p{pjk)dpjk over the 
relative momentum Py* of particles j and k. This approach 
has been invaluable in nuclear physics for high-energy 
neutron scattering by complex nuclei [1] and in atomic 
physics for atom-Rydberg atom collisions at thermal 
energies [2—7] and for electron (ion)-atom and atom- 
atom collisions [8,9] at high energies. It is also valuable, 
for example, in the study of three-body recombination, 
Li + Li + Li —" Lii + Li which limits the density and 
lifetimes of Bose-Einstein condensates at ultralow energy, 
and + p + e*** W + e*** for the formation of 
antihydrogen. 

A classical binary-encounter approximation (BEA) has 
also been formulated for Rydberg collisions [3] and for 
high-energy coilisional ionization [8-10]. Although QIA 
and BEA are conceptually connected, the formal interre¬ 
lationship between them is quite complex and has never 
been firmly established. Previous studies [6,8,11] have de¬ 
pended on fairly complicated theoretical analysis in an ef¬ 
fort to reduce QIA to BEA. 

0031 -9007/00/ 85( 1)/1(5)$ 15,00 


In this Letter, the quantal impulse approximation is pre¬ 
sented in a new form which provides quite naturally “the 
royal road” to classical correspondence. A classical im¬ 
pulse approximation can then be uniquely defined. The 
novel expression yields the first general classical cross sec¬ 
tion for nf - transitions and the standard result for 
— n' transitions. The derived cross sections satisfy 
the optical theorem (which implies probability conserva¬ 
tion) and detailed balance. The method introduced then 
permits BEA formulas to be derived quite succinctly and 
in a unified way. There is substantial renewed [7,10-15] 
interest in the power of classical dynamics in almost all 
fields of modem physics, attributed to the desire [15] to ob¬ 
tain a more thorough understanding of the classical-quantal 
correspondence. This present development provides a pro¬ 
found and important quantal-classicd connection in colli¬ 
sion physics. The new formulation also provides a natural 
origin for development of new semiclassical methods. 

The transition T-matrix for the free-free {pij plj) 
scattering under the two-particle (i - j) potential Vij 
alone is 

where 'PCpj/.r.y) is the exact wave function for i - j 
scattering in the (i - j) center-of-mass frame CM(i,y). 
The momentum transferred from i to j is 

qy = Pi - Pi “ Py - Py “ Piy " Ply • 
where the momentum of particle i is p, in the fixed 
CM{ijk) frame of all three particles, and is p,; in the 
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CMIi.y) frame. The primed (') and unprimed momenta 
denote postcollision and precoliision values. 

The quantal impulse T-matrix for i - (j, k) scattering 
decomposes as the sum T,/ + T,/ of binary 

collision matrices Tlf\ with j — 2 ot 3. Each T// in the 
momentum-space representation (MS) is the overlap 
[1.6.7], 

X iro(p.7.p'y)lfl^/.ii/*)>P,.. (1) 

of the free-free Tij matrix with the initial and final mo¬ 
mentum wave functions 

4,ipjk) = • rji/fi)drjt, 

for the (J,k) subsystem. The momentum transferred to 
j k relative motion is 

/ _ -Jk—n 

q/t = P;* Pj* - ^. + Af t *>> ’ 

where pjk is the relative momentum in the initial CM(y, k) 
frame. After the (i - j) collision, CMC;.*) moves with 
the residual momentum Mjflj/iMj + Alt). The bi¬ 
nary rjf^ matrix has also the configuration space (CS) 
representation 

^^T!/\pi,pl) - <exp(ip- • ti/n)4tf(rjk) 

X |V<7(r,7)|^r(PM*’.-r7t)>r,,^. (2) 

The impulse approximation [6] for the initial total system 
wave function in (2) is 


'Pr " 'Pimr(Pi;r,.r7*) 

- (? j 0 <(P;*)*(Pi 7 .PC * 17 . r*) dpjt , 


where the active and spectator particles arc j and k, re¬ 
spectively. The free-free (i - j) scattering wave function 
is 'Pfp.y.r,;) in the CM{iJ) frame and is 

‘t>(p. 7 .p*;ri 7 .r*) - exp(ipt • rk/n)^{pij,rij) 


in the fixed CM(ijk) frame. With the aid of the 
(r 7 *,r,)-(r, 7 ,ri) transformation equations for the inte¬ 
gration variables in Eq. (2), of the identity pi • rj + 
Pyt • fjk ■ Pt • P* + Py ■ Py» “ 

(2irfi)~^^^ /^(P/*)exp(*P/* * Tjk/h)dpjk, it can be 
shown that the momentum-space and configuration-space 
representations ( 1 ) and (2) of 7</ are equivalent. 

The quantal differential cross section is 




where Mi - AiiiAij + M*)/(Afi Mj + Mk) «the re¬ 
duced mass of the / ~ (7 + k) collision system and Pi “ 


are the momenta of / with v/’ measured with re¬ 
spect to initial CM( 7 , k). Although QIA assumes that each 
scattering is a separate event, the cross section (3) involves 
iTi/P “ which exhibits interference be¬ 

tween the two-body (/ — 7 ) scattering amplitudes. For 
three identical particles ITz/P “ 4|r^^’p. The classical 
representation of Eq. (3), however, involves only |r,/|* = 
which is valid for binary and indepen¬ 
dent scattering (as, e.g., in Rydberg collisions). 

The key idea of this paper is to find a phase-space rep¬ 
resentation for ’P. This is accomplished by writing 

p.OP “ f ( 4 ) 

as the product of the momentum-space ( 1 ) and the con¬ 
figuration-space (2) representations. It can be shown that 
Eq. ( 4 ) with Eqs. (1) and (2) for (i - j) scattering alone 
reduces to 

irS'Vpi.P.OP - (2^«)'/ dp f drp}(r.p + q) 

" X |r,7(p,7,p;y)Pp/(r,p). (5) 

where (Pyt.Pyt.qy*) are now denoted by (r.p.q). The 
internal (j - k) quantal distribution in phase space is 

p(r.p) - (2trfi)-^/V(r)e"‘'’’'^"<^'(p). (6) 

in agreement with that [12] previously defined. The 

standard probability densities p(r) “ /p(r,p)dp == 
ltA(r)P - »» confiMration space and 

p(p) ” /p(r,p)dr — l 0 (p)P “ «^,r(p)|Triii(p)l^ 
momenttim space are recovered directly from ( 6 ), which 
is the standard ordered version [12] of the Wigner distri¬ 
bution [16], In terms of the phase-space distributions for 
the initial and final internal j ~ k states, the differential 
cross section (3), with Eq. (5) for independent scattering, 
separates into (i-2) and (i-3) componenu, each given by 


dp\ 


X {p;(r.p + q)l/.7(Pi7.Po)Pp.(r.P)}. 

(7) 


where the (i - j) and i - (7 + *) reduced masses are 
Aiij and 2 M<, respectively. The scattering amplitude, 
/. 7 (P. 7 .Py) ■ fijig^PUil for free-free (/ - 7 ) colli¬ 
sions is a general function of momentum change qj and 
of the i - 7 relative velocity g v, - vy. It therefore 
depends implicitly via g on the j ~ k relative momentum 
p. Thus, (5) and (7) are the free-free transition probabili- 
ties/diffeiential cross sections, respectively, averaged over 
the initial and final phase space distributions. 

The above novel expression (7) for the quantal impulse 
cross section ( 3 ) has the phase-space representation ap¬ 
propriate for direct classical correspondence, obtained by 



VOLUME 85. NUMBER 1 


PHYSICAL REVIEW LETTERS 


3 JuxY 2000 


replacing the quantal densities pij of (6) by the corre¬ 
sponding classical phase-space distributions [12]; 

p'(£.L;r,p)drdp = {S{p^/2m + V(r) - E)dE} 

( 8 ) 

for a (j k) symmetric interaction V{r). The 
internal Hamiltonian H{r,p) = p^llm + V(r) of the 
j — k subsystem (with relative momentum py* * p 
and reduced mass Mjk * m) and the internal angulr- 
momentum L(r,p) = r X p are conserved quantities in 
phase space. The distribution (8) holds also for bound 
states (n, () but with dE and dL replaced [12] by hvni and 
h, respectively, where the frequency (or inverse period) 
for bounded radial motion is Vni distribution 

(8) is normalized to j p^^ dr d^ * {li + 1) states over 
ail of phase space. 

(a) nf —► transitions. —The angular f(^r.0r) in¬ 

tegration, resulting from Eq. (8) in Eq. (7), involves the 
angular momentum overlap integral, 

Rif(r,q) = ' - [ Sirpsindr - Ijdfcosd,) 

^ 47r J^\ 

rlir 

X / S[rp^sindr “ L']d<i>r. (9) 

Jo 

where cos^r = p • f and cos^r * p' * f with p along 
the 2 axis and with p^ ** p + <1 fixed in the XZ plane. 
Integration via the techniques outlined in [12] yields 

^3 

Rif{r,q)dEdLdE^dV * ^ 

X [Gtfir,q) + C-f(r,q)]. 

( 10 ) 

_ ’..-J 


The classical radial density and the classical radial proba¬ 
bility for continuum states of energy E are 

P£r('’) “ j P«^(r,p)dp = 47rLdEdL/[{27rfi)^r-r]. 

and ^npEfr-dr = (2LdL/n^){2dEdt/h\ respectively. 
The bound state distribution is p^t ” [(2f ^ 
l)/47r](2f/Trte). Each of the contributions, 

G,7(r,^) = [(c?2-Al)(Bl (H) 

must be real, so that q must then lie within the classically 
accessible range given by 

Al(r) * m\r ± r')^ + (L - ^ 

— m^(f ± r')^ + (L + 

For a given momentum transfer the r integration pro¬ 
ceeds over the various [12] radial ranges 'Riq) within 
which the G,y arc real. The initial and final radial speeds 
f and r' determined by 

„ ^(/) ^ . LW/2mr^ 

are functions only of r for given initial nl and final n'/' 
states. Hence, Gff{r,q) and the angular integral Rif{r,q) 
depend only on r and q and not on p. The remaining p 
integration in Eq. (7) therefore involves only the overlap, 

j’5[//(r,p)-£]|/.>[g(p).9;]P 

X 6[/f(r,p + q) - £']<fp, 

of the i - j free-free differential cross section with the 
initial and feal energy distributions at a fixed r. When 
integrated [17] and combined with the angular momentum 
overlap integr^ (9), the cross section for - n't' tran¬ 
sitions resulting from (i - j) collisiqns is 




_27r 

Wi 


>_/•«. f f 1 rr-lM 

PHW('ir.q)dr — j 

fjvf J,,. Jv.(q) 1 ^ JgUr.q) 


\fiji8,qj)W 


[ig^ - gDigl - 


( 12 ) 


The classical probability that the i - j impulsive transition transfers momentum q to the j - 1: system with internal 
separation r in the interval r, r + dr is 


PntMr,q)dr - i^^\A‘irRifr'^dr dE dLdE' dV] - {*7r)^pUir)pn't'(ry ^[C^fir.q) + G-f(r,q)]Q(r,q). 

^ ^ (13) 


The classical limits to the g^ integral in Eq. (12) for fixed [ 
{r,q) are 

gl{r,q) - vf + vj(r) - 2u,V;(r) 

X cos[tf’(r.^) i: d,(q)], 

where, for given transfers 1 ^ E' — E and q in energy 
and momentum, the angles that p/ and p each make with 
q are, respectively, determined by 


CQsdi(q,T) — (2MiX + q])/2piqj, 

coiB'{q,'E) - (2m I - q^)l2pq, 

in accord with energy and momentum conservation, 
implicit in the 5 function and impulse conditions. 
Since Icos^j s 1, the {j,k) relative momentum 
p > po ■ I2mr - q^l/2q| so that the step func¬ 
tion © in Pnt.n't' has value 1 for p(r) a po and zero 


3 




VoLO'ME 85, Number l 


PHYSICAL REVIEW LETTERS 


3 Jlxy 2000 


(14) 


otherwise. Since |cosd,| ^ 1. then q- = Ip, - p,l ^ 
qj s Pi + Pi = q+, as expected. 

The integrated probability for the impulsive (fif - n'C) 
transition, 

PniM'Aq) = f Pnuc(r.q)dr, 

is simply the classical form factor as defined in Ref. [12]. 
Equation (12) therefore illustrates the interesting manner 
in which P„i.nu- becomes deconvoluted over r when 
the general (i - j) scattering cross section <r(g, qj) 
depends additionally on the (i - j) relative speed g. For 
g-independent cr, (12) reduces to the stand^d impulse 
result [1,6] involving the full form factor (14). 

Since the integration limits and /*,./ are symmetrical in 
i and /, the basic cross.section (12) satisfies the detailed 
balance relation pjernweipt) “ “Semi¬ 


Pn(.E'(p>Q)dp “ ^j^^^^[dvR'ifr^drdEdLdE']^ [4»rp^<(p)p^dp], 


quanta!” probabilities may now be defined by replacing 
either or both of the classical radial densities p'^ in (13) 
by their f-independent quantal equivalents, 

pUr) = + l)R„^(r)/4rr, 

m 

for both bound ni and continuum ei states. 

(b) n( — n' transitions .—On integrating Eq. (9) over 
all final angular momenta L\ then 

fir* , \ dEdL 

R'ifdEdL “ jy J 5[rpsind, - L]d(cos^r) '"“-j 

-T-^P,Mr). (15) 

Airmp 

The classical probability for the nf - {£', E' + dE') im¬ 
pulsive transitions in the {j.k) subsystem with internal 
momenta in the interval (p.p + dp) is then 

mdE ', 


since 4np„({p)p^dp ■ Avp„t(.r)r^dr under the energy constraint pHr) - 2m[£ - V(r)] for the initial classical dis¬ 
tribution. Note that m/q - Mj/qj- The classical impulse cross section for n( - (£',£' + dE') transitions is therefore 

2ir f' /■* „ , w f 1 r**!''-’! 

o-i<,£-(p.) “ 7^ / Ijdqj J Pnt.E'(p.q)dpj- J 

MijViJq. J Poiq) {^Jti.iP-<t) 


[{g^ - gDigl - 


The cross section for n( - n' transitions is obtained by substimting hvn> for dE' in Pnt.E‘- By adopting the quantal 
distribution (2f + l)d>h(p) for 4irpJ,(p), then Eq. (16) provides semiquantal cross sections [6]. 

(c) Relationship with Optical Theorem and Classical Binary Encounter Theory .—Integration of Eq. (16) over all pos¬ 
sible final £' states is facilitated by recognizing that the required integral. 


rr- 2g/qjdE' 

Jtl ~ g-)(g+ - g^)Y'^ 

is simply rr. The classical limits T^iq; vj • Vy; v,, vj) to 
the energy transfer I, for the prescribed 0xed arguments, 
need not then be specified. The cross section for ail elastic 
and inelastic transitions is then 

<rni\vi) m ^cr„(j,.{vi) + f (r„e.t(vi)de 

i.e., the rate Vitr^t for transitions to all states is therefore 
the rate for all i - j binary encounters, averaged over the 
initial momentum function fory in the field of k. Expres¬ 
sion (17) simply restates the optical theorem, which, when 
applied to (1), provides (17) directly. 

Expression (17) is the original basis [8] of the stan¬ 
dard classical binary-encounter approximation (BEA) for 
energy-changing collisions, where the p'ij- region of in¬ 
tegration is constrained by the required energy change. 
Previous BEA semiquantal results [3,6] for n( — n' col- 
lisional transitions can then be shown [17] to be identical 
with (16). Generalization [17] of BEA to cover the an- 


Y*-__ 

’ Je. [d - ’ 

^ gular roc.j.ientum changes in n( — n't' collisional transi¬ 
tions is rreJe simply by replacing !pntip) dp in (17) by 
the integral /p„f(p,r) dp dr over phase space. See also 
Refs. [13,14]. 

In summary, the quanta! impulse cross section (3) 
has been presented in a valuable new form (7) which is 
the appropriate representation for direct classical corre¬ 
spondence. The classical impulse cross section has been 
defined by (7) with (8) which yields, in quite a succinct 
fashion, the first general expression (12) for the classical 
impulse cross section for nt - n't' and nt — et' elec¬ 
tronic transitions. The cross section satisfies the optical 
theorem and detailed balance. Direct connection with the 
classical binary encounter approximation (17) has been 
established and the derived nt — n' and nt - e cross 
sections (16) reproduce the standard BEA cross sections 
[3,7]. The present unified method can also furnish 
semiclassical impulse cross sections, obtained simply 
by adopting semiclassical phase-space distributions (at 
present, unknown) within Eq. (7). Although applied 
here to a time-independent formulation of collisional 
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electronic transitions in an atom, the prescribed method 
is general in that it can be applied, where appropriate, 
to atom-molecule rovibrational collisions, to three-body 
recombination, and to explicit time-dependent problems 
as laser-pulse excitation [18]. The method presented also 
helps elucidate, quite succinctly, the role played by the 
quantal-classical correspondence in collision dynamics. 

This work is supported by grants from AFOSR: F 
49620-99-1-0277 and NSF: 98-02622. 
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Abstract A new exact solution of the time-dependent quanial equation is obtained for the full 
array of angular momentum mixing transitions nt -♦ nt in atomic hydrogen induced by collisions 
with charged particles at ultralow energies. Based on this new solution, efficient numerical 
procedures are devised. It is proven that the present (fixed-firame) solution is equivalent to the 
rotating-frame approach described by Kazansky and Ostrovsky (Kazansky A K and Ostrovsky V N 
1996Pfiyr. Rev. Lett. 77 3094) and that it overcomes the difficulties therein. Analytic expressions 
for low quantum numbers n are presented. Numerical resulu for the transition array with n « 28 
are reported. 


Stark mixing occurs when the electron of a Rydberg atom (in a state with principal quantum 
number n) changes its angular momentum without changing its energy, as a result of a 
collision, at large impact parameter b, with a slow massive particle of charge Z\€ moving with 
velocity v. It is important in many areas of atomic physics, as in the Auger (or autoionization) 
process which follows the collision between ions and atoms (Miraglia and Macek 1990), in 
ZEKE spectroscopy (Merkt and Zare 1994), in astrophysics (e.g. Percival 1983), in recent 
efforts (Menshikov and Fedichev 1995) to produce anti-hydrogen at 4 K and for general 
three-body recombination (Flannery and Vrinceanu 1998) at ultralow energies. The first stage 
in ultralow energy electron-ion recombination (Flannery and Vrinceanu 1998) at temperature 
Tc is a very rapid collisional capture into high Rydberg states with high angular momentum t 
and large radiative lifetimes at a rate proportional to Thus the mixing is an essential 
step in producing the low-angular-momentum states required to radiaiively decay at a relatively 
high rate to low levels, thereby stabilizing the recombination. 

On considering the Rydberg atom in a frame rotating with the intcmuclear axis, the Stark 
mixing problem can be reduced to the problem of the Rydberg atom in mixed static fields: 
electric, provided by the projectile ion, and magnetic, produced by the non-inertial (Coriolis) 
forces. In this way, the well known equations, in both classical (Bom 1960) and quantum 
(Demkov et al 1970) mechanics, for the problem of interaction between weak fields and an 
atom, can be adopted to provide, in principle, a solution for the Stark mixing problem. Both 
quantal (Kazansky and Ostrovsky 1996a, b) and classical (Kazansky and Ostrovsky 1996b, c) 
versions of this approach have succeeded only for f = 0 to higher angular momentum t' 
transitions, appropriate to the experiments described in a paper by Sun and MacAdam (1993). 

This letter presents a new exact solution for the Stark mixing process, valid for transitions 
nt nt between any states within the shell of energy £«. The present theory in the 
fixed-frame representation is shown to be formally equivalent to the rotating-frame approach 
(Kazansky and Ostrovski 1996a, b), but, in contrast to it, the full array of transition amplitudes 
can be obtained at once by efficient numerical procedures. Results for transitions within the 
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n = 28 energy shell are presented. Analytical formulae for the transition probabilities are 
possible for low n. 

The trajectory of the projectile, initially moving along the positive Z direction, is assumed 
to be confined in the YOZ plane. In addition to the Hamiltonian Hq = p^ITjn^ — e^/r, the 
angular momentum L = r x p of the unperturbed Rydberg electron and the Runge-Lenz (or 
eccentricity) vector. 


A 



r - (p • r)p - m,e 





(1) 


directed toward the pericentre and normalized to angular momentum units, arc also conserved. 
Here = >/-2me£^. These quantities define the dynamic 50(4) symmetry of the hydrogen 
atom. Because the 50(4) group is isomorphic with the direct product 50(3) 0 50(3) of two 
rotation groups, a special decomposition, L — Af +iV and A = Af— N. permits the dynamics 
of the hydrogen atom to be separated into two decoupled motions. The generators Af and iV 
act independently as angular momenta and arc also conserved quantities for the unperturbed 
Rydberg atom. 

The orbital electron interacts with the time-dependent electric field €{t) generated by the 
passing projectile of charge Zi*^. In the weak-field approximation, this field is constant over 
the spatial extension of the atom. In this ap^ximation, which is the same as the dipole 
approximation, the interaction potential V =er •£ is 


V(r, R) 



Z\e^ d^ . 

r* 

vb dr 


^l^^{ysin^ + zcos<l>) (2) 
vb dr 


where R is the intemuclear vector and is the angle between R and the OZ-axis. The 
angular momentum of relative motion Lni = = fivb^ where fi is the reduced ma^ of 

the projectile-target system, remains conserved since Li^ L, so Lni and L are effectively 
decoupled. A classical trajectory for the relative motion is then valid. 

If the projectile moves very slowly, the orbital electron adjusts itself adiabatically to the 
ion perturbation and no energetic transitions occur. In this limit, Pauli’s replacement rule 
r 5 ^ (r) = -3AI2pn is valid within the n energy shell (sec Vrinccanu and Flannery (20(X)) 
for a detailed explanation). The perturbing potential (2) can then be written in terms of the 
components Ai and Aa as 


V(a) = -a—(A2sin<l> +A3C0S<I>) 

where the Stark parameter a = 3ZxanVnl2bv and a« and u, are the average electron-orbit 
radius and velocity. 

The Schrddinger equation for the time evolution operator U (f. fo) is 


(3) 

ot 

where f/o is the free atom Hamiltonian and V is the interaction potential (2). The position 
operator and hence the perturbation potential (2) commute with the unperturbed Hamiltonian, 
as one can prove directly from the matrix elements of the commutator {r. Ho] between any 
states within the energy £, sheU. The potential in the interaction representation 

is then i dfrtHeai with the potential in the Schrfidinger representation (V| = V). The equation 
to be solved, in the interaction representation, is 


ifi— as -a^(i 42 sin <b + A 3 cos ^)Ui 
dr dr 


(4) 
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where 6/i(r, to) = txp(iHot/fi)U ^xp{-iHoto/h). Because the set of operators {Li, Az, Aj) is 
closed under the commutation operation and generates the rotation group 50(3). the solution 
of equation (4) is obtained in terms of elements of this group as 

t/i(f.ro) =e'*^'/*exp(-i(<t - <t>o)(Li (5) 

Impact angles 4> and <t>o correspond to the position of the projectile at times t and 
respectively. The above solution (5) can be easily verified with the aid of the relations 

, cos A. + A3 sin k 
= /ijcosA. - Li sinA. 
e‘^^</^A3e“^^'''^ = A3 cos A.+ A2 sin A. 
which arc derived from the basic identity 

= B + + + (6) 

and the commutation relations (Li, Aj] = WA3, (Li, A3) * —WA: and [A2, A3] = 'AL\. 
The initial condition U(to, to) = 1 is automatically satisfied. Note, in the limit as a -► 0, that 
UUo. t) = l for all time and no (1, m) transitions can then occur. 

The transition amplitude for a Stark mixing process is 

aj,",’* (n^lf/i(oo,-oo)|na> 0) 

where the initial unperturbed sute |na). at f = - 00 . evolves to the final states In^), at r = 00 . 
Since a and ^ label the sutes within the same energy shell, the superscript n can be omitted 
and all dynamics is restricted to the energy shell defined by quantum number n. The full array 
of transition amplitudes is given by equation (5) in (7) and is feasible and efficient for practical 
numerical applications since it requires only matrix operations. 

The core of solution (5) is the exponential of the operator L\ - aAj. By using basic 
coitunutator algebra, Pauli's replacement and equation (6), this operator can be diagonalized 
as 

e-'w/*(L, _aA3)e‘«'* =yLi 

where g = (2p,/3)arctana and y = Vl +a^. The solution (5) has therefore the alternative 
form 

Ui(to h) = ( 8 ) 

which illustrates quite nicely how the action of the slow distant encounter charged projectile 
incident along the negative Z-axis can be decomposed into successive routions about the X- 
axis and alternating impulsive momentum transfers ±q{a) along the K-axis. In the limit of 
zero impulse q, (8) is unitary and no transitions occur. In the limit of small Stark parameter o 
the above solution (8), for undeflected collisions « (d> — <^o) = —rr. reduces to 

C/, =c“^^^+0{a^) 

in agreement with the impubivc result (Vrinceanu and Flannery 2000). 

It is however interesting to note that by introducing the Pauli replacement directly in the 
potential (2) and by vmiing the Rungc-Lenz vector as A = Af - N. the potential decomposes 
as 

V = + Vyv 

where Vj^ = -afAf: sin ♦ + Af 3 cos and = a{N 2 sin <t> + A /3 cos Because the 
commutators (M,. Nj], [Mi. Ho] and [A/. . Ho] = 0 (for any L j = 1.2,3 combination), the 
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Table 1. The four bases useful for desenbing the quanta! states of the hydrogen atom. 


Quantum Complete set of 

numbers commuting observables Ongin _ 


Orbital \nim)o Ho, Ly 
i>arabolic H\,Hi,Ly 


Algebraic 


Ho, Ay, Ly 
Ho, My, Ny 


Standard for spherical coordinates; describes correctly the sutes 
of the field-free atom 

Separation of Hamiltonian W = Hi //j in parabolic coordinates. 

tanv* yA'• n * «i ♦ "2 ♦ N1 1 
Parabolic basis: describes the Stark states for small electric 
field €, when the tnteracuon ^e£z is diagonal; q mm-ni 
The two rotation groups in which the dynamic symmetry group 
50(4) ■ 50(3) ® 50(3) decomposes; the equivalent angular 
momenoim for bo* 50(3) tepmenuiiooi ii ; - - l)/2; 


problem (3) becomes separable, in exactly the same way as the classical Stark miwng ^nations 
can be decoupled (Vrinceanu and Rannery 2000). The time evolution operator then factonzes 


U = Uh,UuUn 

where, of course. U„, = exp(-i//o(r - fo)/») and both and U, 

mUu/di = VmUu and ihdUN/dt = V^^/^. similar to equaaon (3). With the aid of the 
group theoretical result (5). the solutions for the operators Uu and Un are then 

exp[— i/h(^ — •ho)(Wi “ «Af 3 )]e (10) 


and 

= e'*"'/* exp[-i/h(d) - ♦o)(A'i + aNj)!*"**'"’'*- 
In calculating the ampUtude (7). four interesting basis sets can be chosen fw the one- 
electron hydrogen-Uke atom. Table I summarizes key properties of these bases, -^e orbiW 
basis is useful for describing the field-free atom, before and after the collision, wherw ^ 
algebraic basis appears naturally as a basis where Mj and Nj are diagonal. The solution (9) has 
the simplest expression in this algebraic basis. All four bases in table 1 span^ n degenerate 
energy shell and can be equally adopted to characterize the hydrogen atom. The algebraic basis 
spans a tensorial product of two spaces (Im> ® IW) corresponding to spaces used for matnx 
representation of the product S0(3) 0 50(3). The two spaces have the same dimension 
because - («^ - l)h* and are associated with two angular momenta 

with j ={n — 1)/2. . . u j . _ 

The transition amp litude between the two algebraic stetes is then the product a(<.v).(>»») - 
of two ampUtudes for A# and Af independent actions. Each factor is the matrix eleincnt 
of a ;■ = (n - l)/2-dimensional representation of the routions represented by equations (10) 
and (11). For example, from (10). one obtains = F(a)|i'(*< 

F(o) = (1.0.0)1P<^> (i.O. (1.0.0)) 

where PW.fm.nj.as)] » the Wigner matrix representation 

angle ^ about direction n (see Louck (1996) for the explicit expression). AO is the polw 
angle O - O© swept out during the interaction. The transition probability in the space of N 

is the element v'w of the matrix F(-a). 

Calculation of the transition probability between orbital states requires the exphcit umti^ 
transformation between the orbital and algebraic bases. This can be obtained by direct scalar 
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products of the orbital and parabolic states for which explicit coordinate representations 
are known. The result may be written in terms of hypcrgeometric functions (Tarter 1970). 
However, an equivalent result is provided by the 50(4) ^ 50(3) © 50(3) isomorphism. The 
orbital state, as a combination of two angular momentum states, is 

j 

\nim) = ^ 

where the trar*”formation matrix C^^***^ is given by the standard Clcbsch-Gordan coefficients 
{j^J The transition probability for the Im tm' transition then becomes 

= Yi (12) 

which can-be expressed in matrix form as 

= Trace[C'''"V '■(«)] 

where is the transpose of matrix C. This result (12) is in exact agreement with the solution 
obtained by Kazansky and Ostrovsky (1996a, b). who used the rotating-framc approach. 

The quanta! development above is exquisite in that it follows exactly the same reasoning 
within the exact classical mechanics solution (Vrinceanu and Ranncty 2000). This result 
exhibits the essential power of the 5 O (4) symmetry group for the energy shell of the hydrogen 
atom. The common 50(4) symmetry therefore transcends the chosen formulation (classical 
or quantal) and provides a classical-quantal correspondence at a level more fundamental than 
Ehrenfest’s theorem and the Heisenberg correspondence. 

In practice, the fourfold summation (12) and the use of the Wigner rotation matrices V 
in F{a) arc not very efficient and the difficulty of calculation increases dramatically with n. 
Instead, the solution (5) provides a much simpler approach, because the matrix elements of the 
argument in the exponential have simple expressions directly in the orbital basis. The array of 
transitions is obtained at once, all within one matrix exponentiation of a band diagonal matrix 
for which efficient algorithms are available (sec e.g. (jolub and Loan 1983). 

When the projection of the initial and final angular momentum is not determined, the 
transition probability is 

P#'i(a) = —— ^ ^ \avm\tm\^- 

‘ m—t rn'm^v 

The present theory is now applied in the following paragraphs for low n = 2 and 3, 
when analytic results can be derived, and for n = 28, when accurate numerical results can 
be obtained. A matrix representation for the operator L\ - ais required. Instead of the 
spherical basis \tm), which is difficult to use in this case, we define a new linear basis obtained 
by mapping the (/, m) quantum numben to a unique index lc = ^^+/ + m + l,in such a way 
that (0,0) I, (l.-l) -► 2, (1, 0) 3,(1,1) -► 4, (2, -2) 5 and so on. The inverse 

mapping is given by ^ = floor{\/F—T) and m =1: — — f — 1. The index k counts the 

degeneracy of the energy shell, and runs from I to n^. The matrix element 

= >/(/ + m)(^ - m ^ \)/2SctSm\m^\ + >/(f - m){l + m + \ )/2Stt^m\m^\ 

of LI is non-zero only for = 0 and Am = ±1. which reflects the fact that the cylindrical 
symmetry of the Rydberg atom is broken by the precession of L about the field of the projectile. 
These m-changing transitions are however conditioned by the full structure of solution (5) 
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Figure 1. Matrix rqjrcscntation of Li — aAj for n * 3. 


which shows that such transitions are only in evidence for non-zero a. The matrix element 




{O - m2){ni - 


( 2 f + l)( 2 f- 1 ) 




-/ 


lit 

(2/ + 3)(2/ + l) 


Scui^m 


of the component A 3 = -(2/3rt)z along the fixed Z-axis of quantization is non-zero for 
Af = ±1 and Am = 0 transitions. These dipole transitions only contribute for non-zero a. 
The matrix L\ - a A 3 has therefore the band diagonal structure, as illustrated in figure 1 for 
the special case of n = 3. The transition amplitude for transition k Id \s\htkt matrix 
element of the exponential of the matrix —i A<>(L 1 — or A 3 ), sandwiched between the rotations 
implied by (3). When a % 0, the dipole-forbidden transitions are not possible because the 
cransition matrix *5 exp -iA t> Li still has a band diagonal structure. When a increases, more 
and more off-diagonal clem become populated, resulting in an increasing number of dipole- 
forbidden transitions with ... = ff and Am = m' - m. Efficient algorithms, using Padi 
approximations, are available (Golub and Loan 1983) for matrix exponentiation. The whole 
array of transition probabilities for ail t and ^'is then obtained at once. 

Analytical probabilities for / f' transitions within the n = 2 shell arc listed below: 


n +a^cos(yA<>)]^ 

P® = ^^^t2 + a^+a*cos(yA4>)lsin* 

= (^) cos(x A<>) + cos{ 2 y A^)]. 

The following transition probabilities for the /i = 3 shell are also obtained: 




— 2a^ + + 8 a^ cos(y A <b) + 2a^ cos(2y A<b) 
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FlgUK 2. Stark mixing transition probabiUucs for the 28f -► 28f' transition array, with a = 0.2 
(left) and a ac 0.8 (right). 


[1 +o^cos(yA<>)]^[2 + o^+a^cos(yA<fr)]sin^ 

[2 + a^ + a^cos(KA'J>)l^sin^ 

pW _ t(6 + 30a^ + 4a‘ + 3a*) + 8a*(3 - 2a^ + 2a*) cos(y A<1>) 

+2a*(l 1 - 2a* + a*) cos(2y A<>) + 8a* cos(3y A4>) + a* cos(4y A<^)] 
p(5> = [(20 + + 32a* + 8a‘) + a*(26 + 20a* + 9a*) cos(y A<t>) 

+2o*(4 + a*) cos(2yA<t) + a* cos(3y A^)] sin* ^ 

p^(5> = [(90+ 240o* +318a* + 196a* -63a') 

+8a*(15 + 22a* + 14a*) cos(v 14)) + 2a*(23 + 22a* + 13a*) cos(2y A<l>) 
+8a^ cos(3>'A<>) +a* cos(4>'A<l>)]. 

The detailed balance relation 

is satisfied by the present treatment, so the P/,"’ probabilities for transitions with £' > f can be 
obtained. 

Figure 2 displays the results for calculation of the transition probabilities inside the « = 28 
energy shell. An undeflected path (A<t> = -n) is assumed. For small or = 0.2 only elastic 
or transitions with small angular momentum transfer At have significant probabilities; a band 
along the diagonal is exhibited. As t increases. At increases and then decreases as ^ n - 1. 
As ot increases, the band broadens and larger angular momentum transfers become possible 
foralU. 

The present treatment is valid (a) for weak fields, in evidence for impact parameters 
b > b* = (3Zi/2)*^^a«, which implies Stark parameters o < a* = (3Zi/2)*/^(Un/v), and 
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(b) for adiabatic collisions when the collision frequency <i> is less than the orbital frequency 
of the Rydberg electron, so the Pauli replacement holds. These two conditions 
combine to yield the partitioning u < u* = (3Zi/2)^^^Un and b > b* in (u, i>)-space for Stark 
mixing collisions (Flannery and Vrinceanu 1998) for a slow encounter in a dipole held. The 
limit V < V* defines our meaning of uitralow collision energies. The cross section for Stark 
mixing is 


. 2, f fSW = 2. 


When y < u* and < b*, the Static parameter a > 1. Since the transition probabilities are 
bounded for large a, the contribution to the a-integration is vanishingly small for large a, 
decreasing as and can, in practice, be neglected for a > 1. This implies that the lower 
limit to"^ is At uitralow energies, this limit is always much greater than the 

weak-field limit br At the upper limit v = v* of uitralow energies, b, approaches from above 
the weak-field limit b\ In practical calculations of (13), various physical effects such as Debye 
screening in a plasma, quantum defects and spin-orbit coupling deicmunc an upper limit to 
b and hence a lower limit amm to a. For trajectories with zero deflection {Ad> = —tt), (13) 
varies universally as (Zifl«u«/u)^ Departure from this variation is governed by A4>(i>, v) and 
by the physical limits imposed upon the o-integration. 

In summary, this letter has presented a new form (5) of the exact quantal solution for 
the Stark mixing probabiUtics. Based on this new solution for the whole array of transitions, 
analytical expressions for small quantum numbers n and accurate numerical results even for 
large n can be obtained. Using the rich symmetry group of the hydrogen atom, the relation 
with previously published results (Kazansky and Ostrovsky 1996a, b) has been developed in 
an exquisite fashion. The symmetry group also provides a complete exact classical solution 
(Vrinceanu and Flannery 2000). capable of explaining the main features of the quantal results 
and of providing a quantal-classical (dynamical) correspondence at a level, more fundamental 
than Ehrenfest’s theorem and the Heisenberg correspondence. 


This research is supported by grants from AFOSR: F 49620-99-1-0277 and NSF: 98-02622. 
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Classical Stark Mixing at Ultralow Collision Energies 

D. Vrinceanu and M. R. Flannery 
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Exact solutions of the time-dependent classical equations are obtained for the full array of angular 
momentum mixing transitions ni ^ ni' in atomic hydrogen induced by collisions with charged particles 
at ultralow energies. A novel classical expression for the transition probability Pvt is presented. The 
exact classical results for Pvt(a) as a function of C i* and the Stark parameter agree exceptionally well 
with (exact) quantal results. They complement the quantal results, by revealing essential charactensties 
which remain obscured in the quanta! treatment. 

PACS numbers; 34.50.Pi, 34. lO.+x. 34.60.+z 


Stark mixing occurs-when the electron of a Rydberg 
atom (in a state with principal quantum number n) changes 
its angular momentum, without changing its energy, as a 
result of a collision, at large impact parameter b, with a 
slow massive particle of charge Z\e moving with velocity 
V. It is a subject of broad interest and importance in many 
areas of modem physics [ 1) and chemistry [2], astrophysics 
[3], line broadening [4], Auger processes [5], and for anti¬ 
hydrogen formation by thrcc-b^y recombination [6,7] at 
ultracold temperatures. Although remarkable effort has 
been devoted to obtaining theoretical solutions for Stark 
mixing in Rydberg atoms to various levels of approxima¬ 
tion [1-13], the purpose of this Letter is to point out that 
the problem is capable of an exact solution in the classical 
formulation. The exceptionally rich dynamical SO{4) sym¬ 
metry of H(n, () is the key foundation which allows both 
classical and quantal exact solutions to be constructed (14] 
in a similar and unified way. There is substantial renewed 
[1,15-19] interest in the power of classical dynamics in 
almost all fields of modem physics, attributed to the de¬ 
sire [15,19] to obtain a more thorough understanding of 
the classical-quantal correspondence. Stark mixing by ion 
impact is probably the last problem in collision physics 
which is capable of an exact solution. 

The present new treatment is not an extension of any 
previous theory and is capable of providing the first com¬ 
prehensive classical solution for the full array 
of collisional transitions in A new expression for 

the classical transition probability is defined in a lan¬ 
guage which exploits the dynamical synunctry. The de¬ 
rived probability for the general array € —► of transitions 

has a very simple functional form, can be easily calculated 
for any principal quantum number, and provides physical 
insight and simple geometrical explanations for the behav¬ 
ior of the transition probabilities. Stark mixing probabili¬ 
ties are calculated and compared with the (exact) quantal 
results [14,20]. Considering the Rydberg atom in a frame 
which rotates together with the intemuclear axis, the Stark 
mixing problem can in principle be reduced [11,12] to the 
problem of the Rydberg atom in mixed static electric and 
magnetic fields. This approach is successful only for the 


particular case of f = 0, which is fully recovered by the 
present general fixed frame formulation. 

The target Rydberg atom (with averaged electron orbit 
radius a^, velocity i;„, momentum pn. and angular fre¬ 
quency (t)n * Vn/Cn) is Centered at the origin 0 of a fixed 
coordinate frame. The trajectory of the projectile, initially 
moving with impact parameter b along the positive Z 
direction, is assumed to be confined in the YOZ plane. In 
addition to the energy £, constant along the Hamiltonian 
Ho * P^I2me “ e^/r, the angular momentum L “ 
r X p of the unperturbed Rydberg electron and the 
Runge-Lcnz (or eccentricity) vector. 


A 



XL-" 



( 1 ) 


directed toward the pericenter and normalized to angular 
momentum units, arc also conserved. These quantities de¬ 
fine the dynamic SO(4) symmetry of the hydrogen atom. 
Because die SO(4) group is isomorphic with the direct 
product SO(3) ® SO(3) of two rotation groups, a special 
decomposition, L * ivf + N and A * M “ N, permits 
the dynamics of th' liydrogcn atom to be separated into 
two decoupled motions. The generators M and N act in¬ 
dependently as angular momenta and arc also conserved 
quantities for the unperturbed Rydberg atom. They evolve 
independently [2l]^with time on application of a con¬ 
stant electric field i and preccss about £ with the Stark 
frequencies 0)$ “ ^{3/2)anVn(i/e). For weak fields, 
the Stark splitting ^Es * bcjs < * Afn, the (n —► 

n ± 1) energy splitting. For constant fields, the vectors L 
and A vary periodically with frequency tog* 

The weak field approximation assumes that the time- 
dependent electric field £(r) generated by the passing 
projectile of charge Zie is constant over the atom’s spatial 
extent and the dipole approximation for interaction poten¬ 
tial V — er • £ is valid. Hence 


Tit) 


Z\ek ^ Z^ ^ 
^ vb dt ^ 


where <t> is the polar angle between the intemuclear vector 
R and the Z axis. 
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On assuming that the collision is orbital adiabatic (<I) < 
o)„), i is constant over one period so that the slow rate of 
variation AL/T of the angular momentum over one orbital 
period T = lit/(tin, is the classical average 

— = f^(r X i)dt = -e(r> X i(t). 
dt T Jo 

Since the weak field approximation (tuj < (■>„) also holds, 
the vectors L and A change very little over one orbital pe¬ 
riod and <r) » -3A/2p„ = -3 Ma„v„/2e^). The fol¬ 
lowing set [3], 

- = —tJs^ ^ A, ~y~ “ ~(tisR ^ L, 

dt dt 

of coupled equations^m then be deduced, where both R 
and tt»s = ai now vary with time. The Stark parameter 
a is 3Zia,i/„/2hv. The weak field (tus < w„) and adia¬ 


batic (<i> < v„/a„) conditions combine to yield the par¬ 
titioning b 2 (v/vn)an, when v 2 u* (3Zi/2)*'*v,, 
and b 2 (3Zi/2)'^^a,. when v v*. of {v,b) space for 
validity of the present solutions. For u 2 u’, or remains 
<1. For V < u*. a can exceed unity. When written in 
terms of vectors M and N, the above set of differential 
equations yields the set of decoupled equations, 

^--ai?XM. ^ = (2) 

d(p 

Since the magnitudes ^ = (L^ + A*)/4 = 

n^ti-/4 remain constant throughout the collision, exact 
solutions of Eqs. (2) can then be obtained (14) at general 
angle 4* in terms of finite rotations from the initial 
values M(<l>o) and N(<l>o). via the orthogonal trans¬ 
formations M(4>) “ and N(4>) = 

f/Ar(<t>.<l>o)N(<l>o). For example, if the initial state is 
specified by the vectors (L, A) at <I>o “ ir, then the final 
I state (L', A'), at - 0, is determined by the rule 


- y-2[i -I- a^cos(irr)]Li + ay'* sin(7ry)A2 + ay'^[l - cos(iry)]A3. 
L '2 - -cos(iry)l 2 - r"‘sin(7rr)l3 + ay"‘sin(iry)Ai. 

L'j - y"‘sin(iry)L 2 - y~\a^ + cos(iry)3L3 + ay"^[cos(iry) - l]Ai. 

- y-2[l a2cos(fl-y)]Ai + ay"‘sin(rry)L2 + “ cos(iry)]L3. 

A 2 ” — cos(7ry)A2 ~ y”* sin(7ry)A3 + ay ' sin('iry)Li, 

A '2 - y"'sin( 7 ry)A 2 - y"^[o^ + cos(iry)]A 3 + oy"Hcos(iry) - l]Li. 


Here y * Vl + and the components of the initial and | 
final vectors are defined in the fixed coordinate frame con¬ 
sidered. The above exact solutions (3) are easily verified 
and satisfy the invariant relations 

L' • A' - L • A - 0 (4) 

and 

+ A^ - (5) 


The orbit of the final state (n. L') is confined to a plane 
perpendicular to the final V and the energy is preserved 
(since n does not change). The above constraints, Eqs. (4) 
and (5). define in the (L. A) space a hyoereurface on which 
the in^'ial state, defined by the init’.il angular momentum 
quanm^i number €, is uniformly di'Ltbuted. The volume 
of ihTs hypersurface is therefore given by 


V«< - jj fi(|L| - ne)S(\M - fiy/n^ - f2)5(L • A)dLdA. 


( 6 ) 


Each point from this manifold evolves during the coUision according to the rule (3). so that only a fracUOT of ^ssible 
initial states can produce the final state with angular momentum quantum number (, after the collision. The volume ot 
(L, A) space which overliqis the initial and final states is 


yntr - jf 5(|L| - f«)5(|L'| - €'«)5(|A| - W/i^ - €2)5(L • A)dLdA. 


(7) 


The I-* V transition probability is then, in a geometric 
sense, the ratio of the two volumes 






( 8 ) 


This is a novel result in that Eq. (8) is defined in terms 
of the new (L. A) representative space, being more appro¬ 
priate. than the customary (r.p) phase space, for direct 
expression of the dynamical SO(4) symmetry of 


The six-dimensional integral (6) can be calculated directly, 
while the integral (7) eventually reduces to yield the one¬ 
dimensional integral [14] 

(„) _ 2 ^ t ■_ { 

irfirfi 1 - cosAf } yj(z^ - A) {3 - z^) 

for the transition probability. This, in turn, may be 
expressed in terms of the complete elliptic integral. 
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no. 1. Probability for the 18 — C transition, within the 
n = 28 energy shell, for a given Stark parameter a - 0.4. 
Exact quantal results are denoted by dots. 


K{m) = “ m^in^x) '^'^t/x, as 






f 

0 . 

Xf(3-.a)/3] 


3 < 0 , 

3 > 0. J!l < 0, 
3 > 0, Jil > 0, 
(9) 


where 

Miln.ein^x) 

3{e/nJ'/n’x) 


cos(iii + U 2 ) ^ cos;^:^ 

1 - cosx 

cos(«i - U2) - cosx 
1 cos^ 


( 10 ) 


The angles u\ and U 2 depend only on the initial and final 
states viacosui = - 1 andcosu 2 * 2€'^/n^ - 1. 

The rotation angle x^ which depends only on the Stark pa¬ 
rameter a and the polar angle A4> * 4^0 swept out 

during the collision time interval {to, t), is determined by 


cosy = [U+ a^cos(>/l + a2A4>)]/(l + a^). 



FIG. 2. Density plots of the f — T transition probabilities, 
calculated within the quantal treatment for a -• 0.4 and n — 28. 
The probabilities increase as the color becomes darker. The 
continuous and broken lines represent the classical ^ * 0 and 
3—0 ridges. 


The condition 3 < 0 defines the classical inaccessible 
region. Two types of singular points are apparent from the 
solution (9). At 3 * 0. the transition probability Pf({a) 
has a finite jump (step discontinuity) and at ^ = 0 it 
has a logarithmic (cusp) singularity. These features are 
displayed in Figs. I -5 for representative transitions. For a 
given Stark parameter a and initial angular momentum 
the solutions of the equations = 0 and 3 = 0 provide 
the critical values 



of the final angular momentum, where the cusp or step 
singularities are located. The inequality < 3 is strictly 
fulfilled except for the limiting cases of f or — 0 or 
n when - 3. These limits are readily deduced from 
Eq. (9) to provide the following probabilities: 


sin(Ar/2)-y/sin^(Ar/2) “ (C/n)^ ' 
• sia{x/2)^J{('/n)^ - sitt^(x/2) ’ 

p(")/ \ ^ _ 

sia{x/2)^J(€/ny - sm^(x/2) 


The 0 transitions have zero classical probability. 

The specific case [12] of 0 —» €' transitions is therefore 
directly recovered from our general result (9). 

Figure 1 shows the classical probability for the 
representative array f — 18 —► f' transitions. The quantal 
results [14,20] oscillate about the classical background. 
Figure I also displays the 0 cusp at ffollowed by 
the 3 “ 0 downward step at ('+, as (' is increased from 
0 to n. In the central region .A < 0 and 3 > 0. 

Figure 2 provides the density map of the quantai r'oba- 
bilities [14,20] for the full array of € —► f' ^ansitions 
at a — 0.4. In the classical forbidden regions (where 
3 < 0). the upper left and lower right comers, the quan¬ 
tal probabilities decrease exponentially. The lines of singu¬ 
larities .A “ 0 and 3 — 0 are also shown. Along these. 



no. 3. Probability for the 18 -• 2 transition, within the 
n - 28 energy shell, as a function of the Stark parameter a. 
Exact quantal results: dotted line. 
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0.0 0.2 0.4 0.6 0.8 1.0 


a 

FIG. 4. Density piotToLthe f - 18 -» transition probabili¬ 
ties. calculated within the quantal treatment for n •• 28, as a 
function of a. 

the classical transition probabilities display cusp and step¬ 
like behavior, respectively, and. in their vicinity, the quan¬ 
tal results have local maxima. Figure 1 follows from a 
vertical line drawn through the plot at € * 18. show¬ 
ing a cusp-step variation. Figure 2 predicts step-step and 
step-cusp variations for transitions from € > 19. 

In Fig. 3, the classical and quantal transition probabili¬ 
ties for the 18 —• 2 transition are plotted as a function of 
the Stark parameter a. This plot can also be obtained by 
following the variation of the transition probability along 
a horizontal line with * 2 in Fig. 4. 

Figure presents a density map for the quantal proba¬ 
bilities for transition from the initial angular momentum 
f = 18 to any f' and for any value of ot. As a 0 the 
span of possible final angular momenta is reduced, such 
that only elastic transitions are possible at a ■“ 0. The 
classical singularity lines A ^ 0 and S “ 0 illustrate 
again the correspondence with the quanta! results. The 
cusp-step pattern of Fig. 1 is also explained by a vertical 
line at a « 0.4 of Fig. 4. For low a < 0.4, a step-step 
variation is predicted, i.e., the accessible (' lie within the 
range f'_ < (' < ('+. The results in Fig. 5 confirm this 
prediction. Cusp-cusp patterns occur at higher a 2 0.5. 

In conclusion, the exact solution (3) of the classical 
equations (2) has been obtained, by exploiting the SO(4) 
dynamical symmetry of The novel expression 

(8) provides the general result (9) for the classical proba¬ 
bility Pu- for the full array of Stark mixing transitions 
nf — nC. Since and 2 are symmetrical in (€,€') 
the probability (9) satisfies the detailed balance relation 
- 2('P^tt'- Probability conservation fl P^^lhdC 
=« 1 is also satisfied. Cross sections [14] may be calcu¬ 
lated from Eq. (9). Exceptional agreement is obtained be¬ 
tween the classical and quantal /*«<(») as a function of 
(, (' and the Staric parameter a. The common SO(4) sym¬ 
metry provides this classical-quantal correspondence at a 
level more fundamental than Ehtenfest’s theorem and the 



Heisenberg correspondence. The classical method is also 
complementary in that it reveals very succintly essential 
and valuable characteristics which remain obscured within 
the quantal treatment This reflects the essential power of 
classical dynamics when applied to collision problems. 

This work has been supported by AFOSR Grant 
No. 49620-99-1-0277 and NSF Grant No. 98-02622. 
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Abstract 

An exact expression for the probability of Stark mixing transitions between 
arbitrary angular momentum states t and V within the same energy shell n* as 
a result of a collision with a slow charged projectile at large impact parameters 
is presented. The formula obtained is compact and easy to use for numerical 
evaluations even for very large quantum numbers (/i 100). A classical 

approximation is directly obtained and compared with the exact quantal result 
in the limit of large n. Two distinct sets of quantal oscillations are predicted. 

(Some figures in this article are in colour only in the electronic venioo; see m. lop. ore) 


In a previous letter, Vrinceanu and Flannery (2000a) provided a new quantal theory in a fixed 
reference frame for the probabilities for the full array nt ^ nV of transitions induced by 
collisions between H(/if). or any Rydberg atom, and an ion of charge Ze at ultralow energies. 
The conventional approach to this problem consists in solving a coupled-channel system of 
differential equations. For large n, this task becomes prohibitively difficult even when some 
..simplification can be introduced, such as that, for example, by Beigman and Syrkin (1985). A 
time-evolution operator method (Pfennig 1971) depends on the roots of an algebraic equation 
of degree so that, in general, it is not possible to obtain analytical results and it was only 
applied to transitions between the parabolic states of the n = 2 shell. The alternative method 
proposed by Vrinceanu and Flannery (2000a) exploits the rich dynamical SO(4) symmetry 
of H(nf) to solve the problem exactly in the fixed-frame representation. A similar approach 
by Kazansky and Ostrovsky (1996a, b) was successful only for the special case of nO ni* 
transitions. The solution presented by Vrinceanu and Flannery (2000a) has the advantage 
that the full array P^i of transition probabilities is obtained at once for a given velocity and 
impact parameter of the projectile. Moreover, analytical expressions were obtained for low n. 
Numerical results were possible for low to intermediate n ^ 30. 

The present letter reports a significant advance to what was previously obtained. An 
exact analytical expression for the quantal probabilities P^?^ is now presented for any n. 
The remarkable simplicity of the result allows both accurate numerical results for very large 
n ^ too, and an immediate classical limit 

The target Rydberg atom is described by the Pauli-Lenz vector operator A and the angular 
momentum vector operator I. The projectile with velocity v along the z-axis, intemuclear 
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vector R and impact parameter b is characterized by the Stark parameter a = 3Za„v„t2bv, 
where a, and v„ are the average electron-orbit radius and velocity. The adiabatic dipole 
interaction is V{R) = -Ze^(r) • R/R^, where the averaged vector position of the Rydberg 
electron of mass m, is (r) = -3A/{2m,v„). Under the adiabatic, dipole and classical path 
assumptions, the SchrOdinger equation for the evolution operator in the interaction picture is 
(Vrinceanu and Flannery 2000a) 

= -o((A2sin<b + AjCOsd>)f// 

34 > 

when time t is replaced, via the relation = -vb. by the angle <t> between R and the 
(positive) z-axis of a fixed frame. The previous work of Vrinceanu and Flannery (2000a) used 
the exact solution for U, directly to provide the full array of transition probabilities. However, 
uponihedecomposition L = M + N and i4 = M-N. the solution separates as 

Ui = Uu® Vh 

where the operators Vm and Vn are defined by 


and 

and where A<l>(r) = <l>(f) - ♦(<&)«the polar angle swept out by R within the time inteival 
(2 _ tg). Since the angular momentum-like operators M and JV commute, the corresponding 
evolution operators C/m and Uw act independently as rotations in carrier spaces of dimension 
2; + 1 = /I. 

The probability for transition <-*•<' between states with given angular momentum is 
defined by 


p("> _ * 

"2f+l 


E E I-?---!' 




(I) 


where, the (i, m) -► transition amplitude between angular momentum states within 

the energy shell of quantum number n is 


and where » Wigner’s (2J + l)-dimensional matrix representation of the finite rotation 

The tenfold summation within equation (1) can be contracted to yield the following result: 


E / 

" z-ir'-n I ■' •' 



(3) 


which is compact and involves only one summation. Here {• • •) is the 6-y symbol for coupling 
of three angular roomentt and Hjl is a special matrix element of the irreducible representation 
of the group 0(4). This function (also called the generalized character xi . associated with the 
irreducible representation of the rotation group) is well studied (for example, see the books 
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by Talman (1968) and VarshaJovich etal (1988)). It can be written in tenns of ultraspherical 
polynomials as 


HjlW = (2L)!!, 


/ (2>-k1)(27-L)! 
' ( 2 ; + f,+ l)! 




The angle n is the proper angle of the rotation produced by the successive rotations f/y' and 
Um, and is uetermined from the two-dimensional representation of the product by 

n 1 r n 1 +a^cos>/l +a^A<J> 

cos| = iTrace:[f/«i/^-'] =--• W 

The formula given by equation (3) is capable of providing exact numerical results even for 
large quantum numbers n, in contrast to equation (1) where the number of terms to be summed 
increaMS.dramatically with n. Moreover, a classical limit of equation (3) can be directly and 
easily obtained. Note that the capacity of the projectile to produce angular momentum changes 
is governed solely by the key parameter ij. the argument of Hji, in equation (3). Since A <b varies 
with time f. expression (3) furnishes “ * function of r and the Staric par^ter o. 

within parameter tj. It can be shown that both detailed balance (2t + DP^."* = (2f' + 
and probability conservation Ze satisfied. For the special case of zero initial 

angular momentum, the 6- j symbol is (Varshalovich et al 1988, p 299) 


The probability for 0 ->• f' transitions is then 
f’r"o(x) = 

in agreement with Kazansky and Ostrovsky (1996a,b). Equation (3) may be recast into the 
interesting form 


P<-’(X)=n(2/' + l) ^ 


profx). 


Since the square of the 6-; symbol is the probability of coupling three angular momenta, a 
physical int ->pretation may be given for the above formula. The Stark mixing transition can 
therefore be" ;en as a multi-step process involving partial waves of angular momentum L. 

The proof of equation (3) is based on standard analytical techniques of rotation algebra. 
The main steps are represented by the diagrams in figure 1. The conventions adopted are 
similar to those described in detail by Zate (1988) and by Varshalovich etal (1988). Only the 
basic rules are explained here. Every internal line represents a summation over the magnetic 
quantum number of an angular momentum with the quantum number attached to the line. A 
double line represents the summation (2Z>(2y + 1).. ■) over the angular momentum quantum 
numbers. The coupling of two angular momenu is represented by a vertex with two incoming 
and one outgoing lines and has associated Clebsch-Gordan coefficient 
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(d) 


Since the Clebsch-<iof 'an coefficient is real, the conjugate symbol with one incoming 
and two outgoing lines also represents the same coefficient. The mim 2 matrix element of the 
2y + 1-representation of rotation R is represented by a directed line with a wavy line attached 


* Un>i\R\jmi) = 


> 


;m, yiBj 


The matrix element associated with the inverse rotation R~' has the same representation 
but with the direction reversed. Under this simple notation, the amplitude equation (2), for 
example, has the following representttion: 
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Table 1. Exact expressions for the Siaric mixing transition probabilities inside the shell n = 
1.2. 3, 4 and 5. 


n 

e 


/ 

(2/♦!)/>;." '(I) » / * (1 -i-)''*'’' K ■ (i = cos 17/2) 

1 

0 

0 

1 

1 

2 

0 

0 

1 



0 

1 

1 

1 


1 

1 

I 


3 

0 

0 

1 

9 

(-l+4x*)> 


0 

1 

s 

J 



0 

2 

t 

9 

1 


1 

1 

1 

1 


1 

2 

} 

3 + 2x* 


2 

2 

9 

l9 + 22jr'+4x* 

4 

0 

0 

1 

(-x+2x*)* 


0 

1 

1 

5 

(-l+6x^)> 


0 

2 

1 

4x^ 


0 

3 

4 

J 

1 


1 

1 

1 

B 

2+ I45x^ -396x^ + 324x‘ 


1 

2 

4 

5 

l+9x^ 


I 

3 

il 

4 + 3x^ 


2 

2 

1 

2 - 5x^ ’♦•4X* •►4x* 


2 

3 

i 

11+20x2 + 4x^ 


3 

3 

B 

52 + 75x^4*44x* + 4x‘ 

5 

0 

0 

i 

(l-ilx^+lbxV 


0 

1 

t 

B 

(-3x + 8x5)^ 


0 

2 

B 

(-1^. 4x2)2 


0 

3 


x2 


0 

4 

« 

I 


1 

1 

i 

9 - 162x2 ^ 99^4 _ j792x9 + i024x» 


1 

2 

h 

3+114x2- 384 x^ 4.512x* 


1 

3 

16 

B 

1 4*4x2 


1 

4 

TO 

54*4x2 


2 

2 

9 

27 4* 2 .i 4 u 2 .284x* - 76&X* 4* 1024x* ' 


2 

3 


4 - I3x2 + 2ax*4*16x‘ 


2 

4 


17 4* 38x2 4* 8x* 


3 

3 

is 

59- 52x2 _ ig4x‘4.288x‘4*64x* 


3 

4 

ra 

95 4* 24Ax2 4* |68x* 4- I6x* 


4 

4 

A 

2509 4* 3644x2 4* 3624X* 4* 1184x® 4* 64x* 


The main steps (aHd) in the proof of equation (3), illustrated in figure 1. are outlined 
briefly below. 

(a) The upper left-hand diagram is just the u^sition probability equation (1). The next 

diagram is obtained by inserting the identity operation on line f'. 

(b) The rotation matrix element from one branch of the Clebsch-Gordan coefficient can be 
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Figure 2. Stark mixing probability for tnositioas within the energy shell« * 100 from the initial 
state ^ » 80, with angle ff given by cos n/2 * 0.1. as a ftinction of the final angular momentum 
The exact quanta! result is represented by the dots, the classical limit is represented by the full 
curve and the initial state is indi c ated by the dotted line. 

* redistributed* to the other two branches. As an effect of this, ^ matrix elements of 
the rotations f/v on branches with angular momentum j are annihilated and the Ums are 
multiplied by f/jj* ^ rotation C = The angle of this rotation is r; of 

equation (4). 

(c) This step involves the identity provided by Varshalovich et al (1988) in equation (22), 
p 260. This formula gives the summation of products of four aebsch-Gordan coefficient 
in terms of a double summation over the angular quantum number and the magnetic 
quantum number of two Clebsch-Gordan coefficients multiplied by two coefficients 
represented as tetrahedrons. 

(d) It is easy to prove that the central object, which has two loops linked by the L-summaiion 
(and implicitly an M summation), is a rotational invariant. Therefore, one can choose the 
direction of the z-axis along the direction of the rotation C. As a rotation about the z-axis, 
the matrix representation of G is diagonal and the line that links the two loops disappears, 
since only terms with =: 0 are allowed in the summat'on represented by this line. The 
transition probabi!u/ is then obtained as an L-summ- -ion of products of two identical 6~j 
symbols and two objects which depend only on the ?n^c r;, the principal quantum number 
n and the angular momentum index L. These objects are identical to HjUn) b«:ausc, by 
definition. 

Ml 

A analytical proof of equation (3) has also been derived without utilizing 

the diagrammatic techniques presented here. The transition probability equation (3) is a 
polynomial in the variable x = cosij/2. Only even powers occur in this polynomial and, 
the £>-summation in equation (3) starts with |f — 1'\, this pol)momial has the factor 
(1 -x*)'*"*’'. Table 1 provides the exact analytical probabilities for transitions within the states 
with principal quantum numbern = 1.2.3,4 and 5. They agree with the expressions obtained 
previously by Vrinceanu and Flannery (2000a) for « = 2.3. Expressions for larger values of 
n have a similar structure and the corresponding polynomials can be easily calculated. 
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Figurt 3. Stark mixing probability for transitions within the energy shell n a 100 from the initial 
sute f a 60. with angle n given by cos r\ll a 0.9. as a function of the final angular momentum 
V, The exact quanta] result is represented by the dots, the classical linut is represented by the full 
curve and the initial state is indicated by the dotted line. 



Figure 4. Stark mixing probability for transitions within the energy shell .• a 100 from the initial 
state f a 5 to the final state a 80. as a function of cos r7/2. The exacbqti^^. * result is represented 
by the darker (wavy) curve. The classical limit is represented by the lighter curve. 


In the limii of very large quantum numbers the square of the 6-; symbol reduces (cf Zare, 
1988) to l/(24n'V'). In this case, the volume V of the tetrahedron formed by the angular 
momentum vectors with magnitudes t, t\ L and three vectors of magnitude j as indicated 
graphically in figure 1, is determined by 

>'• = Alt’ (2/V. 2}^ - en - - n‘\ 

The square of the generalized character function Hjiin) has the classical limit obtained by 
Kazansky and Ostrovsky (1996a) as 

j ^ 1 _1_ 

~ 2sinii/2 [sin2^/2 - 
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Using these two approximations and transforming the L-summation in equation (3) to an 
integration, the classical limit of the transition probability is obtained. This limit is identical 
to the probability derived by Vrinceanu and Flannery (2000b) from the classical equations. 
Figures 2-i demonstrate that the classical limit is effective for large quantum numbers. All 
the examples considered pertain to n = 100 and took a few minutes to produce on a PC. 
Figures 2 and 3. for small and large bv. respectively, illustrate that the quantal results oscillate 
about the classical background. Figure 3 further shows rhat the quantal result exhibiU the 
characteristic exponential decay in the classical forbidden regions {t < 20 and t' > 90). The 
classical singularities are related to the fact that the volu..Te V — 0 for certain combinauor- 
of angular momenta. A detailed analysis of the various possible classical panems is a subject 
explored elsewhere (Vrinceanu and Bannery 2000b). Figure 4 illustrates the variation of the 
probability for large momentum transfer with the angle ij. For large bv, when cos rj/2 -*■ 1. 
the transition is classically forbidden and the quantal probability decays exponenoally. Sm^l 
bv implies-cos t)/2 -*• 0 and the quantal results oscillate about the classical background. The 
present theory therefore predicts two distinrt sets of oscillations in /»,,{(i?)—one set evident 
in the variations with V (or 1) and the other in the variation with rj. i.e. with bv and AD. An 
approximation for the 6-j symbol and Hjl, more refined than the classical limits given here, 
should provide closer agreement with the quantal transition probabilities. 

The probabilities PiV are valid in the regions, (ti $ w* * (3Z|/2) Vn.b > b = 

(3Z,/2)‘/^a,]and[v ^v\b^bc = (u/ii,)n,].of(u.h)-space where both the adiabauc and 

weak field conditions are simultaneously satisfied. 

What is possibly the last remaining problem in collision physics, capable of an exact 
solution, has now been solved exactly in closed form. Not only does the present wlution 
reflect the intrinsic mathematical beauty of the problem, it can easily be appUed even for large 
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Classical and quantal collisional Stark mixing at ultralow energies 
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Exact classical and quantal solutions are presented for the full array of intrashell transitions n/’—n/'. 
between any angular momentum states, induced by slow distant collisions with a charged particle. The colli¬ 
sions are adiabatic with respect to the orbital frequency of the atomic electron and the transitions are induced 
by Che weak ion-atom dipole field generated by the ion moving along a classical path. The rich symmetry of the 
problem allows a unified approach and is the source of the excellent agreement, beyond the usual Ehrenfest's 
correspondence principle, between the classical and quantal treatments. A classical transition probability is 
defined. Probabilities for transition between any angular momentum states within a high Rydberg energy level 
are derived in exact analytic forms and are analyzed for a large number of numerical examples. The transition 
probabilities obtained from the three methods—quantal and classical formulations and Monte Carlo classical 
simulations—are directly compared to provide excellent agieemenL 
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I, INTRODUCTION 

The collision of a slow heavy charged panicle with an 
excited atom at large impact parameters induces transitions 
between neighboring angular momentum states of the exited 
atom. For very low velocity of the projectile, the transitions 
with change in principal quantum number are much less 
probable than the quasielastic angular momentum changing 
collisions at large impact parameters. Because these states 
are very close in energy, or are even degenerate as for hy¬ 
drogen and Rydberg atoms, the process is very efficient since 
little or no energy transfer is required. In fact, the cross sec¬ 
tions increase as the energy of the incoming particle is de¬ 
creased. This process is called Stark mixing and is important 
in many problems in atomic physics. For example, Bethe [1] 
analyzed the absorption of low-energy negative K" mesons 
in liquid hydrogen on the recognition that the Stark mixing is 
essential in such processes. Also, Stark mixing is included in 
the calculation of the Auger (or autoionization) process, 
which follows the collision between ions and atoms [2] and 
in zero-electron-kinetic-energy (ZEKE) photoeleccron spec¬ 
troscopy [3]. Stark mixing has also been important in astro¬ 
physics (e.g., [4]), in recent efforts to produce antihydrogen 
at 4 K [5] and for general three-body recombination at ul- 
u-alow energies [6]. The first stage in recombination at ul¬ 
tralow temperatures T, [6] is a very rapid collisional capture 
into high Rydberg states with high angular momentum and 
large radiative lifetimes, at a rate proportional to T,. 
Since the rt<hanging collisions are relatively unimportant at 
ultralow energies, the /-mixing collisions are essential in 
producing the low angular momentum states required to ra- 
diatively decay at relatively high rate to low n levels, thereby 
stabilizing the recombination. 

Experiments [7] on single ion collisions with alkali-metal 
Rydberg atom have measured large /-mixing cross sections 
for slow projectiles, including dipole-forbidden transitions. 
Various theoretical models have been developed to repro¬ 
duce the experimental data. Even though a set of coupled- 
channel equations can be written, their solution becomes im¬ 
practical for the large quantum numbers considered in the 


PACS numberfs); 34.50.Pi, 34.60.+Z, 34.10.+X 

experiment (/i-28). On averaging over the azimuthal quan¬ 
tum number m, the size of the problem becomes much re¬ 
duced and satisfactory results have been obtained [8]. Con¬ 
sidering the Rydberg atom in a frame that rotates together 
with the intemuclear axis, the Stark mixing problem is re¬ 
duced to the problem of the Rydberg atom in mixed static 
fields: electric, provided by the projectile ion, and magnetic, 
produced by the noninertial (Coriolis) forces. In this way, the 
well known results (in both classical [9] and quantum 
[10,11] mechanics) for the problem of interaction between 
weak fields and an atom can be adopted [12] to provide a 
solution for the Static mixing problem. In a remarkable series 
of papers, both classical [13-15] and quantal [15,16] ver¬ 
sions of this approach have been successfully applied for the 
zero to higher angular momentum transitions, by including 
the quantum defect appropriate to the experiments [7]. Clas¬ 
sical trajectory Monte Carlo simulations [17] were also in 
agreement with the experiments [7]. 

All theoretical efforts rely on the impact p.jameter for¬ 
malism, in which the projectile is a classical particle moving 
along a definite trajectory. The dipole interaction has been 
proven to be a good approximation for the projectile-target 
potential because of the long-range Coulomb interactions 
and the decisive role of large impact parameters. For slow 
moving ions. Stark mixing can occur without energy trans¬ 
fer. The dynamics of the Rydberg atom is therefore adia¬ 
batic. The orbit of the Rydberg electron can still be consid¬ 
ered elliptical, but its shape and orientation change slowly 
during the collision time, which is very much longer than the 
orbital time. This classical mechanics picnire translates into 
the quantal description by restricting the dynamics to the 
energy shell, as prescribed by adiabatic perturbation theory. 

In this paper, a unified theory for the general time- 
dependent solution of collisional Stark mixing is presented, 
both in the classical and quantal formulations. The excep¬ 
tional rich dynamic symmetry of the hydrogen atom provides 
the key foundation that enables both the classical and quantal 
solutions to be constructed in a unified way (Sec. II) by using 
group representation theory. This classical-quantal corre¬ 
spondence transcends the well known Ehrenfest’s theorem 
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(as observed in the general case of weak field-atom interac¬ 
tion [15.18]) just because of the SO(4) dynamical group 
symmetry* of the energy shell of the hydrogen atom. The 
agreement, as expected, is very good. It is shown that the 
present quantal solution (developed in Sec. HI A) can be for¬ 
mulated (Sec. Ill B) so as to provide the rotating coordinate 
frame formal result obtained in [12,15,16]. The efficiency of 
the present quantal solution is demonstrated in Sec. VI. A 
new classical solution applicable to transitions between arbi¬ 
trary angular momentum states is derived in Sec. IV A. We 
define the classical transition probability, in a language de¬ 
signed to exploit the dynamical symmetry, as the normalized 
volume of phase space accessible to both initial and final 
states in Sec. IV B. Monte Carlo simulations (Sec. V) arc 
also performed to yield results in agreement with the classi¬ 
cal expression for the Stark mixing probabilities. 



no. l. Geometry of the Stark mixing collision. 

mation, which is the same as the dipole approximation, the 
interaction potential V=er*£ is 


II. KINEMATICS OF STARK MIXING 

Stark mixing occurs when the electron of a Rydberg atom 
changes its angular momentum, without changing its energy, 
as a result of a collision, at large impact parameters, with a 
slow massive particle of charge Z\€. In addition to the en- 
ergy, given by the constant Hamiltonian 



the angular momentum L= rX p of the unperturbed Rydberg 
electron and the Runge-Lenz (or eccentricity) vector 


A= 


p^r-(pr)p~m^e 




/ 


V-2m,£, 


( 1 ) 


,R-r 

V(r.R)»-Z,e^ —= 


Zte^ 


/f3 vb dt 


Rr 


Zie^ d<I> 
ub dt 


{y sin<I>-Hrcos<I>). 


( 3 ) 


The impact parameter b, the impact angle 4>, and intemu- 
clear vector R arc displayed in Fig. 1. The angular momen¬ 
tum of relative motion -/ivb (where is the 

reduced mass of the projectile-target system) remains con¬ 
served since (so that and L are effectively decou¬ 

pled). 

Various frequencies or time scales are important to the 
present discussion of the collision, and are as follows. 

(i) The projectile rotation (collision) frequency 


which is directed toward the pericenter and normalized to 
angular momentum units, are ^so conserved. These quanti¬ 
ties define the dynamic SO(4) symmetry of the hydrogen 
atom with given energy E, which is a subgroup of the global 
SO(4,2) symmetry group. Because the SO(4) group is iso¬ 
morphic with the direct product SC)(3)©SO(3) of two rota¬ 
tion groups, a special decomposition 


d^ bv ^ ^ V I 

fix)m which the collision time Tcqu can be defined, 
(ii) The transition firequency 


L=M+N, A*M-N (2) 

permits the dynamics of the hydrogen atom to be separated 
into two decoupled motions. The generators M and N act 
independently as angular momenta and arc also conserved 
quantities for the unperturbed Rydbe^ atom. They evolve 
independently [9] with time on plication of an electric 
field. 

The orbital electron interacts with the time-dependent 
electric field S(t) generated by the passing projectile of 
charge Z\e. In the weak-field approximation, this field is 
constant over the spatial extent of the atom. In this approxi- 


h 


of the Rydberg electron. For transitions n-^n± 1 between 
neighboring levels, o)n,HZ\ simply » 

the orbital frequency (the Bohr correspondence principle). 
Here v^^vo/n are the averaged orbital radius 

and velocity. 

(iii) The Stark precession frequency 




3 Zia„v^ 

2 R^ 


‘This group contains sufficient generators to enable one to formu¬ 
late the dynamics of the system solely in terms of operations of 
irreducible representations of the group [19]. 


for the precession of A about R provides the prccessional 
frequency of the Runge-Lenz (eccentricity) vector A of the 
Rydberg orbit about the field direction £, 
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(iv) The spin-orbit coupling frequency a>so corresponds 
with the maximum fine-structure splitting and is approxi¬ 
mately [20] a\^oj^ln, where aps is the fine structure con¬ 
stant. 

(v) The quantum defect frequency cuqq is the precessional 
frequency of the electron orbit due to its interaction with the 
polarizable core. This frequency is important when one- 
electron atoms are considered other than hydrogen. The com¬ 
bined polarization of the core due to the orbital electron and 
the charged projectile has to be taken into account. Given the 
quantum defect 8 /. the quantum defect frequency is 

when / is sufficiendy large such that the core penetration 
and relativistic correcfions can be ignored [21]. 

By considering the txpiicjt) factor in time-dependent per¬ 
turbation theory, several types of collisions can be classified, 
as in [22], by comparing the above frequencies. 

The Stark mixing parameter a is defined as the following 
ratio between the Stark and collision frequencies: 

, 4 ) 

(On 2 \ bv I 2 m, 

When a< I, then (ost^ow*^ ^ collision time is much 

shorter than the Stark precessional time so that /-changing 
or Stark sudden transitions are favored. This is in contrast to 
Stark adiabatic transitions where a> 1; the electronic angu¬ 
lar momentum does not change since the atom has sufficient 
time to relax to the Stark effect 
The orbital parameter 13 is the ratio 


^Orbital Sudden ■ Ji Orbital Adiabatic ■ * 

-- ' 1 —1 - 

* . . . Stark Adiabatic ^ T Stark Sudden-* 

no. 2. Partitioning the impact parameter space 6. 

decreases, b^ increases outward and b^ increases inward, 
thereby limiting the extent of the two sudden regions where 
n changes and / mixing occurs. The variation with v can be 
represented by a phase-space diagram partitioned into 
the characteristic regions as illustrated in Fig. 3. For t;>u* 
=(3Zi/2)*S„, the (rt,/) changing and / changing (orbital 
and Stark sudden) shaded regions overlap and expand, in 
direct contrast to ultracold speeds where the orbital 

and Stark adiabatic (clear) regions increase and the shaded 
regions diminish and do not overlap, thereby indicating few 
collisional changes. The region of interest here is the overlap 
of the orbital adiabatic (w/?<w«) region, b>bc 
= {v/v„)a„, with the weak-field (a> 5 <a>^) region, b>b* 
= (3Zi/2)^^^a„,.i.c., the region in Fig. 3 defined by b>bc 
for u>c;*«(3Z|/2)*'^i;« and by b>b* for 

The following formulation assumes that the Rydberg 
atom, during the collision, occupies the same degenerate en¬ 
ergy shell. The main element of the perturbation potential (3) 
is the electron position vector r, which by Pauli’s replace¬ 
ment rule [23] 


r-(r> = 


2p. 


(5) 


^=a,*/a,. = -y 


b • 


is replaced by its average (r), a procedure valid for orbital 
adiabatic collisions (see Appendix A for a detailed explana- 


For an orbital adiabatic collision, when or b>bc^ the 
orbital electron adjusts itself adiabadcally to the slow ion 
perturbation. Since no energetic transitions oc¬ 

cur. The orbital sudden regime, 0>l or b<bc, is associ¬ 
ated with oi^Tcou^l and impulsive n—►n' transitions. 

The product of the a and $ parameters defines the ratio 
u}^lQ)^-{3Zxt2){anlb)^, which depends only on b (and not 
v). For weak fields the Stark mixing splitting 
^ha )„, the (n-^n± 1) energy gap. This also means that the 
intemuclear distance R is much greater than the mean orbital 
radius . In this approximation, the electron’s orbital time 
is then much shorter than any characteristic collisional time 
to cause / changes to the elliptical orbit The vectors A and 
angular momentum L, which are constant for the unper¬ 
turbed motion, then become good dynamical variables for 
the description of the perturbed motion, within the weak- 
field approximation. 

With respect to orbital motion, the collision is sudden or 
adiabatic according to b<bc and b>bc* respectively, 
where bc-(v/v„)a„. With respect to the Stark frequency, 
the collision is adiabatic or sudden according to b<bs and 
b>bs. respectively, where bs~{Vnlv)a„, The impact pa¬ 
rameter b space can then be partitioned as in Fig. 2. As u 



no. 3. Partitioning the u-b phase-space map into regions char¬ 
acterized mainly by (a) energy changes, (b) energy and angular 
momentum changes, (c) angular momentum changes, and (d) no 
changes. Regions for strong and weak field collisions are also 
shown. 
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tion). Since Zxe'vibv is then aA, the perturbing potential 
( 3 ) can be written in terms of the components A 2 and A 3 as 

Vi a) = e{r)‘E(t)= — sin 4>-HA 3 cosd>), 

under the adiabatic, dipole and classical path assumptions. 
‘Moreover, the components {Lj ,A 2 ,A 3 } generate a subgroup 
of the original symmetry group. The solution of the problem 
can then be written in terms of these symmetry-group gen¬ 
erators and the Stark parameter a, which acts as a coupling 
constant. Under the above approximations, the collision pa¬ 
rameters u, and Zj become combined into one parameter 
a. 

The cross section for^tark mixing is 



Departure from this variation is governed by 
l<P{b,v) and by the physical limits imposed upon the a 
integration. 

The discussion above is valid for both quantal and classi¬ 
cal descriptions of Stark collisional mixing, since both in¬ 
volve only kinematics and general dynamical symmetry ar¬ 
guments. It is shown below (in Secs. Ill A and IV A) that 
both quantal and classical dynamics are governed by the fol¬ 
lowing generic equation: 

SU d<t^ 

i— = t:a — {J 2 sin<t>- 1-73 cos<P)U (V 

dt dt ^ ^ 

for the time evolution operator U within the rotation Lie 
group SO(3). The generators {ij group have 

the commutators where €jtn is the Levi- 

Civita antisymmetric permutation symbol for any j,k 
= 1,2,3. The required solution of Eq. (7) is 

This can be easily verified with the aid of the relations 
= cos X-h 73 sin X., 


The probabilities obtained here are defined only in the 
full orbital adiabatic region b>bc for v^v* and in the adia¬ 
batic region restricted by the weak-field condition b>b* for 
(see Fig. 3). When and b<b^, the Stark pa¬ 

rameter a> 1. Since the transition probabilities are bounded 
for large a, the contribution to the a integration is vanish¬ 
ingly small for large a, decreasing as a“ , and can be ne¬ 
glected for a> 1. Cross section ( 6 ) can then be defined at 
ultralow speeds when the lower limit b„i^ to the 

^•integration 's taken as the weak-field limit b*. For higher 
speeds the probabilities determined here arc valid 

only in the full adiabatic region b>bc and do not hold in the 
(oribital sudden) region b*<b<bc required in ( 6 ). In prac¬ 
tical calculations of Eq. ( 6 ), various physical effects such as 
quantum defect, spin-orbit coupling, and Debye screening in 
a plasma determine an upper limit to the b integration and 
hence a lower limit to the a integration. For example, 
the spin-orbit splitting overlaps with the Stark splitting when 
where (3Zi/2)^flo* Similarly, the 

quantum defect comes into effect for the critical 

The Debye radius Ro 
= (ikr/4n^Zie^N)*'^, where Tand Af are the temperature and 
number density of the projectiles, is another viable upper 
limit to the impact parameter [4]. Stray electric fields in the 
collision region can also impose an upper limit to the impact 
parameter. The decision of which limit should be adopted in 
the definition ( 6 ) depends, of course, on the specific problem 
considered. These cutoff procedures are crucial for low an¬ 
gular momentum transfers where the transition probability 
P{a) cannot offset the 1/a^ singularity (cf. Ref. [4]) as b 
—®. The initial and final angles <to and ♦, between which 
the Stark mixing is effective, are also dependent on the spe¬ 
cific cutoff procedure [14]. For trajectories with zero deflec¬ 
tion (A 4 >=-' 7 r), then Eq. ( 6 ) varies universally as 


= cosk-y^ sink, 

=y 3 cos k ■Hy 2 sin k, 

which are derived from the basic identity 
k k^ 

jjlA.B]^ j^lA.lA^B]]-^ • • • (9) 

and the above commutation relations. The net polar angle 
A<I> swept during the collision between to and t is —4>o* 
The initial condition UUqJq)^! is automatically satisfied. 
If no cutoff radius is considered, as and 

4>—►tto (the classical deflection angle) as /-►». For th^ 
simplest case of distant straight line trajectories <I>o=0, and 
the evolution operator is then 

t/(oo,-oo) = exp[i7r(y|±ay3)]e‘*’"*'» 

or, in terms of finite angle rotations 7 ^y,n] b y angle (p 
about direction o and the parameter y- + is 

£/(ao,-oo) = 7^[-yir,( l/y,0,± a/y)]7^7r,( 1,0,0)]. 

in. QUANTAL THEORY 
A. Quantal Intrasbell dynamics 

The Schrodinger equation for the time evolution operator 
UitJo) is 

SU 

ih — ^{Ho^V)U. (10) 

where Hq is the free atom Hamiltonian and V is the interac¬ 
tion potential (3). If the projectile is moving sufficiently 
slowly, adiabatic perturbation theory can be a^jplied and then 
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the whole dynamics of the target atom becomes restricted to 
the initial degenerate energy shell (the orbital adiabatic re¬ 
gion in Fig. 3). This simple fact has two major consequences. 
First, the position operator and hence the perturbation poten¬ 
tial (3) commute with the unperturbed Hamiltonian, as one 
can prove directly from the matrix elements of the commu¬ 
tator [r,//o] between any states within the energy £„ shell. 
The potential in the interaction representation 

is then identical with the ^Ajtential in the Schrodinger repre¬ 
sentation (V/= VO, and the equation to be solved, in the in¬ 
teraction representation, is 

dVi 

( 11 ) 

- at 


where £//(r,ro) = exp(///o//^)£/(r,ro)«xp(-///o/o/^)- Second, 
the components x,y,z of the position operator do not com¬ 
mute between themselves when restricted to the energy shell. 
This follows from the well known Pauli “replacement*’ [23] 
r—*--3rtA/2 (see Appendix A for details). This shows that 
the position vector r behaves, within the intrashell dynamics, 
like an angular momentum and is denoted by {r} when op¬ 
erating only within the n-shell. In fact, the set of operators 
{Lx.-2[y}pJ'h-2{z}pJZ) generates a rotation group and 
Eqs. (11) and (7) are identical when 
-2{y};?„/3 and ^ 3 = “2{z}p„/3. Using the solution ( 8 ) of 
Eq. (7), the exact solution of Eq. (10) of Sec. n for the 
evolution operator, within the adiabatic approximation, is 
then 


U/(r,/o) = ^ 




exw - 




( 12 ) 


This can also be directly verified by substituting Eq. (12) in 
Eq. (10) and using Eq. (9) with the appropriate commuutor 
algebra. 

The transition probability for a general i—►/ transition at 
time t is 

where 'P, is the target wave function, which tends in the 
asymptotic limits (»-»£«) to the unpeiturbed basis set 
<t>y(r.r) = (^y(r)exp(-j£yt/ft). The transition amplitude for a 
Stark mixing process is 

afa^{nP\Ui{^-^)\na). (13) 

where the initial state i-lna) at f= evolves to the final 
states/"Iny?) at r = »; a and ^ now label the states within 
the same energy shell. The superscript n will be omitted, 
since all dynamics is restricted to the energy shell described 
by quantum number n. When a and P label states with a 
given magnitude and projection of the angular momentum, 
the transition amplitude is given by Eq. (12) in Eq. (13) and 
is feasible and efficient for practical numerical applications. 


The core of solution (12) is the exponential of the opera¬ 
tor Li-2ap«z/3. By using basic commutator algebra and 
Eq. (9), this operator is diagonalized as 


where ^ = 2p;, arctan(a)/3 and y=vl"^Q:“. The solution 
( 12 ) has therefore the alternative form 

(14) 

which illustrates very effectively how the action of the slow 
distant encounter charged projectile coming from the nega¬ 
tive z axis is decomposed into successive rotations about the 
X axis and alternating impulsive momentum transfers ('zq) 
along the y axis. 

It is interesting to compare the solution (14) obtained in 
this orbital adiabatic limit with the purely impulsive solution 
presented in sec. Ill C, The evolution operator (14), for un- 
defiected collisions, yields 

In the limit that a—^O. ^“► 2 ap „/3 = Aq/2 and 

exp(iiry£iM)“*>exp(i'rrLi/fi). Then 


C//(®.-») 




[- 0^1 

1! ( 2A 2! I 2A j 


Ih 




which reduces, with the aid of exp(«iTL|/A)y'cxp(-/ir£,i/ft) 
=y’(cos 7 r)"=(-y)*, to 

+ (15) 

i 

This limit merges with the impulsive result Eq, (20) below. 


B. Formal development 

It is interesting, however, to note that by introducing Pau¬ 
li’s replacement directly in the potential (3) and by writing 
the Rungc-Lcnz vector as A=M-“N. the potential decom¬ 
poses as 


where 


= - a(M 2 sin 4> + M 3 cos <l>)<t> 


and 


= a(Ni sin <l> + A ^3 cos <I> )<t>. 

Because the commutators [Mf,A^J = 0, [M,,/fo] = 0. and 
= 0 (for any i,;~ 1 , 2,3 combination), the problem 
00 ) becomes separable, exactly in the same way as the clas¬ 
sical Stark mixing equations become decoupled (see the next 
section). The time evolution operator then factorizes as 
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TABLE I. The four bases useful for describing the quantal states of the hydrogen atom. 


Basis 

Quanrum 

numbers 

Complete set of 
commuting observable 

Origin 

orbital 

\n/m}o 

Hq. l\ Lj 

Standard for spherical coordinates; describes 
correctly the states of the field-free atom 

parabolic 

jn 

Hi. Hj. Lj 

Separation of Hamiltonian // = 1 + //j in 
parabolic coordinates, f=r+z. 77 = r-:, 
tan^>=v/jc; n = + 1 

Stark 

\nqm)s 

Ho. A 3 . £3 

Parabolic basis; describes the Stark states for 
small electric field when the interaction 
—eSz is diagonal; 

algebraic 


Ho. W3.H3 

The two rotation groups in which the dynamic 
symmetry group SO(4)"SO(3)©SO(3) 
decomposes using Eq. (2); the equivalent 
angular momentum for both SO(3) 
representations is; = (n-l)/ 2 ; fi-{m^q)f2 and 
v=(m —^)/2 


U^Uh^UuUn. ( 16 ) 

where, of course. and Uu and 

are the solutions of the equations ihdUM/St=Vi^Uu and 
ihdUsIdt-Vf^UM. respectively. Using the group- 
theoretical result Eq. (8) of Sec. II, the solutions for the 
operators U and are then 

= <,•<«>«.'* expt - «M(4> - <I>o)(A^i - 

(17) 


and 

'*exp[ - «t»o)( Afi + '*■ 

(18) 

In c^culating the amplitude (13), four interesting basis 
sets Oe chosen for the one elr-c^on hydrogenlike atom. 
Table I summarizes key properties of these bases. The orbital 
basis is useful for describing the field-free atom, before and 
after the collision, whereas the algebraic basis appears natu¬ 
rally as a basis where A/ 3 and ^3 are diagonal. The solution 
(16) has the simplest expression in this algebraic basis. Ail 
four bases in Table I span the degenerate energy shell and 
can be equally adopted to characterize the hydrogen atom. 
The algebraic basis spans a tensorial product of two spaces 
(1 m>® I v>) corresponding with spaces used for a matrix rep¬ 
resentation of the product SO{3)0SO(3). The two spaces 
have the same dimension because = = 

-(n^-l)fiV4 and arc associated with two angular mo¬ 
menta with ;■ = («“ l)/2. 

The transition amplitude between the two algebraic states 
is then the product of two amplitudes for Af and N indepen¬ 
dent actions, exactly in the same way in which the classical 
Af and N vectors evolve independently in time, i.e., 
= element of 

a ; = (rt-l)/2 dimensional represenution of the rotations 
represented by Eq. (16). For example, from Eq. (17), one 
gcis = with 


F(a) = I>^^ 


^yir. 





where H<^,(n| ,n 2 ,n 3 )] is the Wigner matrix represenution 
for the rotation 7y^,n] by angle about direction n (see 
[24] for the explicit expression). The transition probability in 
the space of N is the element v' v of the matrix F( - a). 

Calculation of the transition probability between orbital 
states requires the explicit unitary transformation between 
the orbital and algebraic bases. This can be obtained by di¬ 
rect scalar products of orbital and parabolic states for which 
explicit coordinate representations are known. The result can 
be written in terms of hypergeometric functions [25]. How¬ 
ever, an equivalent result is provided by the SO(4) 
-SO(3)©SC)(3) isomorphism. The orbital state, as a com¬ 
bination of two angular momentum states, is (see [26] for 
example) 


|n/m)= 2 C^^”^\nfiv), 

where the transformation matrix is given by the stan¬ 
dard Clebsch-Gordan coefficients {jfijv\/m). The transi- 
tion amplitude for the m* transition is then 




= 2 




,(-a). 


(19) 


which can be expressed in matrix form as 

a/-«../«=Ti{C<'""'’F(-a)C‘^-’^F^(a)]. 

where is the transpose of matrix C. The above result Eq. 
(19) is in exact agreement with the solution obtained in 
[12,15,16] using the rotating frame approach. The quantal 
development here in the fixed frame is exquisite in that it 
follows exactly the same reasoning basic to the exact classi- 
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cal mechanics solution (Sec. IV A). This result exhibits the 
essential power of the SO(4) symmetry group for the energy 
shell of the hydrogen atom. The common SO(4) symmetry 
therefore transcends the chosen formulation (classical or 
quantal) and provides a classical-quantal correspondence at a 
level more fundamental than Ehrenfest*s theorem and the 
Heisenberg correspondence. In practice, the fourfold summa¬ 
tion (19) and the use of the Wigner rotation matrices V for 
F(a) are not very efficient and the difficulty of calculation 
increases dramatically with n. Instead, the solution (12^ .pro¬ 
vides a simpler approach, since the matrix elements of the 
argument in the exponential have simple expressions directly 
in the orbital basis. The array of transitions is obtained at 
once, within one matrix exponentiation of a band diagonal 
matrix for which efficient algorithms are available [27]. 

When the projection of the initial and final angular mo¬ 
mentum is not determined, the transition probability is 

1 

The exact quantal solution can therefore be derived di¬ 
rectly without making use of unnecessary complications of a 
rotating frame and a fictitious magnetic field. The structure 
of the present solution Eq. (14) represents a sequence of 
alternating momentum transfers (in the y direction) and rota¬ 
tions about the x axis. This recognition motivates the follow¬ 
ing section. 


= ( 20 ) 

the inelastic form factor amplitude. This agrees with other 
alternative derivations [28.29] for the probability of an im¬ 
pulsive transition. The momentum transferred in an (impul¬ 
sive) rectilinear collision is 

Aq=(2Zie2/6t;)y = (4aA/3nao)v. 

The transition amplitude, valid in the orbital and Stark sud¬ 
den region (cf. Fig. 3). is therefore 

a}7( a) = (exp - i j “ 1 ), 

which connects with the a—*0 limit Eq. (15) of the adiabatic 
result Eq. (14). Even though this approximation is appropri¬ 
ate for the amplitude for transitions n/between dif¬ 
ferent energy shells, it is a good approximation, in particular 
for very sm^l a [sec Eq. (15)]. if the normalized 



transition probability is adopted for the problem of intrashell 
transitions. This normalization is a consequence of the dif¬ 
ference between the operators {r}, appropriate only to in¬ 
trashell transitions, and r for all transitions. 


C. Impulsive limi t 

In the other extreme situation (the orbital suddcn-Siark 
sudden region in Fig. 3, where a^l), the impulsive limit, 
the collision is very fast and the collision time is much 
smaller than the orbital time. The potential has again he 
same form in both Schrodinger and interaction repr menta¬ 
tions. V= V/. since, in the equality (/]V|i)*(/|Vfl*). V*', 
the exponent [with ai^,-(E/-E,)/A] can be re¬ 

placed by zero. Now the position operator has r the normal 
behavior, in that its components commute between them¬ 
selves. The impulsive transition amplitude is 

a}7=(<^/r)|exp- 


IV. CLASSICAL THEORY 

A. Classical iotrashcll dynamics 

The angular momentum vector L and the Runge-Lenz 
vector defined by 


A = p;‘ 


are constant for the unperturbed classical Rydberg atom. 
Moreover. A L=0 and + = = In the pres¬ 

ence of an electric field of intensity £, the angular momen¬ 
tum L changes at the rate 


— 

It 


erxi. 


as may be directly verified upon using the closure relation. 
Since the force F= - VfV(r,r) acting on the Rydberg elec¬ 
tron is impulsive and imparts momentum Aq, then 


On assuming that the collision is orbital adiabatic (((>< 0 )^ 
and i is constant over one period), the slow change of L 
during the collision is the classical average 


exp- 



V(r./')^/' = exp- 



F(r,/)d/ 


AL 

T 


l—\ 

\dtl^ Tjt-T/2 


(rxt)dt- “e{r)xi(t) 


= exp-^Aqr 

so that the probability amplitude for an impulsive collision is 
simply 


over one orbital period 7. Since the weak-field approxima¬ 
tion holds, the vectors L and A then change 

very little over one orbital period. Using Pauli’s replacement 
rule (r)** -3 A/2p„. the following set of coupled equations 
can then be deduced [9,22]: 
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dA. . dL_ * M(/) = t?M(f''o)M(fo) and a similar equation holds for the 

— - - (i> 5 Rx L, — W 5 RX A. evolution of N. The first equation in Eq. (22) becomes 


where both <t>s= and R vary with time. Under the sub¬ 
stitution 


dO.j 

d<t> 


-ae,„tRnC/tj. 


L+ A 
M , 


N=' 


L-A 


( 21 ) 


the above set of differential equations becomes decoupled to 
yield 

— = —oicRxM, -^ = -t-w5RxN, (22) 
dt dt 

where the magnitudes^Af^=Af^=(L^+A^)/ 4 =n^fiV 4 re¬ 
main constant throughout the collision. The classical analysis 
for constant electric fields is given by Bom [9]. For time- 
independent ws . both M and N process with constant fre¬ 
quency (Os about the (fixed) direction of intemuclear axis R. 
For general time-varying ws. the system of differential 
equations (22) does not have an exact solution. Percival and 
Richards [22] have used classical perturbation theory to 
solve Eqs. (22) and then provided a diffusional theory of 
angular momentum mixing. Bellomo et aL [17] approached 
the same problem by proceeding via the time evolution 
propagator for M and N in the rotating frame, an 

approach that results in formulas too complicated for physi¬ 
cal changes AL and AA to be extracted. A special solution 
for transitions from angular momenta state ^=0 has been 
recently obtained by Kazansky and Ostrovsky [13-15]. 

An exact analytical solution is, however, possible under 
the weak field and orbital adiabatic approximations for a 
classical projectile trajectory and when the magnitude of vec¬ 
tors M and N remains constant during the collision, as for 
the present case of intrashell transitvjos. These vectors are 
then obtained at any moment by orthogonal transformations 
from the initial values MUo) and N(io). Let these transfor¬ 
mations be DuU.h) and respectively. Then 


where e,;* is the fully antisymmetric permutation symbol. 
Because the infinitesimal generators 7* for the rotation group 
arw matrices with elements Uk)ij- ~ ^ijt [30] for ex- 
ami/le) and. since i^ = ( 0 .sin<l>,cos<t>). the above equation in 
matrix form is 

dUu , , - 

= —o(sin<I>72 + COS<t>73)t/M . 
av 

where [7* Jj] - . This equation is the matrix represen¬ 

tation of the group equation (7) and the solution is the 
matrix representation of the general solution ( 8 ). The final 
vector M' is therefore obtained from the initial vector M by 
three successive rotations, 

M'*7?[-<I>,( 1.0.0)]72[yA«*»,( 1/y.O,-a/y)] 

x72[«I>o.(1.0.0)]M. (23) 

The solution for N has the similar form 

N'« 7l[-<*'.( 1.0.0)]7e[ rA<l».( l/r.0,+a/y)} 

X72[«|)o.(1.0.0)]N (24) 

obtained simply by replacing a by - a in the corresponding 
equation for M. The matrix 7Z{*1>,b] is the rotation matrix 
for a vector aitd corresponds with the representation of the 
abstret rotation (specified by the angle and the direction 
D ^rotation) on the three-dimensional vector space. Analyti¬ 
cal expression for and Un can be obtain^ [31] as ex¬ 
plicit functions of a and A*!*. In particular, when <t=0 and 

<|>o=ir, the (3X3) matrix Dm « 


J 




y"*[l+a^cos( 7 ry)] 
ay"' sin(ir y) 
ay"^[cos(iry)- 1 ] 


ay"' sin(fl-y) 
- cos( IT y) 
y"'sin(iry) 


oy"^[l-cos(iry)] " 
-y"'sin(iry) 
-y"^[a*+cos(iry)]. 


(25) 


and On is obtained from 0^ by replacing a with - a. 

As a result of the collision, the initial state of the target atom, specified by 
(L'.A') according to 


the vecton (L.A). changes to the final state 


, . Om-On ^ 

. - 2 2 

(26) 

. Om-On, O^+On ^ 

‘“ 2 2 

(27) 


For the undeflected u-ajectory of the projectile, when 4>=0 and <l>o- rr, explicit results are: 
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L\ = y‘‘[l -fa-cosCiT y)]Li + Qrr“‘ sin(tr 7)^3-^ -cos(7r 7)1^3, 

^2= -cos(7r y)L2- r"** sin(7r y)L3-r ay’* sin('n- y)Aj. 

1,3 = y" * sin( TT y)L2“ y"'[a“-^cos( tt y)]L3 + cry’“'[cos( tt y) - I ]A,, (28) 

= Y”‘[ 1 cos( ?r y)]A i + ay"* sin( tt y)i-2'*'<^y'‘[ I “COs( tt y)]L3. 

A2= -cos( 7 r y)A2- y“* sin( tt y)A3 + ay“* sin{ 7 r y)Li. 

A3= y“* sin(' 7 r y)A2“ y'“[a* + cos(Tr y)]A3-h ay“^[cos(tr y)- l]Lj. 


Here y= vl + < 2 ^ andAhe components of the initial and final 
vectors are defined in the fixed coordinate frame of Fig. 1. 
Similar expressions have also been obtained [31] for general 
A<1>. The above exact solutions are easily verified and satisfy 
the invariant relations 

L' A'=L A=0 
and 

The orbit of the final state (n.L') is confined to a plane 
perpendicular to the final L' and the energy is preserv^ (n 
is not changed). 


I 

the collision. Following the definition (29). the overlap vol¬ 
ume of accessible (L, A) space that contains both initial and 
final states is 

J j S(\L\-/h)S{\V\-/'h)S(\\\-h^n--/'^) 
xS{L-X)dLdA. (30) 

The transition probability is then, in a geometric sense, the 
ratio of two volumes: the volume of the accessible 

states compatible with the required final angular momentum 
and the volume of the acceptable initial states V„/. The / 
transition probability is therefore defined as the ratio 


B. Classical transition probability 

The initial state is defined by the angular momentum L 
and Runge-Lcnz A vectors. Apart from the constraints that 
(i) the magnitude of the L vector is /A, (ii) the magmtude of 
the A vector is h in the given state, and (iii) L and 

A are always orthogonal, the two vectors arc completely ran¬ 
dom in the six-dimensional space {L}®{A}, which is a map¬ 
ping of the usual (r,p) phase space. The initial angular mo¬ 
mentum can have any value between 0 and nh. The special 
case of zero initial angular momentum requires a separate 
analysis, presented at the end of this section. In the following 
discussion, the initial angular momentum is assumed strictly 
positive. 

The hypcrsurfacc in the {L}®{A} space on which the 
initial state is uniformly distributed is restricted by the above 
constraints and has the volume 

V„^=| f S{\L\-h^S{\A\-hyl^i^)S{L A)dLd\, 

^ (29) 

which, upon integration, reduces to 


Each point within this manifold evolves during the collision 
according to the rules (28), so that only a fraction of possible 
initial states can have the final angular momentum after 


P 


(fi) _ ^n//‘ 


( 31 ) 


of phase-space volumes. Transformation to the alternative 
set of vectors M f^id N defined by Eq. (2) facilitates evalu¬ 
ation of the integ^l v30). The Jacobian of this transformation 
is dL dA= 8 dM dN. With the aid of the identities 

d( L • A) = <5( Af 2 - ^2) = <5( A/- A# )/2A/, 


5( I L| - //i) = <5( V2M^( 1 + cos/. MN) - /h) 

<S(|A|-AVn^-/^)=<5(>/2M*(l-cos/.MN)-AVn^-^) 

= S{M-nhf2)fi^Jn^^ldM, 

S(\L'\-/'h)^SU2MHl+cos/LM'ti')-/'h) 

= <S(cos/.M'N'I ))/'*/A^^ 


the accessible phase-space volume is 
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V,// 


■=// 


m^N^dClf^dasdMd!^ 


SiN-M) 

^ 2M ^ AM 
X S COSZ.MN— 

X (Jcos^M'N' - 


/h 

S(M-nhl 2 )— 




-1 


2M^ 

2M^ 


— 


--1 


This finally reduces to the simpler form 



FIG. 4. The proper rotation angle Y “ * function of the Stark 
parameter a for the net polar angle swept A<t> = -rr. -3n-/4. 
- ir/2. and - ir/4. 


fdn„d{if, 

X (5(cosZ.MN— ^)i(cos/.M'N' — 

where the integral is now over only the angular pan of the 
vectors M and N and where the parameters ^ and are 
simply related to the initial and final states by 


the angle x increases up to tt as a increases 
to unity. For smaller values of A<I>. the angle x increases 
monotonically with a. 

The classical transition probability (31) is then 


/3= 



2 /'^ 



(32) 


The final vectors M' = and N' = are given by 
finite rotation [Eqs. (23) and (24)] of the initial M and N. so 
that the relative angles Z.MN and Z.M'N' are independent 
of the specific coordinate frame chosen. Then 

M'-N' (f7„M)-(d7«,N) 
cos4M'N' =--- -5 - 


(f7j^£/MM)-N 


= cosZ.M''N, 


where M'' is obtained by rotation from M using the operator 
U]fOM . Being a product of two rotations, this operator is 
also a rotation about some direction v by the angle x deter¬ 
mined from the trace 

Tii0lDM]^l + 2cmx 

of the rotation operator. The proper rotation angle x depends 
only on the collision (Stark) parameter a and the polar angle 
^(t>s(]>—(pg swept out during collision time interval {IqJ) 
and is independent of the initial or the final state of the Ke¬ 
pler atom. It is determined by 


cosj = [ 1 + cos( Vl + a*A4>)]/( 1 + a^)- 
For small c^l, cosy— 

Y-4asin(A<I>/2) + 0(a^). The plot Fig. 4 of the “univer¬ 
sal” function y(®.^*^) (or the case A<I>= ~ ir shows that x 
has a maximum at a—0.9 and is never greater thra it. When 


X <5(cosZ.MN-/S)<S(cosZM''N-^'). 

where the angle between M and M" is y- The dCl/^ integral 
can be done first if one chooses to work in spherical coordi¬ 
nates with the z axis along the vector v. In doing this, the 
vectors M and M* are fixed and have the coordinates ( 6 . 0 ) 
and ( 6 ,y) as depicted in Fig. 5. The surface area element is 
dClu-dicos 9)dd>, where S and ^ are the spherical coordi¬ 
nates of the vector N. Instead of the {d,d>) system, a new set 
of coordinates can be defined by («i ,M 2 ). the angles of N 
with M and, respectively. Vi*. The surface area element is 
m w anv=dMidM 2 /sin A. -.ler.. A is the angle between the 
N^'? and NM" arcs, as in Fig. 5 A proof of this result is 
derived in Appendix B. 

The dfl/y integral is now simpler to evaluate and yields 



FIG. 5. The geometry and coordinates used in solving the inte¬ 
gral (31). 
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2 /' , 

-r d(cos0)(sin A sinuj sinw 2 )“ . 

where cosWj =2/^/n^“ 1, cosu 2 ^ 2 /^'"/n^^ I* and the fac¬ 
tor 2 arises because the {d, <(>) transformation is 

not single-valued. 

Basic trigonometry (see [32]) applied to spherical triangle 
MNM" yields 

cosF = cosui cosu 2 + sinw| sin cos A 

so that 

sin u 1 sin tt 2 cos A = (cos r - - cos F)(cos r - cos r >), 

where F. = Uj±W 2 are the limits to F as A rotates through 
27r. Since 


where 




cos(U|’^u:)~cos;^f 
I ~ cos 


cos(ui ~ w-*)-cos y 

5(//n,/7n;a) =-;-^-. (33) 

l-cos;^ 


On denoting cos 0 by r. the transition p*^Laoility is the 
one-dimensional integral 


o(") 2/^ i f dz 

irAn* I-cosatJ yj{z--AHB- z') 


cos F = cos 0^ + sin 0^ cos x 
for spherical triangle Mi'M", then 

sin u I sin U2 sin A = (I - cos y)[(cos^0 - A) 
x(5-cos^0)]*'2, 


where the limits of integration are defined by the 
condition of reality for the square-root function. The 
last integral can now be expressed in terms of the 
complete elliptic integral ^(m) = /j'^(l - m sin^x)**^dx so 
that 



if 

S<0 



if 

B>0. 

A<Q 

(34) 

if 

B>0. 

A>Q 



r 


provides the exact* classical probability as a function of /, 
and sin^y2, which combines the Static parameter a 
= (3Z,/2) {a„v„/bv) and the net polar angle swept during 
the collision A<t into one function. The probability (34) sat¬ 
isfies detailed balance 2/P/^/=2/'P//> , where 2/ is the 
classical weight of the state n /. 

Inspecting the definitions (33) reveals more qualiutive as¬ 
pects. First, >1 and 5 are always less than 1, and only B can 
attain I when U| = M 2 * i-®** collisions /=/'. The 

transition probability reduces, for this case, to the simpler 
form 




1 

Aayjn^-/^ 


in the limit of small a, thereby exhibiting the 1/cr singularity 
of elastic transitions. For transitions with and small 

enough a, the factor B is negative and then the transition 
probability is zero. B is always greater than A because 
cos(u,-u 2 )>cos(tti+U 2 ). In the limit of each / and /'-►O 
or rt, then Wi.2=0 or ^ and A-B<1. For /—tran¬ 


sitions, the probability is zero. For the importan' :asc of zero 
initial angular momentum, the transition pruoaoility P in the 
limit of /-♦O is 




/'Khn'^) 


sin(;r/2)x/sin^;r/2)-(/’/n) 


T- 


(35) 


This result was also obtained in [14], When /' =n, the tran¬ 
sition probability is 


^ sin(Af/2)V(‘i^/«)*”Sin^(Ar/2) 

When the argument of the square-root function is negative, 
the transition probability is zero. Similarly, for maximum 
initial angular momentum /'=n, the transition probability 
has the limit 




sin(;^/2) •v'{/'/n)^-sin-‘(Af/2) ’ 
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a=0.2 Qr=0.4 



FIG. 6: Contour plots of the solutions of A 
= 0 (solid line) and 5=0 (dotted line) for various 
values of a. In gray zones 5<0 and the transi¬ 
tion is classically forbidden. 


where* again* the transition probability is zero unless the 
final angular momentum is large enough to make the 
argument of the square-root function positive. If 5>0 and A 
is small, the transition probability has a singularity* typical 
for classical mechanics. Because the complete elliptic K 
function diverges logarithmically when the argument is very 


close to unity, the transition probability has the following 
expansion for small A\ 

P^4=(2ln2-iln>4/5). 

yjB 


I/n = 0.071 //n-0.36 





FIG. 7. Contour plots of the solutions of A 
= 0 (solid line) and 5=0 (dotted line) for various 
values of /. In gray zones 5<0 and the transi¬ 
tion is classically forbidden. 
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When B has small but positive values, the transition prob¬ 
ability has a .finite limit because K(0) = Tr/2. Of course, 
when B approaches zero from negative values, the transition 
probability is zero. Thus, the singularities in the transition 
probability are given by the solutions of the equations A 
= 0 and 5=0. which are ^ 2 = -“i-Y* 

A map of the various zones in the plane of reduced initial 
and final angular momenta (//n./Vn) is displayed in Fig. 6 
for the four values a = 0.2, 0.4, 0.6. and 0.8 of the Stark 
parameter. In the central region, A is negative and B is posi¬ 
tive. Within the lower left and upper right comers, both B 
and A are positive. The transition is classically forbidden in 
the upper left and lower right comers (gray zones) where B 
<0. Along the solid line, .4=0 lines, which represent 
cos’In) + cos" = ;t/2,7r/2 ^ transition 

probability has a logarithmic (cusp) singularity. Across the 
dotted 5=0 lines, which represent cos’\/ln) 
-cos"‘(/Vn) = ;t'^2,7r-;^/2, the transition probability 
jumps from zero (in the gray zone) to some finite value (in 
the central zone). As a—^0, the two inaccessible regions 
(where 5<0) increase until the central region with B>0 and 
A<0 becomes an elongated line strip lying along the diag¬ 
onal /=/'. Only elastic collisions are therefore permitted 
in the limit a—►O. As a increases to unity, the classically 
forbidden zones diminish and the collision becomes more 
and more effective in its ability to induce larger angular mo¬ 
mentum changes. 

Figure 7 presents corresponding maps to Fig. 6. The same 
characteristic regions are now displayed in the plane of final 
reduced angular momentum /7n and the Staric parameter a 
for four values of the initial reduced angular momentum 
//n = 0.071, 0.36, 0.64, and 0.93. Again, the classicaUy for¬ 
bidden regions (gray zones) correspond to the condition B 
<0 in the left upper and lower comers. The elastic /' = / 
transitions are always possible, even when a—^0. Again, 
along the solid (>4=0) and dotted (5=0) lines, the transi¬ 
tion probabilities have cusp and step singularities. When a 
increases, the span of the possible final angular momentum, 
for given angular momentum, increases. Large angular mo¬ 
mentum transfer is only possible for collisions with large 
Stark parameter cr-^ I. Both Figs. 6 and 7 are key to the 
interpretation of the variation of the probabilities F///(a) 
with both and a, respectively. 

The result (35) obtained from the /-♦O limit of the gen¬ 
eral result (34) can be also proven directly. Because L=0, 
the classical orbits arc characterized only by the Rungc-Lenz 
vector A. In this case, the volume occupied by the initial 
state is 


Ko= j ^(|Al-rtA)dA=4irn^/l^. 

The volume of the accessible final states /' is 


V„o/'= I «5(|Al-nyi)5jj (Dn-Os)^ 


-/'h]dA, 


where the orthogonal matrices depend only on the 
Stark parameter a and are defined by Eq. (25). In a spherical 



no. 8. The geometry of the initial state described by the direc¬ 
tions of the angular momentum L and Runge-Lenz A vectors. 

coordinate system where the z axis is along the principal 
direction vof the rotation matrix 0 ^ 0^^4 (as discussed earlier 
in this section), the magnitude is 

j I Af/^A I = n sin( ;^/2) sin , 


where 6 is the angle between (or £7\A) and the direc¬ 
tion R The transition probability is then 



4S(n^sin(>f/2)sin d—/*h)sin Odd - 


_ /7(An^) _ 

sin(;t^) Vsin^(Y/2)”(/7/j)^ 

in agreement with the /-^O limit (35) to the general classi¬ 
cal transition probability (34). 

V. CLASSICAL DYNAMICS SIMULATION 

The present classical Monte Carlo simulations are differ¬ 
ent from the standard classical trajectory Monte Carlo simu¬ 
lations in that the initial state is not specified by the position 
and the velocity of the orbiting electron. Instead, the simula¬ 
tion begins with a random distribution of the initial states in 
the {L}®{A} space, propagates each state according to the 
rule (28), and then performs the statistical analysis of the 
final distribution of states to provide, in the limit of an infi¬ 
nite number of trials, the probability P/*/ for a given tran¬ 
sition. This section describes the correct way to generate the 
initial distribution of states. 

The angular momentum L and the Rungc-Lenz A vectors 
are sufficient to completely specify all characteristics of the 
atom’s orbit, provided they are orthogonal, as for the case of 
pure Coulomb attraction. Only five components are therefore 
independent and two of them, the magnitudes L and A 
= characterize the shape and size of the orbit 

There are then three angles that specify the orientation of the 
orbit in space, for a given assignment of the energy and 
angular momentum. The direction of the angular momentum 
vector L is arbitrary. Two random numbers, the projection 
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L. = cos^i and its azimuthal angle <Pi, generate a uniform 
distribution of L on the unit sphere. The direction of the 
Runge-Lenz vector A is uniformly distributed in the plane 
perpendicular to £, and the angle ip between A and a given 
fixed direction iV in this plane is the third random number 
required to simulate an initial state. In summary, three ran¬ 
dom numbers are required: (i) L; = cos Si uniformly distrib¬ 
uted wiLhin the [- LI] interval, (ii) tpi uniformly distributed 
within the [0,2Tr] interval, and (iii) tp uniformly distributed 
within the [0,27r] interval. The initial L and A vectors arc 
then 


I sin cos <pi ^ I sin ^2^05 (pi 

L-A/l sin^)Sin<pj L . sin^ 2 sin92 

\ cos^i / \ COSS2 
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/ 
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no. 9. Matrix representation of - aA3 for n = 3. 


The spherical polar angles Sj and 92 of A, as illustrated in 
Fig. 8, must now be expressed in terms of the above random 
variables ^1, 91, and 9. Basic trigonometry [32] applied to 
the spherical triangles ZJVA and LzA yields 

cos ^2 = sin 9 sin Si 
and 


ping is given by floor( and m = i — 

The index k counts the degeneracy of the energy shell, and 
runs from I to n^. The matrix element 

+ J +m-f 1) 1 


0 = cos Si cos S 2 + sin sin S 2 cos( 92” i)» 

respectively. On solving these two equations for S 2 and 92* 
the arbitrary Runge-Lenz vector 

A(5i,9it9) 

( -cos Si cos9i sin 9-sin 9i cos 9 
-cos Si sin9j sin 9+cos9i cos 9 
sin sin 9 

is then expressed in term of the random variables. It auto¬ 
matically obeys the constraint requirements A-L=0, A^ 
+ = and A-yV=cos9- 


VI. NUMERICAL EXAMPLES 

In this section, numerical examples for calculation of the 
transition probability between states with given angular mo¬ 
mentum (n/—n/') arc presented. There are three main 
methods used in this paper. The classical Monte Carlo simu¬ 
lation, as described in the preceding section, requires the 
running of a large number of trials to sample the three di¬ 
mensional space of arbitrary parameters. The explicit classi¬ 
cal mechanics expression (34) is used directly. 

The quanta! c^culation is based on Eqs. (13) and (12). A 
matrix representation for the operator L i — aA 3 [where Pau- 
li’s replacement (5) was adopted] is required. Instead of the 
spherical basis |/m>, which is difficult to use in this case, 
we define a new linear basis obtained by mapping the (/,m) 
quantum numbers to a unique index in 

such a way, for example, that (0,0)—► I, (I,"" l)—^2, (1,0) 
(2,-2)—5, and so on. The inverse map- 


of Li is nonzero only for A/=0 and Am = ± 1, which re¬ 
flects the fact that the cylindrical symmetry of the Rydberg 
atom is broken by the precession of L about the field of the 
projectile. These m-changing transitions arc, however, con¬ 
ditioned by the full structure of solution (12), which shows 
that such transitions are only in evidence for nonzero a. The 
matrix element 


'j- 




S /'/-1 


(2/+l)(2/-l) 

(2/+3)(2/+l) 




of the component Aj— — {2J3n)z along the fixed Z axis of 
quantization is nonzero for — 1 and A/n=0 transitions. 
These dipole transitions only contribute for nonzero a. The 
matrix L, - aAs has then a band diagonal struemre. as illus¬ 
trated in Fig. 9 for the special case of n = 3. Explicit analyti¬ 
cal formulas for P^"’/(a,A4>) can be direedy obtained [33] 
for small n = 2,3. 

The transition amplitude for transition k-»k' is the kk' 
matrix element of the exponential of the matrix -/A«I>(L| 
-aAj). sandwiched between the rotations implied by Eq. 
(12). When a—0, the dipole forbidden transitions are not 
possible, because the transition matrix -exptrrZ,|M still 
maintains a band diagonal structure. As a increases, more 
and more off-diagonal elements become populated, leading 
to dipole forbidden transitions. Efficient algorithms, using 
Fade approximations, are available [27] for matrix exponen¬ 
tiation. The full array of transition probabilities is then ob¬ 
tained all at once. 
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no. 10. The Monte Carlo simulation (sicplike lines), the clas¬ 
sical (solid line), and quanta! (dots) transition probabilities 
for given a=0.2 and inhi^ /=4 within the /i = 20 shell. 


The examples, provided by Figs. 10-13, demonstrate that 
the Monte Carlo simulation yields results identical to the 
classical dynamics expression (34). The transition probabili¬ 
ties for given Stark parameters (a=0.2 and 0.6) and initial 
angular momentum (/ =4, 12, 14, and 18, respectively) are 
represented as a function of the final angular momentum ^ . 
This also provides the distribution over the final angular mo¬ 
mentum states, which result from collisions, at given a, from 
an initial population of states with the same initial angular 
momentum /. In this example, n = 20. The present corre¬ 
sponding quantal results are also represented in the same 
graphs by dots. Vertical dotted lines indicate the positions of 
the singularities, corresponding to the ^=0 and S=0 lines 
in Fig. 6. The >1= 0 and 5= 0 singularities produce cusp and 
step variations in / as increases through the singular¬ 
ity (Fig. 6). Four distinct and characteristic classes of varia¬ 
tion of with then emerge. These are displayed in 
Figs. 10-13, where the predicted (cusp,step), (step,step), 
(cusp,cusp), and (step -asp) classical variations are exhib¬ 
ited. These results are fvlly representative and can be ana¬ 
lyzed by vertical cuts through Fig. 6 appropriate to a given 
/. The / and a parameters in Figs. 10-13 correspond to the 
values //n = 0.2 and 0.6 and to //n = 0.7 and 0.9 in the first 



RG. 11. The Monte Carlo simulation (steplike lines), the clas¬ 
sical (solid line), and quantal (dots) transition probabilities P/>Ao) 
for given or = 0.2 and initial ^*12 within the n * 20 shell. 0 for 
all /V 
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RG. 12. The Monte Carlo simulation (steplike lines), the clas¬ 
sical (solid line), and quantal (dots) transition probabilities P/ - /< a) 
for given a=0.6 and initial /= 14 within the n = 20 shell. 0 for 
all 

and third plots of Fig. 6, respectively. The steps indicate a 
classical threshold /I or a classical cutoff for transitions 
to a final value of The cusp-step variation of Fig. 10 
indicates that transitions to are classically allowed 

and the step-step variation of Fig. 11 indicates that classi¬ 
cally accessible are within the range ♦ The 

cusp-cusp variation of Fig. 12 is associated with the fact that 
transitions to all are classically accessible (cf. Fig. 6, plot 
3 ), while the step-cusp variation of Fig. 13 signifies that only 
transitions to are classically possible. For a given 

initial angular momentum / and Stark parameter a, the po¬ 
sition of the cusp singularities, given by the solutions of the 
equation -4=0, is 



Expressions (36) are also solutions of the equation B=0 for 
the step singularities. The threshold value is a cusp (so¬ 
lution of A= 0) provided /<n sin and is a step (solution 
of S=0) provided />n sin;t/2. Similarly, the cutoff is 



RG. 13. The Monte Carlo simulation (sieplike lines), the clas¬ 
sical (solid line), and quantal (dots) transition probabilities P/fA 
for given 0=0.6 and initial /= 18 within the ii = 20 shell. Across 
the first dotted line B changes sign and A is negative on both sides 
of this line. A changes sign across the second line while B remains 
positive. 
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FIG. 14. The classical and quantal Stark mixing transition probabilities within the n*28 energy shell, for a*0.2, as a function 

of the final angular momentum and for various initial angular momenta /. 


a cusp if />n cos or a step if /<n cos "Th® separa¬ 
tion between the singularities is 



when /In'^sinxf2, and is 



when //n^sinx^- The maximum separation of sinx is at¬ 
tained at /♦/n = sin;t/2. This occurs in Fig. 6 where the 
A-0 and 5=0 curves both intersect Hat /In axis at /*ln. 
The transitions have significant quantal probabilities 
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RG. 15. Density plots of the transition prob¬ 
abilities. c^culated within the quantal treatment, 
for four values of a = 0.2. 0.4, 0.6, and 0.8 in the 
(/,/') plane. The probability increases as the 
gray becomes darker. 


only when is in the classical accessible region and exhibit 
the characteristic exponential decreasing behavior in the 
classical inaccessible regions. 

The second set of examples compares between the classi¬ 
cal and quantal transition probabilities within the /i = 28 en¬ 
ergy shell. A set of graphs, for various / and fixed Stark 
parameter a = 0.2, is presented in Rg. 14. The transition 
probability is plotted versus the final angular mn^ icntum. 
The cusp-step, step-step, and step-cusp variations ar^ ill ap¬ 
parent, together with the superimposed quanta! oscillatory 
behavior within the classical accessible region. The oscilla¬ 
tions possibly originate from the fact that there arc two 
points N of intersection of the two arcs (drawn on the sphere 
of Fig. 5) with centers M and M*" with separation r(a) and 
radii «,(/) and Each point of intersection provides 

equal classical contributions to P/»/(o) while the different 
phases associated with each point pnxluce the quantal (semi- 
classical) interference oscillations, exhibited in both the 
variations of P/>/{a) with / and a fixed, as in Fig. 14, and 
the a variation of P/»/{ot) in Fig. 16 for / and fixed. 
Figure 14 also illustrates that the separation (/+-/l)/n 
between the various discontinuities increases to sin at at 
/In - sin ;^/2 and then decreases as / is increased from 0 to 
/I - 1, in accord with Eq. (36). In general. Fig. 14 shows that 
the classical picture is complementary to the quantal in that it 
has the ability to explain the general overall behavior of the 
quantal results and to provide the general firameworic on 
which the quantal results rest It has also provided the vari¬ 
ous regions (between cusps, steps, etc.) that remain obscured 
within the quantal treatment 

The set of plots in Fig. 15 is the quantal correspondent of 


the similar classical set of plots presented in Fig. 6. Density 
maps in the (/,/') plane are shown for the same four values 
of Stark parameter a. The first map for a=0.2 corresponds 
to the results of Fig. 14, where cuts along various values of 
/In are made. The quantal transition probability increases as 
the gray areas becomes darker. The same zones outlined in 
Fig. 6 can be recognized and the boundaries between them 
are in exact correspondence with the classical equations A 
= 0 and S= 0, as discussed at the end of Sec. V. The quantal 
transition probabilities are practically zero over the classical 
forbidden regions, occupying the upper left and lower right 
comers of Fig. 6. The quantal probabilities are maximum on 
the ridge given by the equation ^=0, where the classical 
dynamics results in singularities. 

When a is very small, decreasing toward zero, only elas¬ 
tic transitions are allowed so that only the principal diagonal 
is exhibited, as in Figs. 6 and 15. The quantal calculation 
yields unity for the probability of elastic /=/' transitions in 
the a—^0 limit and the classical result diverges as 1/a in the 
same limit. All other transitions have zero probability. This 
feature is responsible for the well known [4] divergence of 
the cross section (6) for elastic transitions. 

Figure 16 displays the probability for transitions originat¬ 
ing from the initial level /=5 to various final levels, as a 
function of the Stark parameter a. Again, the agreement is 
expectedly very good. For small a, there is always a classi¬ 
cal inaccessible region (S<0) for quasielastic /* i^/^ ex¬ 
cept for the fully elastic transitions /' = /. A threshold step 
at a= aT- is therefore displayed for the probability of transi¬ 
tions with /' = /. This property is fully explained with the 
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aid of the plots in Fig. 7 of /7/i versus a for various values 
of //n. As is increased, determined by the intersec¬ 
tion of /7n with the 5=0 curve increases and a step (5 
= 0) to cusp (>1=0) variation with a is obtained, as exhib¬ 
ited in Fig. 16. For higher /7n values, intersection with the 
>1=0 curve (and therefore the cusp) disappears. Figure 16 
(for the specific n = 28 case considered, or in general Fig. 7) 
also shows that low ► high and high /—► low 


transitions are precluded (5<0) except at high values of a 

— 1 . 

Figure 17 shows that the quantal probabilities are high in 
the classical accessible regions (S>0) of Fig. 7 and are 
more significant in the A<0 region than in the ^>0 region. 
Variation of the transition probability along the horizontal 
line for /' = const provides plots as exhibited in Fig. 16. 
Figure 17 is the quanta! correspondence of Fig. 7. 
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a ct 


no. 17. Density plots of the transition probabilities, calculated 
within the quantal treatment, for four values of the initial angular 
momentum /’=2. 10, 18, and 26. in the (or,/') plane. The prob¬ 
ability increases as the gray becomes darker. 

Figure 18 exhibits variation of the integral factor 


r(« 

/'/ 



p(«) 

a ira: 



which appears in the cross section (6) as a function of the 
final angular momentum for various initial angular mo¬ 
menta ?. It is assumed here that A<l>=-ir. This integral 
does not depend on the projectile properties (velocity, impact 
parameter, or charge) but depends only on the initial and 
final state of the target. Due to the 1/a^ singularity, the cross 
sections for elastic (/=/') and near elastic transitions are 
very much enhanced. Various cutoff (a,nio) procedures can 



no. 18. Classical and quantal comparison of the integral I/*/ 
defining the Stark mixing cross section (6) within the n»28 energy 
shell, for various initial angular momenta / as a function of the 
final angular momentum 


be introduced from physical considerations specific to the 
actual problem discussed. For the results presented in Fig. 
18, an arbitrary cutoff a^n=0.01 was introduced for both 
the classical and quanta! calculations. With the exception of 
elastic transitions, where special care must be exercised since 
the present model fails, the agreement between quantal and 
classical model is much better than one would expect from 
the traditional correspondence principles such as Ehrenfest’s 
theorem and Bob"\> correspondence. The exceptionally nch 
symmetry group SO(4), characteristic for one-electron at¬ 
oms in both the quantal and classical treatments, explains the 
quality of the agreement obtained. 

VU. CONCLUSIONS 

In conclusion, we have presented a case study of the long¬ 
standing problem of the transition array in atomic 

hydrogen produced by distant collisions with slow heavy 
charged particles. Complete formulations from the quantal, 
classical, and Monte Carlo simulation viewpoints have been 
developed. The SO(4) dynamical symmetry of H(n) has 
been exploited to provide exact classical and quantal solu¬ 
tions (28) and (12), under the adiabatic, dipole and classical 
path assumptions. A classical expression (31) for the transi¬ 
tion probability is presented in a language that exploits 
the dynamical symmetry and is used to obtain exact analyti¬ 
cal results (34). A new exact (fixed-frame) representation Eq. 
(12) of the quantal solution has also been presented and di¬ 
rectly appli^ to provide exact quantal probabilities for the 
full array of transitions. This solution is feasible for efficient 
numerical calculation of probabilities for transitions involv¬ 
ing even large n and all / and , in contrast to the previous 
quantal (rotating-frame) version [15,16], which was applied 
only to the transitions. Exact analytical expres¬ 

sions for P/»/ also be obtained [33] for low n. 

The exact results for P^/^ oscillate about the clas¬ 

sical background. By revealing essential characteristics that 
remain obscured within the quantal treatment, the classical 
results complement the quantal results. Further development 
would include a semiclassical analysis capable of reproduc¬ 
ing the oscillatory structure in the quantal P^/^ without sac¬ 
rificing the physical transparency of the classical model. 

Although the numericid results presented here are associ¬ 
ated with a straight line projectile path (i,c., A<1>= - tt), the 
theory is suitable for all polar deflection angles A<l> swept 
out during the collision by the classical trajectory. Specific 
problems (nonhydrogenic atoms, stray fields, the Debye ra¬ 
dius in plasmas, etc.) may impose further restrictions on the 
range of impact parameter h. which results in restrictions on 
a. The basic theory presented can still be applied to these 
cases with minor adjustments. The classical trajectory appro¬ 
priate to a rotating dipole can be determined to account for 
the influence of the target on the projectile’s motion. 

In summary, both classical and quantal solutions of Stark 
mixing have been presented in a compact form reflecting the 
mathematical beauty of the problem as well as their prag¬ 
matic value. It is probably one of the last remaining prob- 
lerns in collision physics capable of an exact solution. 
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Note added in proof. A general concise analytical expres¬ 
sion for the quantal probabilities for all n has since 
been obtained [34] in a form (a) which is easy to use for 
numerical evaluation, even at very high n~100. (b) which 
yields compact analytical results' for n = 2-5 and (c) which 
naturally provides the classical limit (34) obtained here. 
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APPENDIX A: PAULI’S REPLACEMENT 

Pauli’s replacemeiu [23] for the classical one-electron 
atom follows from the straightforward calculation of the av¬ 
eraged projection of the position vector on the Runge-Lenz 
vector, 

. 3 3 1 


ducing M=(L+A)/2 and M=(L-A)/2. Each M and N 
operators generate separately a so(3) subalgebra, such that 
SO(4)»SO(3)eSO(3). 

Following [26], it is useful to prove the following. 
Result, The matrix elements of the operators H defined by 


are 2icro between any states wiih the same energy, i.e., ft 
»0 within the energy shell. 

Proof. The commutator oetween the Pauli-Lenz vector 
and the position vector can be written as 

(AD 

Because the commuuuor [W.X] is zero for any Hermitian 
operator X. when restricted to the energy shell {the hyper- 
virial theorem [35]), commutator (Al) is simply 


where e is the eccentricity of the orbit corresponding with 
the energy E. Because the perpendicular component of r on 
A averages to zero, the following rule is valid, on averaging 
over an orbital period: 


3 A 


r** — 


2Pn 


provided A docs not change significantly in this time. Here 
Ph- “ >/ 2m,£ is the characteristic orbital momen¬ 
tum. 

The classical Runge-Lenz vector has the symmetrized 
(Pauli-Lenz) quantal form 


A= 


1 

-(pXL-LXp)-m,e^r 


/ 


alternative to the form given by Eq. (1). The operator A has 
the following properties: 


[A^,//]=0 (a conserved quantity), 
[Lj,At]=iEejb.A^ (also a vector). 


•''*] “ 2p^ 

in the Hilbert subspace of degenerate states with the same 
energy. Thus, the vector ft commutes with the Pauli-Lenz 
vector and has the commuutois 


[Af^.n^] —[Ny,flt]= 

with the vectors M and N obtained by decomposition (2). 
With the aid of the basic Jacobi identity for commuutors. 

[A.[B.C]]+[B.[C,Ar+[C.[A.B]]=0. 

it follows directly that 0=0. For example. 02=0 ^ 

= A/ 3 , B=Ni, and C=fl 3 . 

This result provides a concise proof that Pauli's replace¬ 
ment (5) is valid whenever the dynamics is constrained to the 
constant energy manifold (cf.. [18] for a lengthier proof). 

APPENDIX B: AREA ELEMENT FOR SPHERICAL 
BIFOCAL COORDINATES 


[A; .A*] = IA (its components do not commute). 

A L=L-A=0 (vector A is orthogonal on L), 

+ DA^ (constant for intrashell transitions). 

These commutation relations define the SO(4) dynamic 
symmetry group for the restricted motion of the orbital elec¬ 
tron to the energy shell. The Hamiltonian W is an invariant 
and can be used to label matrix representation.^ for this 
group. The correct energy levels for the hydrogerr atoms re¬ 
sult from the symmetiy without solving any differential 
equation. The SO(4) operators can be disentangled by intro¬ 


The spherical bifocal coordinates are given by the angles 
u, and M2 between a point (P) and two fixed foci (A| and 
A 2 ) on the unit sphere. Of course, these coordinates are 
unique only on the half sphere. 

Theorem. The area element in spherical bifocal coordi¬ 
nates is 

</u, dU 2 
dS*^ —r-r—, 
sm A 

where A is the angle between the arcs joining a given point 
on the sphere with the two foci: 
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A = Z.A,PA 2 . 

Proof, Here are two proofs of the theorem. One is short 
and intuitive. The other one is longer but a bit more rigorous. 

Short Proof The length of the arc described by the point 
P on the sphere for an infinitesimal change 
when u, is kept fixed, is dsi. Basic spherical triangle geom¬ 
etry provides 

du2 
sin A 

and a similar relation 


where 6-/LPD and <t>-CAiD, Applying again spherical 
irigonomccry to the spherical triangles A^PD and AjPD. 
then 

P := (cos u (, V cos* Z. PD - COS' u sin/.PD), 

where 

sin/.PD=sinui sinu 2 sin A/sin r, 

where, by definition, A = Z.A;PA 2 , and T 

= /.A|A 2 . The goal is to express all coordinates of P in 
terms of M| and «2 and to do the corresponding derivatives in 
Eq. (Bl). Spherical trignometry yields 


^ 5^ = sin^ u 1 sin^ M 2 A 

sin A =[C 0 S(M| “M2)~COSr][cOS r-COS(M| +M2)]. 


for the infinitesimal arc dji- The elementary surface area is On defining 
given by the cross product of the two arcs dSiXds^. Since 
the angle between these arcs is, in the first order of approxi¬ 
mations, again A, then 


5i = sin^[r±(Mi + M2)], 


dS^ 


dM| dui 
sin A ' 


as advertised. 

Direct Proof A direct proof calculates by brute force the 
area element using 


J-.sin-[r±(M,-u2)]. 
the angle A is then determined from 

S yjrjlsjl 

;- =2 —:-^- 

Sin Ml sinM2 sin Ml sin M2 


dS^ —X-— dMi^Mj. (BI) 

duy dU2 

The length of the calculation depends, of course, on the par¬ 
ticular coordinate frames of choice. For example, one can 
choose Aj and A 2 in the equatorial xOy plane, with Ai on 
the Ox axis. Let D be a point on the arc A 1 A 2 such that the 
arc PD is perpendicular to A 1 A 2 (like a meridian passing 
through P). The coordinates of the point P are then 

P:=(cos dcos<f>, cos ^ sin <t>, sin 8), 


and the coordinates of F are now 

s^S^-s.S- 


‘ sin r 

sinF 

Then Eq. (Bt) finally reduces to 


sinui sinu. 

dui dui 

^5= --- duidui- 

Sind ‘ 

S J 


which proves directly the above theorem for spherical bifo¬ 
cal coordinates. 
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